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PREFACE 

There  doesn’t  seem  to  be  any  officially  agreed  upon  definition  of  precision 
engineering  but  for  purposes  of  this  monograph  the  following  statement  could  serve. 

Precision  Engineering'.  The  application  of  higher  technology  and  science  to  the 
design  and  manufacture  of  components  or  subsystems  calling  for  geometric  and 
position  control  to  within  fractions  of  a  micron. 

An  effective  graphical  definition  of  the  objectives  of  precision  engineering  in  a 
specific  discipline  is  by  Fig.  P.l  from  the  work  of  Professor  Niro  Taniguchi  of  the 
Tokyo  Science  University.  While  this  illustration  focuses  on  machinery  accuracy  it  is 
nevertheless  an  excellent  guideline  to  the  comprehension  of  the  implicit  goals  of 
precision  engineering  in  general.  Currently  in  the  U.S.,  there  is  a  belated  but  growing 
interest  in  precision  engineering.  The  primary  motivation  for  this  change  is  the 
growth  of  high  technology  both  for  defense  and  commercial  needs.  Such  products  are 
characterized  by  miniaturization  and  very  close  tolerances.  There  is  another 
compelling  force  which  should  activate  an  increased  interest  in  precision  engineering. 
It  is  the  loss  of  markets  for  the  kind  of  manufactured  hard  goods  which  might  be 
roughly  classified  as  the  top  curve  of  Fig.  P.l.  Second  and  even  third  world  nations 
are  now  capable  of  performing  at  this  industrial  level  but  at  less  cost  to  our 
disadvantage.  Product  lines  or  methods  of  manufacturing  which  exploit  more 
sophisticated  engineering  such  as  precision  engineering  could  develop  a  more 
competitive  edge  for  us.  It  is  no  coincidence  that  nations  such  as  Switzerland, 
Germany  and  Japan  each  have  a  balance  of  trade  which  we  envy.  The  common 
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denominator  among  the  three  nations  is  a  history  of  training  and  education  in  PINE 


MECHANICS  or  its  current  identity  PRECISION  ENGINEERING.  Several  Japanese 


universities  now  have  departments  of  precision  engineering.  They  have  recently  coined 


the  expression  "mechatronics"  giving  further  evidence  of  ferment  in  activities 


associated  with  precision  engineering.  There  is  even  a  Japan  Society  of  Precision 


Engineering.  Clearly  the  advantage  to  Japan  is  growing  and  might  become 


unreconcilable  for  the  U.S.  if  the  current  level  of  effort  here  is  not  stepped  up. 
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Fig.  P.l 


Educational  programs  in  the  ways  of  precision  engineering  aimed  at  not  only 
student  engineers  but  also  practicing  engineers  transferring  to  this  field  is  very  likely 
the  quickest  avenue  for  the  catch  up  program.  A  strategy  for  initializing  this  activity 


is  the  development  of  a  series  of  short  texts  or  monographs  for  beginning  graduate 


students  and  engineers  in  practice  desiring  retraining.  These  texts  would  also  keep 
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faculty,  with  a  little  formal  background  in  precision  engineering  design,  to  transfer 
their  skills  to  this  movement.  The  conventional  comprehensive  design  text  authored 
by  one  or  two  generalists  is  probably  unsuited  to  the  program  contemplated. 
Monographs  each  concentrating  on  a  topic,  within  the  field  and  composed  by  an 
authority  in  that  field  will  serve  better. 

Fig.  P.2  illustrates  a  possible  matrix  of  topics  on  which  a  series  of  monographs 
could  be  based.  This  report  is  a  possible  model  for  the  very  first  item  of  the  matrix. 
As  such  it  is  an  exploratory  experiment  in  the  suggested  program. 


Fig.  P.2  -  Topics  &  Contents  Matrix  for  Precision 
Engineering  Design  Monographs 
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1.1 


I.  INTRODUCTION 

A  SENSE  OF  PRECISION 

The  engineering  design  of  precision  machinery  and  devices  in  this  monograph 
dwells  on  the  practice  and  the  analytical  tools  necessary  for  the  development  of 
mechanical  systems  whose  performances  are  defined  by  their  ability  to  control  and 
measure  dimensional  phenomena  in  orders  of  less  than  two  microns  (2//m).  It  might  be 
surmised  that  the  intuitive  comprehension  of  minute  dimensions  or  displacements  by 
the  man  in  the  street  falls  off  rather  completely  at  resolutions  below  1mm.  For  most 
engineers  the  limits  of  intuitive  comprehension  of  smallness  falls  off  at  about  0.1mm. 
In  this  series  the  focus  is  on  minuteness  less  than  one  micron  (.001mm).  The 
"Dimensional  Spectrum  of  Phenomena",  Fig.  1.1  exhibits  the  typical  dimensions, 
tolerances  etc.  of  various  items  of  engineering  and  scientific  phenomena.  Graphical 
juxtapositionings  such  as  this  help  the  practitioner  appreciate  the  magnitudes  of 
accomplishment  needed  to  achieve  various  levels  of  precision.  While  proficiency  in 
engineering  analysis  and  experimentation  are  paramount  to  the  creation  of  successful 
new  designs  and  their  development,  one  can  not  overlook  the  power  of  subliminal 
information  and  patterns  in  the  engineering  of  precision  systems.  More  often  than  we 
realize  problem  definition,  which  precedes  analytics,  is  initiated  by  a  subconscious 
ordering  of  past  observations  and  associations  of  clues.  This  subliminal  process  is 
variously  identified  with  such  terms  as  "practical  experience",  "insight",  "creativity", 
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"detectivity ",  etc..  As  nothing  comes  of  nothing,  insights  and  problem  definition  must 
come  from  accumulations  of  information,  analogs,  comparisons,  contrasts,  etc..  To 
this  end  the  practitioner  benefits  from  familiarization  with  such  spectra  as  that  of  Fig. 
1. 1.  He  might  well  develop  the  practice  of  adding  new  data  as  they  become  available. 
A  cultivated  sense  of  the  micron  enables  the  diagnostic  or  design  engineer  "feel"  what 
is  taking  place  in  a  precision  system  even  before  the  hardware  is  produced. 
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Figure  1.1  Dimensional  Spectrum  of  Phenomena 


A  MODELING  LANGUAGE  &  SYMBOLS 


Much  of  the  analytical  prediction  of  machine  component  and  systems  behavior  is 
based  on  such  engineering  sciences  as  heat  transfer,  fluid  mechanics,  elasticity, 
dynamics,  electrical  and  magnetic  field  theory.  Each  of  these  in  its  evolution 
developed  its  own  languages  and  uniquely  named  laws,  principles  and  property  states. 
The  specialist  in  any  one  of  the  fields  by  virtue  of  single  minded  activity  is  quite 
comfortable  with  specific  names.  On  the  other  hand  project,  design  and  the 
administrative  engineers  must  perform  across  lines  of  disciplines  and  need  a  more 
universal  language  to  facilitate  understanding  and  accurate  communication. 
Fortunately,  despite  differing  vocabularies,  analogues  can  be  established  among  the 
engineering  sciences  so  that  a  unifying  language  and  a  mathematical  symbology  can  be 
generated.  The  analytics  of  these  monographs  will  be  presented  in  one  such  language 
the  principles  of  which  are  herewith  described. 

Generic  Variables 

Common  to  the  engineering  sciences,  cited  above,  is  the  concept  of  a  Through 
Variable.  This  is  usually  a  basic  extensive  quantity,  such  as  mass,  energy,  charge,  etc. 
which  is  being  transmitted  along  a  path  from  one  point  or  state  to  another  point  or 
state.  In  this  monograph  the  generic  symbol  for  a  Through  Variable  is  q„  H. 

Definition:  Let  q„  ,  be  a  vector  or  scalar  quantity  defining  the  transmission  or  flow 

rate  of  a  Through  Variable  from  point  or  state  a  towards  point  or 
state  J.  Because  the  transmission  is  from  the  prior  subscript  towards 
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the  posterior  subscript  the  sequence  of  subscripts  is  critical.  It  follows 
that: 

Qn../ 

The  Through  Variable  above  is  driven  from  a  towards  (i  by  an  Across  Variable 
typified  by  such  intensive  quantities  as  temperature,  pressure,  voltage  etc..  The 
generic  Across  Variable  will  be  symbolized  as  e„/,,. 


Definition:  Let  ely/t)  be  a  vector  or  scalar  value  of  a  generic  Across  Variable  at  cr 

less  the  value  at  //.  It  follows  that: 

®,y/o  ///• 


Some  specific  examples  of  Through  and  Across  Variables  are  illustrated  by  TABLE 


I. 


TABLE  I 


Engineering  Science 

Through  Variable,  q 

Across  Variable,  e 

Heat  Transfer 

Heat  Flow;  Q 

Temperature  0 

Fluid  Mechanics 

Mass  Flow;  M 

Pressure  p 

Electric  Fields 

Current;  I 

Volts  v 

Mechanics 

Force;  F 

Displacement  s 

5 
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Constitutive  Law 

Through  and  Across  Variables  are  related  by  Constitutive  Laws  perhaps 
the  most  familiar  example  is  Ohm’s  Law.  The  generic  form  of  these  laws  is; 

=  en/li  EI.l 

Proportionality  term  K„_y  is  the  average  property  of  the  media  or  element 
connecting  points  <*  and  ft  hence  these  subscripts  apply  to  K  as  well  as  to  q  and  e.  In 
most  engineering  problems  it  may  be  assumed  that  proportionately  K  is  approximately 
passive  and  constant  such  as  an  electrical  conductance  or  a  spring  constant.  In  such 
cases  the  order  of  ex  and  ft  subscripture  of  K  is  inconsequential.  Recall  that  on  the 
other  hand  the  order  of  the  subscripts  when  written  as  a.ft  and  (x/ft  is  critical  to 
polarity.  Attention  is  drawn  to  the  fact  that  while  <►  and  ft  appear  in  all  three 
quantities  of  EI.l  there  are  distinctly  different  separating  punctuations,  (.)  (/)  and  (-), 
for  the  three  distinctly  different  quantities  shown.  More  on  this  a  little  later  but  at 
this  point  it  can  be  stated  that  because  of  the  definitions  given  above  there  is  the 
happy  consequence  that  if  the  prior  subscripts  and  posterior  subscripts  appear  in  the 
same  order  on  both  sides  of  the  equation,  polarity  for  the  equation  is  automatically 
correct.  TABLE  II  illustrates  some  specific  constitutive  laws. 


'*V'vv.4: 


TABLE  II 


Engineering  Science 

Constitutive  Law 

9n  j  e„/j — 

Heat  transfer 

Fourier’s 

Q n.J  = 

Fluid  Mechanics 

D’Arcy’s 

M„..J  ~  K„_j  P 

Electric  Fields 

Ohm’s 

l  .1  =  K„_j  E„yj 

Elasticity 

Hook’s 

./  —  K,. -j_S^ 

Conservation  and  Compatibility  Principles 

Finally  two  other  principles,  Conservation  and  Compatibility  are  needed  to  model 
the  problem  into  differential  equations.  The  two  principles  invoke  the  effect  of 
geometry  of  related  points  on  the  derived  equations. 

Definition:  Conservation  Principle  is  a  summation  law  for  all  Through  Variables 

transmitted  from  surrounding  points  towards  a  common  vertex.  The 
generic  form  of  this  summation  in  terms  of  the  subscripts  defined 
above  is: 

L'  <fi  j  =  j  +  <i2  j  +  <ik  j  +  ' ' '  i  j  +  qn j  =  0  EI-2 

k  =  l 

Correct  polarity  is  preserved  if  the  points  1  through  n  surrounding  j  appear  as  the 
prior  subscripts  while  the  vertex  j  appears  as  the  posterior  subscript  in  all  terms. 
Kirchoff’s  Current  law  is  a  familiar  example  of  a  constitutive  law. 

Definition:  Compatibility  Principle  is  a  summation  law  for  Across  Variables 

relating  all  points  connected  in  a  chain.  Kirchoff’s  Voltage  law  is  a 


familiar  example. 
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The  generic  form  of  this  summation  is: 

e„/n  =  e„/i  +  e,/2  +  e2/3  +  *  •  ’  +  en/n  =  0 


El. 3a 


e"/b  —  e«/i  +  ei/2  +  e2/3  +  ’  *  '  +  en/;y  EI.3b 

In  the  form  of  I.3b  the  compatibility  principle  serves  to  generate  a  new  value  en/,y 

from  the  known  values  en/lt  etc.  in  a  connected  chain.  TABLE  III  illustrates  some 

specific  examples  and  names  of  these  two  principles. 

TABLE  III 


Engineering  Science 

Conservation  Principle 

Compatibility  Principle 

Heat  transfer 

Conservation  of  energy 

No  name 

Fluid  mechanics 

Conservation  of  momentum 

No  name 

Electric  fields 

KirchofFs  Current  Law 

KirchofFs  Voltage  Law 

Dynamics 

Conservation  of  momentum 

Polygon  of  Velocities 

The  algorithm  of  subscripts  in  EI.2  is  entirely  different  from  that  of  EI.3.  Because 
of  this  it  is  imperative  that,  while  an  adjunct  pair  of  Through  and  Across  Variables 
have  the  same  subscripts,  the  Through  quantities  have  the  dot  subscript  separators  (.) 
while  the  Across  quantities  have  the  slash  subscript  separators  (/).  Equation  EI.2 
having  the  confluent  algorithm  of  Through  Variables  must  contain  only  terms  with  a 
(.)  separator.  Equation  EI.3  on  the  other  hand  deals  only  with  Across  Variables.  It 
has  a  zig-zag  algorithm  to  its  subscripts  and  must  contain  only  terms  with  a  (/) 
separator. 


i\V 

V*' 
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Implicit  Subscript  (*) 

In  addition  to  the  physically  apparent  <*,.</  reference  to  nodes,  terminals  or 
points  in  the  field  of  a  machine  component,  the  language  of  modeling  requires  the 
recognition  of  a  particular  point  or  state  not  explicitly  apparent.  That  is  an  absolute 
reference  point  henceforth  identified  by  {*).  Electrical  sciences  make  explicit  use  of 
such  an  absolute  datum  point,  it  being  "ground"  and  identified  by  the  graphic  symbol  ( 
).  Curiously  the  other  engineering  sciences  fail  to  provide  a  symbol  for  an  absolute 
datum  even  though  the  concept  prevades  the  language  of  the  sciences  by  references  to 
absolute  temperatures,  coordinates  etc.  By  establishing  an  implicit  terminal  (*)  to 
represent  an  absolute  datum,  ambiguity  is  diminished  and  consistency  in  notation  is 
served.  The  lamentable  omission  of  a  graphical  symbol  for  a  ground  or  absolute  datum 
in  mechanical  sciences  makes  itself  felt  even  in  the  most  rudimentary  models  such  as 
the  following  observation. 


The  Absolute  Point  of  Reference  * 


There  are  three  fundamental  structural  properties  to  every  machine  component. 
These  are  flexibility  (which  implies  an  ability  to  store  potential  energy),  damping  (the 
agent  of  energy  dissipation)  and  mass  the  basis  of  kinetic  energy.  A  flexible  member 
can  be  symbolized  by  Fig.  1. 2a  that  is  a  spring  with  ends  1,  2.  There  is  a 
corresponding  constitutive  equation  with  the  ends  1,  2  playing  the  roles  that  <>,/? 
demonstrated  in  the  introductory  discussions  and  El. I.  The  damping  role  of  an 
element  can  be  symbolized  as  in  Fig.  I.2l»  by  an 


Fig.  1.2  Symbolic  Representation  of  Potential  Dissipative  and  Kinetic  Energy  Elements. 
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elementary  dash  pot  with  terminals  3,  4  and  a  constitutive  equation  where  force  where 
force  and  velocity  effects  are  in  equilibrium.  In  both  cases  there  is  a  consistency  of 
subscripts  based  on  reference  points  on  both  sides  of  the  equation.  Polarity  is 
rigorously  self  defined  by  the  order  of  these  subscripts.  Consider  now  the  case  of  a 
simple  mass  as  an  inertial  element.  Fig.  1. 2c  is  the  conventional  graphical 
representation  and  the  familiar  constitutive  equation.  There  is  only  a  single  reference 
point  on  this  conventional  representation,  that  is  point  5,  which  might  be  a  point  on 
the  surface  of  the  mass  at  which  an  external  mechanical  force  acts.  The  related 
constitutive  equation  would  show  only  that  one  subscript  which  not  only  fails  to  show 
th<  double  subscript  consistency  with  the  other  two  types  of  structural  elements  bti t 
even  makes  the  nature  of  the  acceleration  ambiguous.  Fig.  1. 2d  is  the  suggested 
symbolic  notation.  An  implicit  connection  of  the  mass  to  (*)  is  indicated.  This 
immediately  suggests  that  acceleration  be  rewritten  as  S5/<  and  so  establish  the 
acceleration  as  an  absolute  quantity.  As  such  the  symbol  (*)  evokes  the  Newtonion 
inertial  frame  of  reference  as  a  rigorous  analysis  should.  Taking  a  cue  from  the  first 
two  examples  of  elements  it  follows  that  mass  would  be  identified  by  two  subscripts 
i.e.  M5_*.  Subscript  *  now  reallirins  that  the  property,  mass,  of  the  envelope  5  is 
indeed  an  absolute  quantity  invariant  throughout  the  universe.  In  this  way  the  r.h.s. 
of  (he  constitutive  equation  is  established.  The  l.h.s.  is  made  to  assume  the  two 
subscripts  in  the  sequence  established  on  the  r.h.s.  by  the  acceleration  term.  Even  this 
"forcing"  of  subscripts  to  conform  results  in  a  useful  interpretation.  Consider  the 


inverse  force  F •  5.  By  the  basic  definitions  and  the  role  of  subscripts  in  a  through 
variable,  F .  r,  would  be  interpreted  as  the  force  transmitted  from  the  reference  point 


11 


1.11 


w, 


' 


;$l 

■S$ 


'sVi 

Si 


Vi1 


$ 


-v 

•  s 


(*)  to  the  body  (5).  This  interpretation  squares  up  perfectly  with  the  intent  of  the 
vaguely  defined  term  "force  of  inertia".  Finally  this  format  of  subscripts  reaffirms  the 
equivalent  roles  of  a  Newtonion  equation  and  a  corresponding  D’Almbert  equation. 
Consider  four  bodies  1,  2,  3,  4  exerting  forces  on  body  5.  Newton’s  equation  becomes 


Fi  5  +  F25  +  F3.5  —  m5_*  S5/*  =  F5  .  =  — F*  5 


This  can  be  rearranged  as 


Fi.s  +  F2.5  +  F3  5  +  F4  5  +  F*  5  =  0 


This  equation  is  D’Almbert’s  form  of  the  dynamic  problem  wherein  the  F*  5  plays  the 
role  of  an  "inertia  force"  applied  to  the  system  of  bodies.  The  reader  will  note  that 
the  D’Almbert  form  assumes  exactly  the  general  format  of  equation  El. 2.  Not  only 
does  (*)  serve  in  solid  mechanics  and  electric  theory  but  it  also  remedies  an  omission  in 
thermal  and  fluid  mechanics.  In  both  these  disciplines  absolute  quantities  (such  as 
temperature  and  pressure)  are  verbalized  but  never  symbolized  differently  from 
relative  quantities.  An  absolute  reference  state  (*)  permits  one  to  express  absolute 
pressure  at  state  (1)  as  P^*  while  gage  pressure  would  be  P,/a  wherein  subscript  (a) 
identifies  the  atmosphere  as  a  reference  state.  Similarly  absolute  temperature  at  state 
I  could  be  identified  as  T\/*  while  Tj/f  might  be  the  temperature  on  a  Fahrenheit 
scale.  Extended  tables  of  Through,  Across  Variables  and  Constitutive  laws  are  given 
in  Appendix  I. 
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Influence  Coefficients 

The  above  discussions  dwell  on  the  relationships  between  adjacent  terminals  in  a 
field.  Complementing  this  tool  of  analytics  is  the  mutual  influence  concept  for 
modeling  linear  systems.  In  this  concept  a  variable  at  a  particular  point  in  a  field  is 
said  to  be  influenced  by  a  combination  of  variables  at  remote  points  in  the  field,  that 
is: 


em/o  ^-*m,2fi2/o  "T  ^'in.r/’n/o 


Ef..r 


>a 


or 


em/o 


N 


^m,j  ej/o 

j“l 


El. 5b 


Definition:  Gmj  Is  an  influence  coefficient  defining  the  effect  on  a  quantity  at 

location  m  by  a  quantity  at  location  j. 


ej/0:  The  value  of  an  Across  variable  at  point  j  relative  to  a  common 
reference  point  o. 


Equation  El. 5a  is  still  one  more  form  of  summation  useful  to  modeling.  It  usually 
appears  in  the  more  ambiguous  form  em  =  V  Gm j  Cj  in  the  literature.  Coefficient  Gm j 
is  another  doubly  subscripted  property  but  now  with  a  (,)  subscript  separator.  Like 
the  (.)  and  (/)  subscripted  variables  the  order  of  the  subscripts  in  this  case  is 
consequential.  The  influence  coefficients  themselves  must  be  generated  from  solutions 
of  equations  modeled  by  the  constitutive  laws  described  earlier  or  from  some  extremum 
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principle.  Across  and  Through  variables  can  also  be  expressed  as  the  mutual  influence 
on  each  other  i.e.. 

em/o  —  or  ^.o  ~  Xi  ^s,jej/o 

j 

An  array  of  such  simultaneous  equations  can  be  expressed  succinctly  in  matrix 
form 

(em/o)  =  lGmjl{«lj.o)  iO  =  EI.6 

The  similarity  of  the  matrix  form  and  the  constitutive  equations  illustrated  in 
TABLE  II  is  a  telling  and  convenient  modeling  tool.  The  similarity  of  forms  is  usually 
not  recognized  because  the  datum  subscript  (o)  is  conventionally  dropped  in  the 
matrix  format  as  in  EI.7. 

(*»}  =  |G.jl  {q,l  bi.7 

ADDITIONAL  FIELD  MECHANICS  SYMBOLS  AND  FORMATS 

In  the  discussions  above  the  subscripts  took  their  identities  from  specific  points  in 
bodies.  Considerations  of  machine  kinematics,  metrology,  dynamics  and  even 
thermodynamic  systems  will  call  for  additional  indentifiers  such  as  bodies,  coordinate 
systems,  states  etc..  The  following  is  a  glossary  of  useful  expressions  and  their 
shorthand  forms  used  in  the  text  that  follows.  Once  this  nomenclature  is  mastered  it 
will  serve  to  reduce  ambiguity  in  verbalized  statements  and  mathematical  expressions. 
The  notation  suggested  earlier  and  in  the  glossary  below  will  also  expedite  the  creation 
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of  algorithmic  formats  which  help  in  developing  mathematical  expressions. 

NOTATIONAL  SYSTEM 

Numerals  0,  1,  2,  3,  ...  or  letters  a,  b,  c,  ...  will  identify  points  on  bodies  or 
thermodynamic  states. 

Roman  numerals,  □,  I,  II,  III....  and  letters  A,  B,  C,  etc.  will  represent  individual 
bodies  and  coordinate  systems  or  thermodynamic  media. 


=  R(8/7).  A  position  space  vector  of  point  8  relative  to  point  7.  For 

purposes  of  kinematic  analysis  the  prefix  R  is  superfluous,  since 
positioning  is  the  nature  of  the  subject,  hence  it  will  be  dropped.  On 
the  other  hand  when  thermodynamic  analysis  is  involved  a  prefix 
describing  the  property  will  be  retained  ie  h(8/7)  —  Value  of 
enthalpy  at  state  8  relative  to  its  value  at  state  7. 

(877)h  =  8/7  expressed  in  coordinate  system  II.  If  thermodynamic  analysis  is 

involved,  II  would  identify  the  medium  possessing  the  property  whose 
change  is  being  described. 

Xu, Yu,Zn  =  X,  Y,  Z  components  of  a  vector  expressed  in  coordinate  system  II. 


Xu 

V„ 

Zu 


Matrix  form  of  (8/7)h 


J8/7 
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(3/2) 


j  3/2 1  = 
|3/2|  = 


(372)  = 

X,„  Yj|,  Z„  — 


Velocity  vector  of  point  3  relative  to  point  2.  If  a  system  is  static 
the  same  form  will  be  used  for  virtual  displacements  of  3  relative  to 
2.  It  will  also  be  used  to  represent  static  errors  in  the  position  of  3 
relative  to  2.  The  context  of  the  problem  should  be  the  guide  to 
choice  of  interpretation, 

The  length  of  vector  (3/2). 

Rate  of  change  of  length  of  vector  joining  point  3  to  point  2.  This  is 
a  scalar  quantity.  If  system  is  static  the  same  symbol  will  be  use  to 
represent  a  virtual  change  in  the  length  of  (3/2).  The  same  notation 
will  also  be  used  to  represent  an  error  in  length  !°/2j. 

Unit  vector  parallel  to  (3/2). 

Unit  vectors  of  coordinate  system  II. 


n(III,II)  =  Point  n  in  space,  such  as  a  pin  connector,  shared  by  systems  III  and 

II. 


(372)  - 

*  ___ 

(37*)  = 


Acceleration  vector,  point  3  relative  to  point  2. 
Absolute  reference  point.  Also  ground. 
Absolute  acceleration  of  point  3. 


Angles  and  angular  motion  will  be  represented  by  Greek  letters  and  the  identity  of 
the  bodies  separated  by  the  angle. 
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|  311/*)  = 


w(ii/m)  = 


aM.  #11.  7u  = 


BpCiAn)  = 


*tX,/Z„)  = 


C(q,s)  = 


[C,,k 


Absolute  angular  velocity  of  body  II  or  angular  velocity  of  body  II 
relative  to  a  nonrotating  body.  This  notation  will  also  be  used  to 
represent  the  the  absolute  virtual  angular  displacement  of  a  static 
system.  Alternatively  this  symbol  will  be  used  to  represent 
absolute  angular  errors  in  the  position  of  body  II. 

Related  to  above  but  non-absolute  motion  is  implied.  Specifically 
it  is  the  angular  velocity  (or  virtual  displacement)  or  system  II 
relative  to  system  III.  It  follows  that:  cI(II/HI)  —  ..(II/*)  - 

-tm/*) 

X,  Y,  Z  components  of  c3(Il/*). 

Note  ft jj  =  o.(ny«).  Yj|. 

Angular  vector  position  of  axis  Xj  relative  to  axis  Zn.  Illustrated 
on  Fig.  1.3. 

Time  rate  of  changes  of  F(Xj/Zji).  Alternatively  virtual  change  or 
error  in  F{X[/Z|j). 

Influence  coefficient  of  a  property  applied  to  point(s)  on  the 
property  created  at  point  (q). 

An  array  of  many  influence  coefficients,  such  as  C(q,s),  in  a  matrix 
form. 
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COORDINATE  TRANSFORMATION 

Perhaps  the  most  frequent  application  of  El. 6  and  E.I.7  in  machine  tool  and 
precision  engineering  is  for  coordinate  transfers  from  one  subassembly  to  another. 
This  equation  appears  again  as  a  basic  of  the  finite  element  method  of  analysis.  To 
illustrate  the  application  of  modeling  forms  described  in  this  chapter  consider  the 
derivation  of  the  familiar  coordinate  transfer  matrix.  As  an  example  assume  that  one 
is  given  a  point  P  in  the  spatial  coordinates  of  system  I,  Fig.  1.3,  and  it  is  desired  to 
restate  the  position  of  P  in  terms  of  system  II.  Furthermore  the  relationship  of  the 
two  space  systems  is  given  in  the  form  of  the  rotations  of  the  axes  of  II  relative  to  the 
axes  of  1  (such  as  #(Xu/Y])  of  Fig.  1.3.  Also  the  position  of  the  origin  (Ojj)  of  system  II 
relative  to  the  origin  (0|)  of  system  I  is  known.  That  is  the  vector  0n70|  is  known. 
However  the  coordinates  of  this  vector  are  given  in  terms  of  system  II. 

Referring  to  Fig.  1.3,  the  rules  for  subscripture  of  across  variables,  El. 3b  and  the 
NOTATIONAL  SYSTEM  the  statements  above  are  converted  to  the  following 
symbolic  forms. 


Fig.  1.3  Coordinate  Systems  I  and  II  Describing  Position  of  P. 
Relative  Position  Expressions  Illustrated. 
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GIVEN:  Yj  ...  The  vector  P/Oj  in  system  I  coordinates 

h  P/o. 


...The  vector  Oji/Oj  describing  the  position  of  origin 
of  system  II  relative  to  the  orgin  of  system  I. 

It  is  expressed  in  system  II  coordinates. 


H  Jo.,/0, 


0(X„/X() 

...  The  angles  that  the  coordinates  of  system  II 

Q(y  /y  \  °  J 

1  11/  l)  make  relative  to  the  coordinates  of  system  I. 

fl(Zu/Z,) 


DESIRED:  P/Ojj  that  is  Yu 

Z"  Jp/0.. 


SOLUTION:  P/0„  =  P/Oj  +  0^0,1  =  P/fy  -  0n/0, 


that  is: 


=  Y, 


-  Y, 


J”  Jp/o„ 


JI1  lp/0, 


11  J0u/0| 


Matrix  El. 9  is  a  more  complete  expression  of  El. 8  since  X|It  Yjj,  Zn  emphasize  the 
requirement  that  the  equation  be  written  in  system  II  coordinates  exclusively. 
Meanwhile  the  subscripts  at  the  base  of  the  column  matrix  inform  one  that  the  terms 
within  the  matrix  deal  with  the  relative  distances  between  the  two  points  cited.  For 
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the  column  on  l.h.s.  of  EI.9,  the  footnote  indicates  that  elements  of  the  matrix  are 
actually  the  coordinates  of  P  relative  to  the  origin  0tl. 

Part  of  the  given  data  of  this  example  appears  as  the  last  vector  of  El. 9.  On  the 
other  hand  the  first  vector  on  the  r.h.s.  of  EI.9  is  not  among  the  given  data.  What  is 
actually  given  is  P/0,  in  system  1  coordinates.  The  required  vector  can  be  synthesized 
by  a  coordinate  transformation  in  the  form  of  El. 7  which  in  this  case  becomes 

Xu  X! 

Yn  =  [cn,il  Yj  El. 10 

Zll  p/o,  ZI  P/o, 

A  particular  element  of  the  matrix  [Cni],  Cj  j,  is  the  cosine  of  the  angle  that  axis  i 
of  system  II  makes  relative  to  axis  j  of  system  I.  That  is  Cs j  —  cos  0(i/j).  Pig-  1.3 
illustrates  the  physical  meaning  of  some  terms  0(i/j).  For  compactness  the  following 
notation  will  be  used: 

Let  C  i/j  =  cos  0  (  axis  i  relative  to  axis  j) 


Hence  C  (Xj/Zji)  =  Cosine  of  the  angle  that  axis  X|  makes  relative  to  axis 

Zip 

Similarly  S  (Xi/Zn)  =  Sine  of  the  angle  that  axis  X,  makes  relative  to  axis 

zll 

If  the  direction  of  the  relative  angular  positions  is  given  graphically  such  as 
^(Xu,Yj)  of  Fig.  1.3  the  prior  subscript  of  0( i/j)  is  the  axis  touched  by  the  arrow  or  the 
arc  vector.  The  posterior  subscript  is  the  axis  touched  by  the  tail  of  the  arc  vector. 


m 
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where: 


a 


The  vector  of  the  offsets  of  the 
coordinate  origins  relative  to  each 
other. 


b 

c 


-x„ 

-Y„ 

-Zji 


Jo„/o, 


Infinitesimal  rotations  of  vectors  are  common  in  kinematics  and  error  analysis  of 
machine  systems.  With  infinitesimal  rotations  the  transfer  matrix  El.  11  is  simplified  to 
EI.13. 


i  -,(ii/i)  -Mn/i) 

fl(II/I)  -a(II/l)  1 


EI.13 
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APPENDIX:  AI 

"Modelling  Aida  and  Analogues  With  Doubly  Subscripted  Terms" 


Across  type  of  variables  are  identified  by  a  slash  separating  the  two  subscripts 
identifying  terminals  at  the  end  of  an  element.  Through  type  variables  are  identified 
by  a  dot  separating  these  two  subscripts.  A  dash  separates  the  two  subscripts  of  the 
property  of  the  element  between  the  two  terminals  TABLES  A.I.a,b,c,d  organizes  some 
of  the  basic  Across,  Through  variables  and  element  properties.  These  are  related  by 
the  appropriate  constitutive  laws  for  electrical,  mechanical,  fluid  flow  and  heat 
transfer  mechanics.  The  relationships  of  TABLES  A.I  are  useful  for  modeling  a  class 
of  problems  generally  expressed  by  the  partial  differential  equation  i 


d  .  dA.  d  ,,  3  A  » 

aT(a  aT> +  ^(b^> + 


3  ,  3  A  .  -  ,  ,  3A  tf  . 
-(c— )  +  d  A  =  e—  +  f(x,y,i) 


where  A  =  a  generic  symbol  for  across  variables.  j 


When  a  continuum  field  is  approximated  by  a  network  of  discrete  elements  such  as 
those  from  TABLE  I,  Finite  Difference  equations  are  directly  generated  without  first 
having  to  go  through  the  process  of  determining  the  exact  partial  differential  equation 
of  the  process  and  then  decomposing  it  into  the  finite  difference  form.  The  property 
terms  such  as  Ka_b  must  be  determined  for  the  specific  case.  For  instance  for  heat 

n 

conduction  through  a  thin  plate  made  of  a  material  with  a  thermal  conductivity  < 

< 

k(Btu/deg  sec  in)  and  a  conducting  area  Aa_i,  between  nodes  a  and  b,  conductance  | 


would  be  expressed  as  Ka_t,  = 


kB.t 


for  a  cartesian  net. 


-  24  - 


AI.2 


where:  B.t  =  Aa_b,  mean  cross-sectional  area  of  the  conducting  branch  between  a 
and  b 

L  =  L(a_b)  the  length  of  the  branch  connecting  points  a  and  b. 

The  same  strategy  can  also  lead  to  a  LINEAR  FINITE  ELEMENT  solution  of  field 
problems.  This  method  and  the  determination  of  Ka_b  for  an  approximating  mesh  of 
irregular  triangles  is  described  in  Appendix  A.IV.3. 

TABLES  AI.l  offer  iconic  symbols  for  the  elements  of  the  different  disciplines. 
Iconic  representation  in  network  drawings  enhance  communication  of  concepts  and  has 
long  successful  history  in  electrical  engineering.  There  seems  to  be  no  such  standard 
of  representation  in  a  field  such  as  heat  transfer.  TABLE  A.I.d  offers  a  possible  set  of 
iconics  however  it  is  more  convenient  to  borrow  the  symbols  from  analogous  electrical 
elements  for  depicting  heat  transfer  nets.  Fig.  A.I.l  is  an  analog  wheel  which  helps 
establish  the  relationships  between  different  physical  systems.  Both  mobility  and  clas¬ 
sical  analogs  as  well  as  dualities  are  depicted. 

The  mobility  analogue  (nodes  go  to  nodes  and  loops  go  to  loops)  is  generated  by 
taking  the  symbols  along  a  bold  radial  line  of  Fig.  A.l  and  transforming  them  into  the 
symbols  on  another  bold  radial  line.  As  an  example  a  mechanical  element  connecting 
ends  a  and  b  has  an  electrical  analog  between  electrical  terminals  a  and  b.  The 
mechanical  terms  on  the  radial  line  at  11  o’clock  are  replaced  with  the  electrical  terms 

at  3  o’clock.  Thus  spring  rate  K  is  transformed  to  the  reciprocal  of  inductance,  — . 

1 . 

The  classical  analogue  (nodes  transform  into  loops  and  loops  into  nodes)  is  done  by 
taking  the  symbols  along  a  bold  radial  line  of  the  analogue  wheel  and  transforming 


-  25  - 


AI.3 


them  into  the  symbols  along  a  light  radial  line  for  the  other  physical  system.  As  an 
example  a  mechanical  element  connecting  points  a  and  b  has  an  analog  in  the  form  of 
an  electrical  element  bordered  by  loops  a  and  b.  The  mechanical  elements  on  the  bold 
radial  line  at  11  o’clock  now  have  their  analogues  elements  on  the  light  radial  line  at  5 

o’clock.  Thus  the  mechanical  spring  with  rate  K  is  transformed  to  — ,  the  reciprocal 

C 

of  electrical  capacitance.  The  easy  manipulation  of  analog  and  dualities  is  not  only  a 
very  valuable  tool  for  the  numerical  modeling  of  physical  phenomena  but  it  is  also  a 
useful  aid  in  consulations  with  specialists  outside  one’s  particular  expertise. 


I 


'.l,*  I 

'»,*1 


k  c  ill 

— —  “ED"-  “O  * 


MECHANICAL 


CLASSICAL  ANALOGUE  C 


DUAL 


^K\W 

Vc  N  CLASSICAL  ANALOGUE 


c  hi  K  e  I  l 


DUAL 


DUAL 


FLUID 


ELECTRICAL 


-nsw^-L 


CLASSICAL  ANALOGUE 

a : node 
X  :  loop 


Fig.  AI.l  Analogue  and  Dualities  Wheel 
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II.  GUIDE- WAYS 

The  essence  of  the  quality  of  a  machine  tool  or  a  mechanical  actuator  is  the 
precision  of  its  guide  ways.  The  machine  tool  of  Fig.  II.  1  illustrates  the  implicit  role  of 
several  interconnected  guide  ways  and  sliding  tables.  The  slides  are  driven  in  a 
fashion  as  to  move  the  workpiece  and  cutting  edge  in  three  dimensional  space.  The 
expectation  is  that  point  7  of  the  tool  comes  in  in  contact  with  a  predetermined  point 
8  of  the  work  piece.  If  that  is  actually  the  case  the  relative  position  vector  is  8/7  — 
0.  This  defines  ideal  conditions.  In  reality  each  of  the  interconnected  bodies  suffers 
errors  which  combined  produce  a  mismatch  and  therefore  an  error  in  the  work  equal  to 
the  error  vector  *(8/7).  The  error  (8/7)  is  itself  a  chain  rule  summation  of  error 
contributions  by  all  the  connecting  members  of  Fig.  II.  1.  It  is  therefore  necessary  to 
examine  each  of  the  component  errors  in  turn. 


To  illustrate  the  nature  of  localized  kinematic  errors  consider  the  limited  case  of 
bodies  I  and  II  sliding  on  each  other  without  elastic  or  thermal  distortions.  Figure  II. 2 
is  a  rudimentary  representation  of  the  pair  with  Fig.  II.2a  representing  an  ideal 
geometry  and  Figs.  II  2b,  c,  d  depicting  various  malformed  guide  ways.  Consider 
measurements  taken  with  table  II  tixed  and  table  I  positioned  along  the  X  axis  of  II  by 

1  Refer  to  NOTATIONAL  SYSTEM  of  Section  I  for  interpretation  of  symbols  and  abbreviations. 


REFERENCE  POINT 


Name  of  bod 


Svmbol 


Main  frame 
X  Slide 
Y  Slide 
Z  Slide 
Work  Piece 
Tool  or  Probe 
Spindle 


Point 

Host 

0 

□ 

1 

□ 

2 

III 

3 

□ 

4 

II 

5 

I 

6 

S 

7 

8 

T 

W 

Fig.  II.  1  Machine  With  Composite  Guidance  By  Components. 
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Resultant 
Rotational 
Error 

Resultant 

Error 

Transmitted 


Fig.  11.3  Errors  Vectors  of  Table  II  Relative  to  Table  I  due  to  a  Malformed  Guide. 


the  actuator  F.  Control  imperfections  of  this  feed  create  the  first  observable  error, 
that  being  along  X  as  in  Fig.  II. 2b.  Points  4  and  5  are  reference  points  which  would 
serve  as  gaging  points  for  measuring  this  error.  The  observed  error  is  therefore 
identified  as  XII(5/4).  It  is  the  error  component  along  the  X  axis  of  II  due  to  a  pure 
translatory  error  between  I  and  II.  The  actual  points  5  and  4  would  normally  be 
chosen  as  convenient  stations  at  which  to  mount  the  two  parts  of  an  experimental 
measuring  gage.  While  analytical  convenience  might  suggest  other  reference  points 
the  physical  limitation  of  locating  instruments  must  prevail. 

Figure  II. 2c  shows  the  error  imposed  by  a  skewed  guide.  The  error  is  a  rotation, 
around  the  Z1I  axis,  of  II  about  I,  that  is  ^II(I/II). 

The  curved  guide  of  Fig.  II. 2d  causes  a  sideways  translation  of  table  II  hence  an 
error  Yll  (1/11).  In  real  world  examples  the  guides  of  Figs  II. 2b  and  .2c  become  one  in 
the  form  of  a  skewed  curved  guide  as  in  Fig.  II. 3. 

For  convenience  only  a  view  along  the  Z  axis  was  illustrated  in  Fig.  II. 2.  Similar 
studies  would  need  to  be  made  by  sights  along  the  X  axis  of  body  II  as  well  as  along 
its  Y  axis.  Combining  the  three  studies  would  enable  one  to  extract  the  six  com¬ 
ponents  which  would  describe  the  total  translatory  and  rotary  errors  between  the 
members  I  and  II  due  to  the  geometric  errors  in  the  mutual  guide  way.  These  errors 
are  XI  1(5/4),  YU(5/4),  ZII(5/4)  and  r'rII(l/II),  ?fll(l/ll),  ^11(1/11).  If  similar  studies  are 
made  of  all  the  pairs  of  bodies  of  Fig.  II. 1  the  results  could  be  combined  in  a  chain 
form  as  suggested  by  eqtn  El. 3b  to  arrive  at  the  total  error  (8/7)  described  by  cqtn 
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(8/7)  -  (8/2)  +  (2/1)  +  (1/3)  +  (3/4)  +  (4/5)  +  (5/6)  +  (6/7)  EII.l 

Actually  several  other  types  of  error  sources,  besides  the  as-manufactured-errors 

just  described,  must  be  considered  for  eqtn  Eli.  1 .  A  listing  of  this  type  of  error  and 

some  of  the  more  important  agents  of  error  is: 

1.  Geometric  or  rigid  member  kinematics  errors.  These  arise  from  imperfect  fabri¬ 
cation  and  positioning  of  assembled  members.  These  imperfections  manifest 
themselves  as  out  of  straightness,  skewed  and  twisted  guides,  displaced  con¬ 
straints,  non  orthogonality  of  axes  etc..  These  errors  are  essentially  invariant 
with  time. 

2.  Wear  and  backlash  are  temporal  factors  whereby  errors  resulting  from  a  mode  of 
operations  are  superimposed  on  the  geometric  errors. 

3.  Thermal  drift  represents  the  most  conspicuous  form  of  temporal  errors.  The 
kinematic  constants  now  vary  with  operating  time.  Errors  arise  from  the  chang¬ 
ing  thermal  expansion  and  contractions  of  the  components  involved. 

4.  Elastic  and  or  plastic  deformations  of  components  due  to  static  and  dynamic 
loads  caused  by  the  weights  of  the  system  members  and  forming  loads  applied  to 
work  pieces.  These  again  change  the  kinematic  constants  with  time  and  operat¬ 
ing  sequence. 

5.  In  addition  to  distortion  of  "solid"  members  implied  in  4  above,  there  are  fluid 
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Shapes  and  Precision  of  Linear  Guideways 

Commercial  linear  guideways  take  on  various  profiles  as  in  Fig.  II.4.  The  funda¬ 
mental  objective  of  such  slideways  is  to  provide  only  one  degree  of  freedom  that  being 
unrestricted  and  pure  translation  along  the  X  axis  (refer  to  Fig.  II.2a,b  and  Fig.  11.4). 
Existence  of  any  movements  in  the  other  five  coordinates  are  recognized  as  errors.  For 
illustrative  clarity  the  guideways  shown  in  Fig.  II.  1  are  indicated  as  squared-off  sur¬ 
faces  (surfaces  1,3,4).  Other  possible  commercial  forms  are  depicted  in  Fig.  II.4.  The 
squared-off  profile  (a)  is  a  very  popular  choice  for  horizontal  installations  being  a  good 
compromise  between  accuracy  and  fabrication  costs.  The  double  "V",  Fig.  II. 4e,  offers 
the  highest  accuracy  possibilities  but  with  the  highest  fabrication  costs  of  the  series 
shown. 

As  depicted  in  Fig.  II. 4  only  shapes  (b)  and  (c)  seem  to  have  the  necessary  five  con¬ 
straints  needed  for  the  ideal  linear  guide.  Tables  (a),(d),(e)  have  the  theoretical  free¬ 
dom  to  lift  off  in  the  Z  direction  as  well  as  to  rotate  in  the  (3  mode  hence  only  three 
restraints  would  remain.  However  the  practical  application  of  these  profiles  is  such 
that  the  external  force  W  preloads  the  surfaces  sufficiently  to  force  contact  at  all 
times  so  that  live  constraints  are  actually  preserved. 

Earlier  the  agents  of  error  in  guideways  were  listed.  Let  us  use  some  of  the  listings 


to  compare  the  designs  of  Fig.  II. 4  by  way  of  the  score  sheet  TABLE  II.  1. 


TABLE  II.  1 


Relative  Potentials  for  the  Minimization  of  Error  Sources 

(5  =  Best  1  =  Poorest) 

Guideway  Design 

Error  Source 

a 

b 

c 

d 

e 

Geometric 

3 

2 

1 

4 

5 

Wear 

5 

2 

I 

4 

4 

Play 

3 

2 

1 

5 

5 

Flexibility 

5 

5 

1 

5 

5 

Play  (Back  Lash)  Guidance  Error 

The  sliding  interface  surfaces  of  Fig.  II.4  are  identified  by  the  bolder  lines.  In 
design  (d)  the  load  W  forces  intimate  contact  at  these  surfaces.  By  contrast  W  will 
not  close  up  the  inherent  clearance  on  the  sides  profile  of  (c).  There  is  therefore  an 
opportunity  for  undesirable  movement  in  the  Y  and  7  modes.  The  apparent  one 
degree  of  freedom  of  design  (c)  is  actually  three  degrees  of  freedom  even  with  load  W 
applied.  Hence  the  poor  score,  of  1,  for  play  control  in  design  c.  Consider  the 
squared-off,  design  (a),  which  must  also  have  an  assembly  clearance  between  its  side 
walls.  That  and  the  tolerance  on  parallelism  of  the  side  walls  will  result  in  some  Y 
and  ")  play.  To  raise  the  performance  of  (a)  to  a  score  of  3  a  gib  (tailored  to  a  partic¬ 
ular  pair  of  members  I,  II  as  to  minimize  the  clearance  and  out  o"  parallelism)  is 
inserted  on  one  of  the  side  walls.  Dove  tail  guide  (b)  has  the  same  problem  just 
described.  Forming  a  quality  dove  tail  and  a  possible  gib  is  obviously  a  more  expen¬ 
sive  design  than  (a)  but  a  dove  tail  guide  is  a  necessary  alternative  to  (a)  if  the  load 


W  is  insufficient  to  force  full  contact  on  the  horizontal  surfaces. 
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Flexibility 

Rail  designs  a,b,d,e  impose  no  particular  limitations  on  the  stiffness  of  the  support¬ 
ing  bed  II.  Guide  design  (c)  is  however  very  much  more  flexible  in  bending  because  the 
cylindrical  rails  are  supported  only  at  their  ends.  The  free  span,  approximately  the 
length  of  carriage  travel,  is  given  to  excessive  bending.  Supporting  these  rails  along 
their  entire  lengths,  as  in  the  supplementary  detail  to  Fig.  II. 4c,  helps  to  reduce  flexure 
in  the  vertical  and  horizontal  planes  but  the  latter  mode  is  still  severely  compromised 
by  the  thinness  of  the  fin. 

Achieving  Geometric  Accuracy 

W.  R.  Moore^  gives  an  excellent  discourse  on  the  mechanics  of  fabrication  of  accu¬ 
rate  guide  ways.  Very  persuasive  arguments  are  offered  for  the  superiority  of  guide 
profiles  (d)  and  (c).  Quality  of  linear  guides  stems  from  the  accuracy  of  master  sur¬ 
faces  and  straight  edges  which  arc  used  as  gauges  and  laps.  Accuracy  starts  with  the 
generation  of  master  flat  surfaces  by  the  three  plate  method.  In  this  three  nominally 
flat  plates  are  abrasively  rubbed  against  each  other  in  a  systematic  cross  matching 
and  successive  90°  turning  on  each  other.  The  plates  in  a  natural  fashion  work  each 
other  into  highly  accurate  plane  surfaces.  Flatness  accuracy  is  enhanced  still  further 
by  reducing  residual  hills  and  valleys  by  hand  scraping  and  local  lapping.  The  objec¬ 
tive  is  to  create  a  surface  flat  to  approximately  0.6  //rn  over  a  span  of  1000  mm.  From 

(1)  Cited  references  are  listed  at  the  end  of  each  section. 
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these  master  surfaces  four  sided  master  straight  edges  are  derived  Fig.  II.5.  The  four 
sided  (box)  geometry  possesses  several  critical  advantages  for  attaining  and  maintain¬ 
ing  straightness  and  parallelism.  The  axi-symmetric  form  insures  maximum  flexural 
stiffness  in  all  positions  as  well  is  minimum  thermal  distortion.  The  process  of  genera¬ 
tion  also  results  in  very  accurate  90°  corners.  By  setting  the  box  straight  edge  on  the 
master  surface  plate  the  eight  gage  surfaces  can  be  effectively  checked  for  straightness, 
squareness  and  parallelism  with  respect  to  the  master  surface  and  therefore  with 
respect  to  each  other.  Finally  the  90°  corner  and  straightness  are  self-proving  in  the 
generation  of  the  masters  and  the  machine  ways.  This  cannot  be  said  of  straight 
edges  with  non  square  rectangular  cross-sections  or  of  V  shaped  straight  edges. 
Because  the  accurate  90°  corners  derive  naturally  from  matching  the  box  straight  edge 
to  the  master  surface  plate,  these  corners  are  then  the  ideal  tool  for  marking  and  lap¬ 
ping  female  V  guides  Fig.  II.5.  By  tipping  the  box  master  straight  edge  end  for  end  as 
well  as  rotating  sequentially  through  90°  while  checking  for  high  spots  or  lapping,  a 
female  V  guide  can  be  formed  which  is  truer  than  any  of  the  individual  eight  working 
surfaces  of  the  box  straight  edge.  So  it  is,  because  it  is  naturally  derived  from  gaging 
and  forming  surfaces  inherently  superior  in  accuracy,  the  90°  V  track  is  a  superior 
guideway.  The  double  V  of  (e)  has  a  slight  advantage  over  (d).  The  V  and  flat  tracks 
wear  at  somewhat  different  rates,  hence  twist  error  is  somewhat  less  likely  in  symmetr¬ 
ical  arrangements  such  as  e. 


I 


Cast  iron  is  an  excellent  material  for  masters  and  machine  guideways  because  it 
responds  well  to  scraping.  Furthermore  cast  iron  machine  beds  possess  good  internal 


vibration  damping.  On  the  other  had  cast  iron  sliding  on  cast  iron  wears  many  times  i 

faster  then  cast  iron  on  steel  (W,  Moore^  observes  a  10:1  wear  ratio).  Hardened  and  } 
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ground  inserts,  Fig.  II.6,  add  wear  resistance  to  the  latent  accuracy  of  the  V  guide 


concept.  In  this  arrangement  the  female  V’s  of  the  C.  I.  base  are  scraped  to  a  male 


double  V  master  to  about  2.5  jim  straightness  and  parallelism.  Although  these  V’s  will 


not  be  used  as  final  slide  ways  the  accuracy  will  minimize  possible  distortions  created 


by  bolting  the  steel  inserts  to  the  bed.  The  upper  C.I.  sliding  table  is  prepared  in  the 


same  fashion  as  described  for  the  bed  V'ways.  After  the  steel  rails  are  ground  true  in 


assembly,  the  upper  table  is  lapped  into  the  steel  guides  achieving  accuracy  by  match- 


Upper  C.I.  Slide  Table 


Ground  Steel  Inserts  (Rails) 


C.I.  Base 


Fig.  II. 6  Guideway  With  Precision  Rails 


Geometric  Calibration  of  Precision  Mechanical  Systems 


Superior  accuracy  of  each  component  is  fundamental  to  accurate  positioning  by 


mechanical  systems.  Yet  the  accumulation  of  errors,  as  suggested  by  Eli.  1  flaws  the 


over-all  accuracy.  On  the  other  hand  over-all  repeatability  from  an  assembly  of  accu¬ 


rate  components  could  be  very  high.  The  principle  objective  of  precise  and  accurate 


performance  can  be  served  by  calibrating  a  full  machine  system.  Several  calibration 
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techniques  are  described  in  publications  (2, 3, 4, 5, 6). 

Using  the  minimal  number  of  reference  points,  the  N.B.S.  method^  involves  two 
interferometer  metrology  stations  A  and  D,  Fig.  II. 7.  A  third  metrology  station  G 
(representing  work  piece  W  of  Fig.  II.  1)  contains  gage  points  such  as  g  one  of  which 
could  be  the  surrogate  for  point  8  of  the  work  piece  W.  Metrology  posts  A  and  B  vec¬ 
tor  onto  each  other,  the  gage  point  g  and  onto  the  probe  p.  The  latter  is  the  surro¬ 
gate  for  the  tool  tip  7  of  Fig.  II.  1.  There  are  three  coordinate  systems  A,  B,  G  and 
their  origins  are  at  a,  b,  g  respectively.  For  analytic  simplification  the  three  coordi¬ 
nate  systems  are  set  parallel  to  each  other  except  for  infinitesimal  rotations. 

The  metrology  bases  A  and  B  vector  on  each  other,  on  (g)  and  on  (p)  to  generate 
position  vectors  (p/a)A,  (b/a)A  and  (g/v)B  throughout  the  range  of  motions  of  tables  I, 
II,  III  and  point  8  of  Fig.  II.  1.  The  range  of  these  observations  is  contained  in  a  rec¬ 
tangular  parallel  piped  envelope  whose  volume  is  divided  into  a  cubic  lattice  with  a 
mesh  size  of  about  20%  of  the  smallest  dimension  of  the  envelope.  Position  vectors  are 
recorded  for  all  the  nodes  of  this  mesh  field. 

The  N.B.S.  calibration  method  aims  to  express  the  vector  (p/g)  in  the  coordinates 
of  the  gage  piece  G.  Employing  the  summation  principle  of  El. 3b  the  required  vector 
(p/g)c  would  be  synthesized  from  the  measured  vectors  as  in  Eli. 2 

(p/g)o  =  (pA)g  +  (a/b)o  +  (b/e)c 


or  (p/g)c  =  (b7a)c;  “  (g/b)c 


Ell. 2 


Fig.  II. 7  Machine  of  Fig.  II. 1  With  Calibration  Stations. 
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Because  the  coordinate  axes  of  the  three  metrology  bases  are  off-parallel  from  each 
other  by  minute  rotations  the  measured  vectors  must  be  transformed  to  the  base  G  as 
follows 


IpA}c  =  Icg,aI  {pAJa  =  Icg,b]  [Cb,a]  {pA}a 
{b/a}G  =  [Cg,a!  t  b  A  }a  =  [Cg.bI  [cb,a1  {^A)a 

{g/Mc  =  (Cq.bI  {g/b^B 

Substituting  these  transformations  into  EII.2 

(pA)g  -  lCG.n]  |Cb.aI  I(pA)a  -  (BA)a1  -|cg.b]  (g/b)B  Eli. 3 

where: 


(p/g)o  “ 

^p/g 

YP/g 

*  (pA)a  = 

Yp/a 

Yp/a 

f 

Zp/g 

G 

ZP/a 

A. 

^b/a 

Yg/b 

(b/a)A  = 

Yb/a 

»  (g/b)B  = 

Yg/b 

Zb/a 

A 

Zg/b 

B 

I 

l 

I 
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[CB,a]  - 


—  'JG/B  1  <yC 

*  e 

Pg/B  -(yG/B  1 


1  Tb/a  >n/j 

*)B/A  1  <VB/A 

/^B/A  — °B/A  1 


Equation  Eli. 3  implemented  with  the  data  from  metrology  posts  A  and  B  yields  the 
coordinates  of  the  probe  point  p  (or  tool  point  7)  relative  to  a  reference  point  in  the 
gage  G  (or  a  reference  point  in  work  piece  W).  This  information  is  in  the  coordinates 
of  the  work  piece.  These  data  are  stored  with  corresponding  machine  measurements 
of  the  positions  of  tables  I,  II,  and  III.  This  record  is  then  the  calibration  of  the 
machine  implicit  in  which  are  the  actual  built  in  errors  described  in  Figs.  II. 2,  3  and 
equation  Eli.  1 .  This  calibration  is  then  used  to  automatically  compensate  the  machine 
controls  during  the  production  mode. 


‘ '/  V ■  VJ*  '>*  ♦  *  ' 
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III.  FRICTION  OF  THIN  FILM  SLIDEWAYS 


Precision  guideways  lapped  with  box-straight  edges  described  earlier  Fig.  II. 5  can 
be  made  straight  to  within  a  .001mm  camber  in  a  0.4  m  length^.  This  then  sets  a 
standard  by  which  operating  distortions  to  guideways  may  be  judged.  Thermal 
distortions,  structural  flexibility  back  lash  and  wear  are  the  most  serious  detractors 
from  precise  straight  line  guidance.  As  an  example  of  the  magnitude  of  distortion 
attributable  to  thermal  gradients  consider  the  case  of  a  cast  iron  guideway  represented 
by  a  prismatic  bar  100  mm  deep.  Assume  that  along  the  length  of  this  bar  there  is  a 
portion  0.4  m  long  in  which  there  is  a  temperature  drop  of  1°C  from  the  top  to  bottom 
surfaces.  Such  a  thermal  gradient  will  have  produced  a  1.8  //m  camber  in  the  .4  m 
length  and  a  slope  of  18  /v  rad  on  each  end  (Appendix  A. 111).  That  is  the  1°C 
temperature  difference  introduced  errors  1.8  x  the*  geometric  errors  inherent  in  the 
fabricated  guideway. 

A  potent  direct  and  or  indirect  source  of  such  thermal  distortions  is  the  frictional 
drag  of  sliding  surfaces.  Not  only  does  friction  cause  local  temperature  gradients  in 
the  sliding  surfaces  but  it  is  also  reflected  in  the  power  demands  on  the  driver  which  in 
turn  becomes  a  source  of  heat  contributing  to  the  thermal  distortion  of  the  structure. 
Vibration  related  to  friction  dynamics  is  another  error  generating  phenomenon. 
Patently  the  mechanics  of  friction  must  be  understood  for  rational  design  and 
diagnostics  of  precision  mechanisms. 

The  mechanics  of  drag  between  two  bodies  depends  on  the  details  of  the  interfac.ial 
zone  between  them.  Excluding  macro  electro-mechanical  effects,  five  interface 
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structures  are  normally  referred  to  in  the  study  of  sliding  surfaces.  Figs.  III. la  -  e 
show  these  basic  structures  and  the  coefficients  of  friction  roughly  associated  with 
them. 

Friction  Model  (a);  Fig.  III.  1 

The  identification  DRY  for  this  model  implies  that  there  is  no  deliberate  effort  to 
provide  a  fluid  lubricant  to  the  interface.  This  is  also  the  most  rudimentary  model  of 
the  five  cases  given.  Its  simplicity  gives  promise  of  lowest  first  cost  and  maximum  sup¬ 
port  rigidity.  On  the  other  hand  it  suffers  by  comparison  because  it  is  associated  with 
the  highest  expected  friction  and  wear.  While  the  purest  form  of  model  (a)  is  not  a 
promising  basis  for  precision  guides,  the  friction  mechanics  involved  is  worth  examin¬ 
ing  because  it  exists  at  times,  to  a  degree,  in  all  the  other  models  of  Fig.  III.  1  as  well. 

Many  theories  explaining  friction  between  pairs  of  metal  surfaces  prevade  past  and 
current  literature.  The  earliest  concept,  credited  to  Amontons  (1699),  pictures  drag 
resulting  from  interlocking  asperities  of  the  mating  surfaces.  Later  when  instruments 
were  developed  capable  of  measuring  the  microscopic  topography  of  smooth  surfaces 
the  results  gave  apparent  support  to  the  interlocking  concept.  Even  the  finest  finishes 
magnified  show  roughness  as  in  Fig.  111. 2b.  Such  apparent  serrations  certainly  suggest 
interlock  resistance  to  relative  motion  of  the  mating  elements.  On  the  other  hand  Fig. 
III. 2b  is  somewhat  misleading  because  this  form  of  record  is  usually  based  on  unequal 
scales.  The  horizontal  scale  is  often  shown  as  1%  of  the  vertical  scale.  Pig.  III. 2c 
shows  the  same  profile  with  equal  scales.  The  true  surface  contour  of  a  well  finished 
guideway  now  appears  in  the  microscopic  view  as  a  gently  rolling  meadow.  Rather 
than  breaking  through  a  barrier  of  slender  bristles  it  would  now  seem  more  plausible 
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Classes  of  Lubrication  by  Interface  Structures. 
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a)  Visually 
Smooth 


b)  Instrument  c)  Equalized 
Profile  Record  Scaling  of  b 


Fig.  III. 2  Effect  of  Choice  of  Scales  on  Comprehension 
of  Surface  Roughness. 


a)  Initial  Contact 


Contact  Area 
Mno  of 
Atparitlas 


b)  Bowden -Tabor 
Model 


c)  Edwards  -  Hailing 
Model 


Fig.  III. 3  Asperity  Contact  Geometry. 
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that  motion  between  mating  surfaces  could  be  accommodated  by  a  hill  climbing  and 


dropping  action.  Such  motion  is  however  energy  conservative,  whereas  friction  decid¬ 


edly  a  nonconservative  action.  An  alternate  concept  of  friction,  motivated  by  the  pic¬ 


ture  of  stubby  asperities,  is  the  adhesion  model  advanced  by  Bowden  and  Tabor^. 


The  B  &  T  model  is  the  basis  on  which  are  built  current  analytical  refinements  to  the 


theoretical  mechanics  of  friction  between  sliding  solids. 


Adhesion  Mechanics  of  Friction 


When  two  flat  visually  smooth  surfaces  are  pressed  together  the  nominal  pressure 


pa  equals  N/Aa  where  N  is  the  compressing  load  and  Aa  is  the  apparent  contact  area 


of  the  bodies  involved.  In  reality  the  load  is  supported  on  a  multitude  of  peaks  as 


depicted  in  Fig.  III.3b.  If  A'  is  the  area  of  contact  of  a  typical  asperity  pair  and  n  is 


the  number  of  such  contacting  pairs  then  the  real  total  area  of  contact  is  Ar  —  n  A'. 


It  follows  then  that  N  =  p  Ar  where  p  is  the  mean  pressure,  acting  on  the  n  contact 


areas,  sufficient  to  sustain  load  N.  As  the  load  N  is  increased  the  existing  contact 


points  soon  reach  plastic  flow  conditions.  Having  reached  their  load  carrying  limit,  the 


junction  points  offer  no  further  resistance  to  deformation.  Further  overall  deformation 


caused  by  an  increase  in  N  brings  new  asperities  into  load  sharing  contact.  There  is 


an  increase  to  the  total  real  area  in  proportion  to  the  increase  in  load.  This  hypothesis 


and  analogous  electrical  conductivity  tests^  on  pairs  of  flat  plates,  indicate  that  the 


total  real  area  Ar  increases  in  proportion  to  the  increase  in  applied  force  N  and  is 


independent  of  the  apparent  area  Aa.  Furthermore  the  actual  total  contact  area  at  a 


typical  full  load  never  reaches  more  than  a  minute  fraction  of  the  apparent  area. 


That  is  maximum  A.  ^  - A.. 
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Bowden-Tabor  First  Analytical  Approximations 
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In  this  initial  analytical  approximation  to  the  coefficient  of  friction  it  is  assumed 
that  the  stress  distribution  within  the  asperity  is  one  dimensional,  essentially  a  simple 
bar  in  compression.  Hence 


N  =  pcAr  pc  =  Critical  stress  observed  in  compression 

loading  of  a  simple  bar. 

N  —  Compressive  load  normal  to  plane  of  the 
interface  between  the  bearing  surfaces. 

T  =  sc  A,.  sc  =  Failure  shear  stress  of  a  bar  in  simple 

tension,  sc  —  0.5  rrc.  Where  ac  is  the 
tensile  stress  in  a  bar  undergoing  large 
plastic  flow  due  to  simple  tensile  loading. 

T  =  Traction  force  applied  along  the  plane  of  the 
interface  i.e.,  a  shearing  load. 


But  Pc  3  rrc  ,  sc  =  1/2  <rc  =  j  x  ^-pc 

T  scAr  1 

Since  f  =  —  then  f  =  — -  =  — . 

N  PcAr  6 

It  seems  that  a  simplified  analysis  leads  to  a  coefficient  of  friction  f  =  l/6.  Actual 
tests  of  friction  in  dry  metal  pairs  indicate  coefficients  ranging  from  f  =  0.5  to  f  — 
1.5.  Apparently  this  mathematical  model  based  on  simple  compression  needs 
refinement. 
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Combined  Stress  Model,  Junction  Growth  and  Surface  Contamination  Effects. 

If  a  shear  stress  is  indeed  applied  across  the  junction  area  the  combined  stress  pat¬ 
tern  is  now  such  that  the  pressure  component  available  to  sustain  load  N  is  actually 
reduced.  This  suggests  that  a  compensatory  increase  in  contact  area  will  take  place. 
This  area  growth  increases  the  required  shear  load  T  to  be  associated  with  sliding. 
This  alone  will  increase  the  coefficient  of  friction  from  the  earlier  estimate  of  l/6 
towards  values  more  like  those  experienced  in  practice.  Bowden  and  Tabor^  observed 
in  the  laboratory  that  two  similar  and  thoroughly  decontaminated  metals  will  exhibit 
coefficients  of  friction  f  >  100.  This  behavior  was  attributed  to  local  welding.  Sliding 
motion  creates  high  temperature  at  the  contacting  asperity  peaks.  This  speculation  is 
reinforced  by  laboratory  tests  during  which  sparks  (incandescent  temperatures)  were 
observed  under  a  microscope.  At  these  temperatures  pure  materials  will  spot  weld. 
The  welds  offer  shear  resistance  with  almost  negligible  call  for  a  squeeze  load  N,  hence 
the  observed  coefficient  of  friction  approaches  infinity  in  some  cases. 

As  in  normal  welding,  the  micro-welding  described  above  is  diminished  by  contam¬ 
inants  especially  oxides  and  moisture  which  are  ever  present  in  normal  sliding  condi¬ 
tions.  The  analytical  influence  of  less  than  clean  conditions  can  be  approximated  by 
calling  up  equation  EIII.l  which  is  a  reduced  plastic  failure  criterion  for  the  case  of 
combined  stresses. 


Sf  =  apparent  shear  failure  stress  of  contaminated  bearing  surface. 


Sf  =  c  \  c  =  cleanliness  coellicient  0  <  c  <  1 

T_  „  =  fL 

N  pAr  p 


Substituting  these  equivalents  into  Eqtn.  III.  1  leads  to 


f  = 


a(l— c2)* 


EIII.2 


If  the  surfaces  are  perfectly  clean,  c  =  1,  equation.  EIII.2  indicates  a  theoretical  f 
~  oo.  If  c  =  .85,  then  f  =  .51  and  if  c  —  .50.  f  =  .18.  Dry  metal  pairs  under  normal 
atmospherics  conditions  should  be  oxidized  sufficiently  to  expect  a  cleanliness  quotient 
c  of  about  .85  which  would  lead  from  EIII.2  to  a  f  expectation  of  about  0.5. 

Equation  EIII.2  implies  that  a  completely  oxidized  surface  (c  —  0)  should  enjoy  a 
coefficient  of  friction  f  =  0.  No  such  conclusion  can  be  gleaned  from  any  test  results. 
This  conflict  between  expectation  and  reality  is  more  easily  resolved  if  c  were  termed 
an  adhesion  coefficient.  What  seems  to  take  place  is  that  the  protective  oxide  being 
very  brittle  cracks  easily.  Under  the  high  local  pressures  of  the  asperities,  the  virgin 
metals  diffuse  through  the  fissures  in  the  cracked  oxide  film  to  meet  each  other  to  form 
welds.  Since  the  fissure  conduit  area  is  a  fraction  of  asperity  area  A'  the  apparent 
traction  strength  is  less  than  a  perfectly  decontaminated  surface  but  more  than  a 
completely  oxide  protected  surface. 

A  Refinement  to  the  Adhesion  Theory 

Edwards  and  Halling^’^  invoke  the  asperity  interlock  concept  and  interpose  it  on 
the  Bowden-Tabor  model.  In  this  case  the  supporting  mesa,  is  on  the  average  inclined 


at  some  mean  angle  0  to  the  horizontal  as  in  Fig.  III. 3c  as  contrasted  to  0  =  0  of  Fig. 
III. 3b. 


Moving  against  the  buttress,  defined  by  6,  increases  the  resistance  to  gross  motion. 
The  coefficient  of  friction  is  now  defined  as 

fox  "b  0  Ia\ 

fgjj  =  - —  —  coefficient  of  friction  by  Edwards  &  Hailing'  K  EIII.3 

(l-Offix) 

where  <f>  is  a  function  of  0  and  work  hardening  effects  on  the  junction  areas.  At  0 

=  0,  d>  =  0. 

fnT  =  - — -  =  coefficient  of  friction  calculated  by  the  Bowden-Tabor 

—  c  ) 

model  where  0  ~  0  ie,  equation  EIII.2. 

To  gather  a  sense  of  the  influence  of  asperity  interlocking  consider  a  case  of  0  —  10° 
and  fBx  —  -3  the  Edwards  Hailing  formula  yields  fEH  —  .45  or  a  50%  increase  in  the 
coefficient  of  friction  because  the  effect  of  interlocking. 

If  dissimilar  metals  are  run  dry  on  each  other  welding  is  resisted  somewhat  better. 
Phosphor  bronze,  brass  or  cast  iron  running  dry  against  steel  are  the  better  bearing 
combinations  with  dry  friction  coefficients  of  about  0.4.  While  this  figure  verges  on 
being  acceptable  for  guideways  the  problem  remains,  namely  dry  metal-metal  friction 
is  not  smooth  or  steady  enough  for  precision  systems.  There  are  however  non-metallic 
solids  which  not  only  have  a  somewhat  lower  to  very  low  values  of  f  and  behave  much 
more  smoothly  and  are  less  prone  to  seizure.  These  dry  solids  interposed  between  the 
metal  surfaces  add  lubricity  to  the  interface.  Resulting  coefficients  can  vary  between 
.03  and  .35.  Solid  lubricants  coine  in  inorganic  and  organic  forms.  The  two 
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classifications  have  entirely  different  mechanisms  to  arrive  at  their  dry  lubricity. 
Inorganic  Solid  Lubricants 

Fig.  III. lc  symbolizes  solid  dry  lubrication  as  a  process  whereby  flakes  of  a  solid 
slide  easily  on  each  other.  Graphite  is  the  original  solid  lubricant.  Its  lubricity  comes 
from  its  natural  laminar  structure.  It  is  a  crystalline  form  of  carbon  whose  atoms  are 
arranged  in  regular  layers.  While  these  layers  posses  high  bending  and  breaking 
strengths  the  bonds  between  them  are  very  weak.  The  evident  low  friction  is  a  man¬ 
ifestation  of  this  weakness  which  allows  the  layers  to  slide  on  each  other  rather  easily 
with  the  result  that  .05  <  f  <  0.15.  The  lower  limit  of  f  occurs  at  high  pressures  and 
the  upper  limit  at  low  pressures.  This  performance  is  maintained  up  to  200°C.  The 
sliding  is  smooth  and  steady.  Graphite  adheres  well  to  metal  surfaces  and  also  embeds 
well  into  carrier  materials  such  as  polymers. 

Graphite  has  some  idiosyncrasies  which  must  be  recognized  however.  Low  friction 
and  low  wear  are  realized  only  in  the  presence  of  water  vapor.  Graphite  does  not 
work  well  in  vacuum  and  dehydrated  conditions.  Furthermore  the  edges  of  the  gra¬ 
phite  crystal  planes  are  abrasive  hence  it  is  important  to  orient  the  crystals  parallel  to 
the  bearing  surface.  Fortunately  such  orientation  is  enhanced  by  the  lamellar  struc¬ 
ture  which  encourages  the  formation  of  flat  flakes  with  a  strong  tendency  to  stick  to 
flat  surfaces  in  the  desired  fashion.  The  problem  with  explicitly  introducing  a  form  of 
solid  lubricant  such  as  graphite  between  metal  surfaces  is  that  the  thin  film  is  eventu¬ 
ally  lost  by  wearing  away.  Fortunately,  unlike  wear  of  metals,  the  debris  from  solid 


■ 


lubricants  such  as  graphite  and  other  successful  solids  does  not  cause  rapid  wear  until 
the  surface  is  rubbed  down  to  raw  metal.  Nevertheless  the  life  of  the  film  is  an 
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engineering  concern  which  currently  is  assessed  in  terms  of  its  limiting  PV  factor 
(apparent  pressure  x  sliding  velocity).  This  term  if  augmented  with  an  appropriate 
coefficient  of  friction  suggests  power  dissipation  per  unit  area  of  the  bearing.  This 
connection  in  turn  suggests  that  the  local  temperature  generated  is  strongly  related  to 
the  level  of  PV  at  which  the  bearing  is  operating.  These  connections  are 
oversimplified  but  a  relationship  between  the  dynamic  strength  of  the  film  at  particu¬ 
lar  temperatures  is  roughly  integrated  into  limiting  PV  values  observed  in  life  tests. 
Ideally  more  complete  data  on  solid  lubricant  would  include  limiting  pressure  recom¬ 
mendations  and  thermal  conductivity  of  the  film.  Graphite,  having  a  superior  thermal 
conductivity,  enjoys  one  of  the  higher  allowable  PV  ratios  among  the  list  of  viable  dry 
solid  lubricants. 

Molybdenum  disulfide,  M0S2,  powder  is  the  other  important  inorganic  lubricant.  It 
is  normally  assumed  to  be  a  later  discovery  than  graphite  but  it  might  have  been  in 
use  as  long  as  graphite  but  not  recognized  as  such  since  it  was  confused  to  be  a  special 
graphite  form.  Again  low  friction  results  from  failure  along  the  cleavage  faces  this  time 
between  M0S2  layers.  This  weak  bond,  parallel  to  the  sliding  plane,  mated  with  great 
strength  normal  to  the  sliding  plane  is  the  basis  for  low  coefficients  of  friction  (.03  <  f 
<  .2  depending  on  load).  Unlike  graphite,  M„S2  does  not  need  water  vapor  or  oxida¬ 
tion  to  serve  its  purpose  hence  it  is  an  excellent  alternative  to  graphite  in  such 
environments.  MQS2  has  a  PV  limit  5  x  better  than  graphite.  Furthermore  M0St  has 
even  better  adherence  to  metals  than  graphite  and  a  load  carrying  capacity  28  x 
better  than  graphite.  Graphite  does  have  a  distinct  superiority  if  electrical  and  ther¬ 
mal  conductivity  is  desirable.  Except  for  the  latter  operating  condition,  M„S2  films 
will  be  superior  to  graphite  films  for  dry  precision  guideways. 
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Molybdenum  films  can  be  deposited  on  metal  surfaces  by  spraying  a  mist  of  M0S2 
in  a  volatile  fluid  carrier.  After  the  carrier  has  evaporated  from  the  surface  M0S2 
remains.  The  surface  is  then  rubbed  with  an  appropriate  cloth  burnishing  the  M„S2 
into  the  metal.  The  resulting  film  has  a  very  low  coefficient  of  friction  and  a  thickness 
of  only  0.1//m  which  insures  dimensional  control  so  necessary  to  a  precision  guideway. 
Increased  durability  can  be  achieved  by  increasing  the  film  thickness  up  to  10/im 
through  repetitions  of  the  process  described  above.  The  larger  film  thickness  does 
detract  from  the  quality  of  a  guideway  precision.  Sputtering  is  an  alternative  process 
for  laying  down  very  thin  yet  very  durable  films  for  precision  guides.  Another  con¬ 
venient  method  for  laying  down  a  film  is  to  spray  M0S2  in  a  resin  carrier  onto  the 
metallic  surface.  The  carrier  again  is  allowed  to  evaporate  and  then  the  film  is  cured 
at  a  specific  temperature  to  give  a  hard  varnish  like  finish.  Such  a  film  has  a  very 
good  life  but  the  resulting  thickness  is  many  times  greater  than  that  converted  with  a 
burnishing  or  sputtering  method.  Again  the  thickness  of  resin  bonded  films  might  not 
be  compatible  with  precision  guideways. 

Organic  Solid  Lubricants 

The  other  approach  to  dry  film  lubrication  is  by  way  of  polymers.  The  two  most 
successful  are  Nylon  and  PTFE  (polytetrafluoroethylene)  commercially  known  as 
Teflon,  Fluon  etc..  PTFE  is  the  more  interesting  polymer  because  its  coefficient  of 
friction  can  be  as  low  as  .03.  Its  properties  and  mechanics  of  sliding  are  unlike  other 
materials.  A  very  low  coefficient  of  friction,  chemical  inertness  and  a  high  melting 
point  are  the  characteristics  which  have  made  PTFE  such  a  desirable  and  popular 
material.  But  there  are  prominent  negative  characteristics  that  have  been  difficult 


barriers  to  bearing  applications.  Its  normal  form  is  acquired  by  moulding  but  it  is  not 
easily  machined  because  of  excessive  thermal  distortion  created  by  machining,  l'inalty 
the  very  nature  of  its  low  friction  is  due  to  its  poor  mechanical  strength.  T  he  attrac¬ 
tiveness  of  the  low  coefficient  of  friction  and  its  stability  at  high  and  low  temperatures 
(300°C  to  -200°C)  has  been  a  motivator  for  development  of  methods  to  overcome  the 
stubbornly  negative  qualities. 

For  reasons  that  are  not  clear,  despite  chemical  inertness  and  poor  mechanical 
strength  PTFE  forms  good  smooth  strongly  adhering  films  to  metal  surfaces.  Good 
smooth  films  can  be  produced  by  a  transfer  process.  Solid  PTFE  is  made  to  slide 
under  pressure  on  a  clean,  dry  metal  surface  and  so  transfer  itself  to  the  metal.  But 
the  film  so  created  has  insufficient  wear  resistance  for  bearings. 

PTFE  in  solid  form  in  addition  to  poor  mechanical  strength  also  suffers  from  creep 
under  load  and  poor  wear  resistance.  To  help  circumvent  these  deficiencies,  PTFE  is 
currently  impregnated  into  a  reinforcing  matrix  of  such  materials  as  glass  fibers,  car¬ 
bon  fibers,  M0S2  and  porous  lead-bronze.  The  reinforcing  carriers  add  strength  and 
increase  the  modulus  which  in  turn  reduces  the  deformation.  Finally  wear  rate  is 
improved  10,000  fold.  But  the  coefficient  of  friction  of  the  composite  is  raised  to 
approximately  0.25  from  the  promise  of  f  =  .03.  Several  commercially  available  com¬ 
posites  offer  PTFE  incorporated  in  the  pores  of  a  sintered  bronze  strip  which  itself  is 
bonded  to  a  steel  backing  strip.  Another  widely  used  commercial  composite  with  PTFE 
is  in  the  form  of  fibers  interwoven  with  glass  fiber.  The  weave  is  impregnated  with  a 
phenolic  resin  the  whole  bonded  to  a  steel  back  strip.  These  composites  can  sustain 
maximum  PV  values  of  100,000  psi  x  ft/m.  The  stability  of  PTFE  composite  reinforce¬ 
ments  have  good  prospects  for  guides  especially  those  which  must  operate  in  non-fluid 


and  vacuum  environments. 


Nylon  (Dupont)  is  another  polymer  with  considerable  utility  in  mechanical  systems. 
Its  mechanical  strength  and  resistance  to  creep  is  superior  to  that  of  PTFE. 
Nevertheless  like  PTFE  to  acquire  sufficient  strength  and  stability,  Nylon  must  also 
composited  for  satisfactory  guideways.  Because  of  growth  due  to  moisture  absorption 
Nylon  at  present  is  probably  inadequate  for  precision  guides. 

Delrin  (Dupont)  is  an  acetal  resin  filled  with  PTFE.  The  acetal  resin  lends  its  high 
strength,  low  creep  and  stiffness  to  the  composite  while  the  PTFE  provides  the  low 
friction.  The  static  and  kinetic  coefficients  of  Delrin  500  AF  are  .08  and  .14  respec¬ 
tively. 

Boundary  Lubrication,  Fig,  III. 2b 

The  coefficient  of  friction  and  wear  life  of  dry  bearing  surfaces  can  be  improved 
considerably  by  wetting  the  interface  of  the  solids  with  an  oily  fluid.  The  mechanics 
of  the  easy  sliding  is  identified  as  Boundary  Lubrication.  Although  the  introduction  of 
a  fluid  might  lead  one  to  expect  the  phenomenon  of  friction  to  be  vastly  different  than 
that  described  in  the  previous  sections  the  fact  is  that  the  principles  are  remarkably 
similar  that  is: 

a)  A  contaminant  impairing  the  welding  of  asperity  junctions  is  introduced. 

b)  The  contaminant  reduces  the  shear  strength  of  the  asperity  junctions. 

c)  Finally  the  contaminant  increases  the  ratio  of  crushing  strength  to  shear 
strength  of  the  thin  film  involved. 

The  oily  film  involved  in  boundary  lubrication  goes  beyond  the  obvious  role  of  a 
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weld  impairing  contamimant.  Hardy^  provided  the  earliest  fuller  explanation  of  the 
part  that  certain  liquids  play  in  boundary  lubrication.  Bowden  and  Tabor^2'  built  on 
Hardy’s  concepts  to  develop  a  model  of  the  mechanics  of  boundary  lubrication.  Hardy 
observed  that  successful  fluid  lubricants  possess  a  long  chain  highly  polarized  molecu¬ 
lar  structure.  The  polar  head  of  such  molecules  induce  a  strong  physical  adsorption  to 
the  bearing  surfaces.  The  long  chain  molecules  standing  on  end  with  heads  adhering 
to  the  metal  surfaces  appeared  somewhat  like  brush  bristles  which  would  separate  the 
surfaces  yet  allow  the  bearing  surfaces  to  move  lateraly  on  each  other  by  elastically 
flexing  to  the  displacement  and  snapping  back  into  a  supporting  mode.  The  interfer¬ 
ing  asperities  concept  is  not  present  in  Hardy’s  model.  The  inconsistences  between  the 
mechanistic  model  and  observed  behavior  caused  Bowden  and  Tabor  to  modify 
Hardy’s  model. 

The  B  and  T  model  is  depicted  in  Fig.  III.3b.  The  concept  of  squeeze  load  support¬ 
ing  asperities  is  retained  but  now  this  load  is  shared  by  Hardy’s  bristle  like  oily 
molecules  filling  the  voids  between  asperities.  The  real  area  of  support  is  increased 
from  Ar  described  in  the  early  sections  to  A  of  Fig.  III. 4.  The  asperity  function  area  is 
a  fraction,  <*,  of  A.  The  area  supported  by  the  molecules  introduced  by  the  oil  is  the 
remaining  fraction  1— <».  Because  of  this  cooperative  support  of  the  load  the  asperity 
junction  area  is  much  reduced  and  of  course  the  presence  of  the  fluid  essentially  elim¬ 
inates  asperity  welding. 


Fig.  IIL4  Long  Molecule  Support  of  Sliding  Bodies. 


It  has  been  known  a  long  time  that  animal  fats  and  vegetable  oils  make  good  boun¬ 
dary  layer  lubricants  while  mineral  oils  and  alcohols  are  poor  lubricants.  Some  key 
ingredients  responsible  for  good  boundary  lubrication  are  stearic  acid  (animal  fats)  or 
palmitic  acid  (cotton  seed  and  palm  oil)  or  oleic  acid  (olive  oil).  Such  ingredients  are 
instrumental  in  creating  strong  adherence  of  the  fluid  molecules  to  the  metal  surface. 
It  is  poor  adherence  to  the  working  surfaces  that  marks  pure  mineral  oils  as  poor 
lubricants.  On  the  other  hand  if  only  1  to  5%  by  volume  of  stearic  acid  is  added  to 
mineral  oil  it  is  converted  into  a  good  lubricant.  Strangely  stearic  acid  and  the  other 
favorable  acids  by  themselves  are  not  good  lubricants  but  rather  play  a  chemical  role. 
The  presence  of  these  acids  causes  a  chemical  reaction  with  the  oxides  of  some  of  the 
more  reactive  bearing  materials.  The  result  of  this  reaction  is  the  creation  of  a  new 
substance,  a  metallic  soap  film  chemically  bound  to  the  bearing  metal.  Copper,  zinc, 
tin,  etc.  are  good  reactive  materials  for  boundary  lubrication.  Boundary  lubrication  is 
poor  on  the  noble  metals  because  their  chemical  inertness  giving  witness  to  the  obser¬ 
vation  that  it  is  the  chemical  attack  of  the  favorable  oils  that  makes  for  good  boun¬ 
dary  lubrication.  Typically  reactive  metals  will  yield  a  coefficient  of  friction  f  =  0.1 
whereas  with  unreactive  metals  .3  <  f  <  .7. 
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The  metallic  soap  under  high  pressure  becomes  a  solid  lubricant  sharing  the  load 
with  the  asperities.  Under  favorable  boundary  lubrication  conditions  the  asperity 
junction  area  of  Fig.  III.4  is  very  small  compared  to  the  area  supported  by  the  metallic 
soap  (a  is  small).  Therefore  the  friction  is  predominantly  that  of  the  metallic  soap.  If 
wear  is  defined  as  the  breaking  off  of  the  tips  of  the  asperity  junctions,  the  reduced 
junction  area  of  Fig.  III.4  makes  for  smaller  wear  particles.  In  this  way  the  support 
contributed  by  the  metallic  soap  results  in  lower  ware  rates. 

Boundary  lubrication  is  normally  discussed  in  terms  of  a  thin  oil  film  between  slid¬ 
ing  metal  surfaces  as  was  the  case  above.  Applications  of  oil  to  organic  solids  is 
beneficial  although  the  chemical  generation  of  soaps  is  more  problematical  in  such 
cases.  Denton^  studied  the  efficacy  of  two  commercial  types  of  PTFE  composites 
operating  with  liquid  lubricants.  The  composites  were  in  the  form  of  .060  inch  skived 
tapes  intended  as  bearing  surface  overlays  for  linear  slidcways  of  machine  tools.  Slid¬ 
ing  tests  were  conducted  with  bearing  loads  generating  an  apparent  pressure  of  200 
psi.  Kinetic  coefficients  of  friction  and  wear  were  observed  under  pv  =  3333  psi  x 
ft/m  operating  conditions.  Static  coefficients  of  friction  were  measured  at  the  termi¬ 
nation  of  each  sliding  run.  By  way  of  a  control  identical  tests  were  conducted  with 
lubricated  cast  iron  slideways  surfaces.  TABLE  111.1  is  or  an  abstract  of  test  results 
pertinent  to  this  discussion.  Friction  coefficients  are  given  for  sliding  on  freshly  ground 
hardened  steel  surfaces.  Initially  the  readings  are  given  for  a  very  short  run  (D  < 
1.0").  Run-in  operations  are  represented  by  the  results  observed  after  a  sliding  pas¬ 
sage  of  10s  inches.  The  short  run  tests  with  the  two  PTFE  composites  types  show 
kinetic  coefficients  of  friction  fk  .035  and  fk  -  .031.  These  values  are  essentially 
equal  to  the  kinetic  coefficient  for  bulk  PTFE  on  steel.  On  the  other  hand  after  a 
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run-in  period,  D  =  105  inches,  the  friction  coefficients  have  increased  substantially  to 
fk  =  .131  and  fk  =  .097  for  the  two  PTFE  types.  This  pattern  of  results  suggests 
that  at  the  on  let  of  run-in  the  contacting  asperities  are  predominantly  PTFE  against 
steel  but  as  the  run-in  progresses  the  soft  PTFE  asperities  are  worn  down  exposing 
more  of  the  filler  to  the  hardened  steel  mating  surface.  The  coefficient  of  friction  at 
the  end  of  the  run-in  period  is  the  resultant  of  PTFE  and  filler  contributions  to  fric¬ 
tion.  Near  the  end  of  the  run-in  (D  —  105  inches)  the  coefficient  of  friction  for  the 
mineral  filler  composite  has  risen  to  fk  --  .131  while  the  bronze  filler  composite  has 
risen  to  fk  =  .097  or  a  coefficient  of  friction  35%  less  than  that  of  the  mineral  filler 
model.  Apparently,  since  bronze  is  favorably  reactive  to  the  lubricating  oil,  good 
boundary  lubrication  occurs  on  and  around  the  bronze  asperities.  On  the  other  hand 
the  mineral  asperities  not  being  reactive  to  the  lubricant  do  not  develop  beneficial 
boundary  lubrication.  It  could  be  said  then  that  the  presence  of  a  lubricant  does  very 
little  to  reduce  friction  of  PTFE  composites  unless  the  filler  is  a  ‘good’  bearing  metal 
that  is  one  that  is  associated  with  effective  boundary  lubrication.  It  should  be  remem¬ 
bered  however  that  a  liquid  lubricant  also  serves  as  a  coolant  and  anti-corrodent  so 
that  at  higher  "Tiding  velocities  the  PTFE  mineral  composite  could  enjoy  some  benefit 
from  these  secondary  roles  of  lubricants  in  bearing  systems. 


The  cast-iron  tests  show  performance  to  be  much  more  independent  of  length  of 
run-in  distance.  This  could  be  because  the  surface  of  cast-iron  can  be  ground 
smoother  than  is  possible  for  PTFE  composites.  In  a  way  cast  iron  is  also  a  composite, 


being  a  distribution  of  graphite  in  the  body  of  iron.  Wearing  down  of  a  cast  iron  sur¬ 


face  merely  exposes  new  graphite  spots  so  that  coefficients  of  friction  remain  fairly 


constant  from  the  onset  of  sliding  through  run-in,  as  indicated  in  TABLE  III.  I. 
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Meanwhile  the  exposed  iron  being  reactive  to  the  oil  lubricant,  boundary  lubrication 
prevails  as  well. 

The  ratio  of  fs/fk  is  a  rough  indicator  of  stick-slip  (chatter)  tendencies  of  a  pair  of 
bearing  surfaces.  If  this  ratio,  called  the  "stick-slip  coefficient"^,  is  equal  to  or  greater 
than  unity  there  is  risk  of  chatter.  From  TABLE  III.  1  after  run-in  (D  =  105in.)  the 
PTFE  mineral  composite  is  running  dangerously  close  to  chattering,  while  the  PTFE- 
bronze  composite  and  the  cast  iron  have  almost  identical  safe  fs/fk  ratios.  Once  again 
the  existance  of  boundary  lubrication  shows  its  beneficial  hand. 

TABLE  III.l 


Surfaces  and  Lubricants 

D 

1* 

2* 

3* 

f. 

.027 

.027 

.110 

0-1  in 

fk 

.035 

.031 

.126 

f.Ak 

.77 

.87 

.87 

fs 

.128 

.084 

.112 

lO^in 

fk 

.131 

.097 

.127 

f./fk 

.98 

.87 

.88 

lO^in 

K  wet 

5804 

1757 

1916 

K  dry 

11220 

30400 

... 

105in 

K  wet 

30 

39 

13 

K  dry 

37 

138 

... 
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1*  Mineral  filled  PTFE,  Compound  Way  oil  againsl  steel. 

2*  Bronze  filled  PTFE,  Compound  Way  oil  against  steel. 

3*  Cast  iron,  Compound  Way  oil  against  steel. 

fg:  static  coeff.  of  friction 

f^:  kinetic  coefficient  of  friction 

fs/fk:  stick  slip  coefficient. 

K  ( 1 02):  (in.3  min. /ft.  lb.  hr.)  x  10-10  after  102  units  of  travel. 

K  (105):  (in.3  min. /ft.  hr.)  x  10-10  after  10s  units  of  travel. 

D:  Sliding  distance  traveled  at  observed  data. 

Summary 

Table  III.  1  demonstrates  that  PTFE  composite  surfaces  could  work  on  ordinary 
machine  tool  slideways  but  the  wear  rate  and  geometrical  instability,  due  to  moisture 
absorption,  makes  them  poor  candidates  at  this  point  of  development  for  precision 
guideways.  Dimensional  growth  of  almost  10%  due  to  moisture  absorption  as  not  com¬ 
patible  with  precision  tolerences  of  less  than  l//m.  A  wear  rate  3  -  4  x  more  rapid 
than  a  lubricated  cast  iron  surface  can  be  accommodated  by  frequent  replacement  of 
the  composite  tapes.  While  this  is  an  acceptable  procedure  for  reconditioning  an  ordi¬ 
nary  machine  tool  it  is  not  so  for  precision  surfaces  since  each  replacement  will  intro¬ 
duce  errors  well  beyond  that  tolerated  for  precision  systems.  The  same  assesment 
could  be  made  of  M„S2  composites.  The  synthetics  are  at  present  not  likely  substi¬ 
tutes  for  lubricated  cast  iron  precision  slideways. 
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G'r  p.n  a  lapped  guideway  with  a  residual  error  in  straightness  described  as  a  .001 


mm  camber  in  a  sector  400  mm  long.  Assume  that  the  curved  path  is  part  of  a  circle 
then  the  given  path  geometry  can  be  expressed  as  a  radius  of  curvature  R  and  slope 


M  Fig.  AIII.l. 


Thermal  Distortion,  1  C 
Variation  from  Top  to  Bottom 

Fig.  AIII.l  Geometry  of  Curved  Path 


Camber  d  =  .001  mm  where  d  =  —  R  A#2 

2 


Sector  length  L  =  400  mm  and  —  =  RA# 

2 


Slope  A B  <  .01  rad. 


Hence  R  =  —  =  — — —  =  2  x  107  mm 
8d  8x.001 


A<9  =  — -  =  - -  =  20  //  rad 

2x2x1 07 


Ther  ral-Bowing  of  a  Straight  Guidewat 


G  .en  a  guide  rail  in  the  form  of  a  prismatic  bar  of  depth  h.  Along  the  length  of 
the  rail  there  is  a  sector  L  units  long  wherein  the  difference  in  temperature  between 
the  upper  and  lower  surfaces  is  AT.  The  strains  developed  by  this  temperature  gra- 


V'} 

’ 
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dient  will  cause  the  rail  to  bow  as  in  Fig.  AIII.I  to  a  curvature  of  radius  R'  and  a  slope 
A#  rt  each  end  of  the  thermal  sector.  The  thermal  strains  e  at  the  upper  and  bottom 
surfaces  are 


<*AT 


e  =- 


Furthermore  e  = 


a  —  coefficient  of  thermal  expansion  for  the  rail, 
h  h 


or  R'  = 


2R'  oAT 

Given  R',  the  corresponding  camber  d'  and  slope  A ff  are  determined  by  the 
geometrical  transformations  given  earlier  i.e. 


T  *  1  T  * 

d'  1=  —  «AT 

8R'  8  h 

A("  =  =  —  <»AT 

2R'  2h 

For  purposes  of  illustration  consider  a  cast  iron  guideway  with  the  following 
specifications: 

o  =  9  x  10“6,  L  =  400  mm,  AT  =  1°C 


It  follows  then  that: 


d'  = 


_L  4002  a  (9  x  10  6)  x  1.0 
8  X  100 


=  1.8  // m 


A&  _  — 152 —  m  x  10  8)  x  1.0  =  1.8  //  rad 
2  x  100  v  ' 
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IV.  THICK-FLUID  FILM  LUBRICATION 

The  preceding  chapter  dealt  with  two  bodies  wherein  sliding  on  each  other  was 
eased  by  a  solid  protective  film  as  little  as  .005  microns  thin.  The  friction  and  wear  of 
these  thin  films  are  not  negligible  but  in  the  light  of  the  simplicity  of  such  bearing 
systems  they  do  provide  effective  lubrication  for  many  commercial  applications.  On 
the  other  hand  aside  from  metal  to  metal  boundary  lubrication  wear,  friction  and 
dimensional  instability  of  these  thin  films  are  not  acceptable  for  precision  systems 
especially  at  the  higher  rubbing  speeds.  An  alternative  strategy  for  bearing  systems  is 
to  create  a  fluid  film  which  is  pressurized  sufficiently  to  push  the  bearing  surfaces 
apart  so  that  surface  asperities  do  not  contact  or  interlock  on  each  other.  When  this 
separation  is  accomplished  wear  is  completely  eliminated  and  friction  is  reduced  to 
viscous  shearing  of  the  fluid  film.  Such  attractive  behavior  is  indeed  realized  in 
externally  and  self  pressurized  (EX.P  &  S.P.)  bearings  wherein  the  friction  is  less  than 
10%  of  the  best  thin  film  experience.  While  the  fluid  space  separating  the  metal 
surfaces  in  these  bearings  is  very  small  (less  than  2.5  /y m  in  some  instances)  it  is 


nevertheless  hundreds  of  times  thicker  than  that  of  solid  films  discussed  in  section  III. 


Hence  by  contrast  these  bearings  are  referred  to  as  thick  film  bearings.  Alternatively 
«un<  <  the  frid  ion  experienced  is  so  much  less,  this  form  of  lubrication  is  sometimes 
referred  l<>  as  perfect  lubrication.  The  friction  and  wear  characteristics  of  the  thin 
film*  of  tertion  Ill  were  a  consequence  of  the  rheological  (mechanical)  properties  of  the 
*'  '  i£*  ii<  rated.  On  the  other  hand  the  viscosity  of  the  fluid  is  the  property  dictating 


ti  •  *  1  i  ><!  film  lubrication. 
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t 

|  EX.P  and  S.P.  lubrication  is  essentially  a  process  whereby  fluid  inertia  and  viscos¬ 

ity  forces  generate  hydraulic  pressures  capable  of  lifting  one  surface  off  its  mate.  The 
analytical  mechanics  describing  this  action  is  based  on  the  fundamental  Navier-Stokes 

i 

Equations  of  Appendix  ATV.l.  A  completely  rigorous  solution  of  the  N-S  equations  is  a 
formidable  challenge  even  with  the  fastest  computers.  Fortunately  there  are  reason¬ 
ably  realistic  engineering  approximations  appropriate  to  bearings  which  can  reduce  the 
N-S  equations  to  tractable  forms.  Stripped  of  the  less  effective  terms  as  is  done  in 
Appendix  AIV.l  the  remaining  key  elements  clarify  the  physics  of  pressure  generation 
and  distribution.  These  simplifications  lead  to  useful  quantitative  determinations. 


Appendix  AIV.l  demonstrates  the  process  of  reducing  the  N-S  equations  to  the  par¬ 
tial  differential  equations  of  EX.P.  or  S.P.  lubrication.  The  fundamental  simplifying 
assumption  stems  from  the  observation  that  the  film  thickness  of  practical  fluid  film 
bearings  is  infinitesial  compared  to  the  other  two  flow  path  dimensions.  It  is  therefore 
further  assumed  that  pressure,  viscosity,  and  density  variations  across  the  film  are 
negligible  compared  to  these  variations  along  the  other  two  flow  paths.  The  slit  like 
geometry  of  the  flow  path  also  encourages  the  assumption  that  the  Reynolds  number 


in  most  practice  will  be  sub-critical.  That  is  the  inertia  forces  of  the  N-S  equations 
are  negligible  compared  to  the  dissipative  (viscous)  forces.  These  initial  relaxations  of 
the  complete  N-S  equations  lead  to  lubrication  equation  HIV.l.  This  equation  will  be 
referred  to  here  as  the  generalized  Reynolds  equation. 


-  (—  I'  )  +  ~~  (—  /'  Tj^)  =  I2(r.h.s.) 
(h c  it  dx  <>y  it  <•  y 


EIV.l 
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where:  h(x,y)  =  fluid  film  thickness 

/i(x,y)  —  fluid  viscosity 
(>{x,y)  —  mass  density  of  fluid 
p(x,y)  =  fluid  pressure  of  the  film 


r.h.s.  =  p 


U1(f),-Ull(f)ll+V1(f,1 


+ 


)  +  -^(phw  )  + 

<rx  (jy 

-  (V  +  +„) 

where:  «  „d  W  =  ^  +  and  <*(''/»)  + 

2  2 


Subscripts  II,  I,  and  0  refer  to  points  on  the  upper  surface, 

lower  surface  and  stationary  reference  point. 


d'lj  f  and  d>i  f  —  fluid  mass  flux  rate  into  interface  clearance  film  through 
surface  II  and  I  respectively. 

Furthermore  additional  reductions  are  allowed  if  the  following  observations  are  also 
applicable 


1.  Bearing  surfaces  are  not  porous.  That  is  there  is  no  fluid  flow  through  the 
bearing  surfaces  d’n  f  =  r  =  0. 


ih) 


2.  Fluid  is  incompressible  that  is  density  is  invariant  (-7-  =  0). 


3. 


Viscosity  is  constant  throughout  the  film. 
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4.  Only  surface  II  is  moving  and  that  only  along  the  x  axis. 


5.  Surfaces  do  not  pulsate  {-7^-  =  0  = 

v  (h  Di 


6.  Surfaces  are  not  stretchable, 


dUn 

dx 


=  0. 


7.  Surface  I  lies  in  the  x,  y  plane  ie  — )i  =  —  0. 

(h.  ay 


These  further  reductions  to  the  N.S.  equations  lead  to  the  familiar  lubrication 
equation  EIV.2a. 


(h3 


D_ 

Dx  y  Dx 
where:  U 


9£L+  A 


(h3^-)=- 
Dy  Dy 


*  TT 

6  fJ  V  —  ... 
OX 


C(II/0) 


EIV.2a 


Even  before  equation  EIV.l  is  actually  solved  for  p(x,y)  it  does  convey  the  roles 
played  by  the  parameters  of  the  physical  model.  The  right  hand  side  (r.h.s.)  identifies 
the  terms  responsible  for  the  fluid  flow  (pumping)  sources  within  the  clearance  between 
the  bearing  surfaces.  The  left  hand  side  (l.h.s.)  of  EIV.l  not  only  contains  the  pres¬ 
sures  one  would  seek  from  a  solution  but  also  the  resistivity  to  the  flow  through  the 
clearance  h(x,y).  Referring  to  the  rudimentary  model  and  its  equation  EIV.2a  the 

r.h.s.  indicates  that  for  self  generation  of  flows  to  occur,  the  product  •  ~~  must 

exist  in  most  regions  of  the  film  domain.  That  is  one  surface  must  have  a  sliding  velo¬ 
city  and  the  clearance  h  must  vary  in  places  i.e.  a  hydraulic  wedge  or  step  must 
occupy  some  portion  of  the  film  domain.  In  the  case  of  a  linear  slideway  X(u/0)  does 
occur.  If  one  surface  were  allowed  to  rock  relative  to  its  mate  Dh/Dx  would  also  occur 
but  such  cyclic  rocking  is  not  acceptable  for  precision  systems.  In  precision  slideways 
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— —  =  0  except  for  errors  in  the  guiding  surfaces  hence  the  r.h.s.  of  EIV.2a  is  essen- 

(>X 

tially  zero  and  the  lubrication  equation  becomes  EIV.2b 


0'3J£)  -  o ... 

<*y 


EfV.2b 


If  atmospheric  gauge  pressure,  p  =  0,  prevails  as  the  sole  boundary  condition  to 
EIV.2b  there  will  be  no  pressure  generated  in  the  film.  Thick  film  lubrication  therefore 
fails  to  develop  in  this  model.  The  bearing  system  retreats  to  the  thin  film  of  boun¬ 
dary  lubrication  described  in  section  III. 


It  is  sometimes  argued  that  hand  scraping  of  linear  guideways  creates  pockets  with 

sloped  sides  such  that  local  =/=  0  and  therefore  the  r.h.s.  of  EIV.2a  could  exist  in 

dx 

spots  throughout  the  film  domain.  But  pockets  due  to  scraping  or  error  undulations  in 
the  surfaces  produce  a  surface  profile  of  odd  symmetry.  The  resulting  pressures  gen¬ 
erated  will  also  be  in  odd  symmetry  throughout  the  film  domain.  Hence  the  sum  total 
of  the  localized  self  generated  lifts  is  zero.  This  neutralizing  symmetry  could  be  des¬ 
troyed  if  the  sliding  velocity  X(n/0).  were  sufficiently  high  to  cause  cavitation  on  the 
trailing  side  of  each  pocket.  The  velocity  needed  to  incur  cavitation  however  is  well 
beyond  the  sliding  velocities  of  linear  guides.  Finally  since  stop-go  motion  is  the  usual 
mode  of  action  for  linear  guides,  the  sliding  velocity  near  each  end  of  travel  decays 

causing  the  r.h.s.  of  the  Reynolds  equation  to  decay  as  well  even  if  an  average  did 

dx 

exist.  In  short  it  could  be  said  that  self  pressurized  (S.P.)  lifting  of  the  surfaces  is  not 
a  realistic  expectation  for  precision  linear  guideways.  Unless  external  pressurization 
(EX.P.)  is  provided  precision  linear  guideways  will  operate  on  only  the  thin  film  of 
boundary  lubrication. 
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Externally  Pressurized  Guideway  (The  Hydrostatic  Bearing) 

If  the  r.h.s.  of  EIV.2a  is  negligible,  as  argued  above  with  respect  to  linear  precision 
guideways,  this  equation  can  be  reduced  to  the  homogeneous  form  ElV.2b.  If  further 

one  considers  a  flat  linear  guideway  —  0)  this  equation  can  be 

chc  ()y  (h 

reduced  further  to  the  Laplacian  form  EIV.3. 


ax2 


a2P 

+  Tj  -  0 

<hr 


EIV.3 


The  determination  of  the  pressure  p(x,y)  distribution  from  EIV.3  depends  on  the 
x,y  plan  view  of  the  domain  I)  occupied  by  the  working  fluid  and  the  pressures,  pB, 
existing  at  the  boundaries  of  the  domain.  A  trivial  solution  of  EIV.3  occurs  when  pB  is 
constant  all  along  the  boundaries.  Then  by  observation  p(x,y)  =  pB  —  a  constant.  If 
only  ambient  gauge  pressure  {pB  ==  0)  acts  on  the  boundaries  domain  pressure  p(x,y) 
will  also  disappear.  Hence  there  is  no  lift  force  as  was  already  suggested  in  the  earlier 
discussion  of  EIV.2a.  To  generate  a  working  pressure  inside  the  domain  it  is  necessary 
that  the  boundary  pressures  impressed  vary  around  the  boundaries  of  the  lifting 
domain.  Such  is  the  condition  which  prevails  in  KX.l*  Hat  bearings.  Fig.  IV.  1  is  an 
effective  model  of  a  basic  hydrostatic  linear  guideway.  The  upper  sliding  member  II  is 
a  precision  ground  surface  with  a  series  of  pockets  machined  into  the  surface  Fig. 

IV.la.b.  The  pressure  distribution  working  on  the  sliding  member  11  at  sample  sections 
A-A  and  B-B  is  also  illustrated.  The  integral  of  the  pressures  at  such  sections  becomes 
the  generated  lift  W.  If  laminar  flow  is  assumed  for  the  model  of  Fig.  IVl.a.b  there  is 

an  analogous  electrical  resistance  net  work  (see  Appendix  1  for  analogs).  Hence  the  ] 

< 

established  iconics  and  terminology  of  electrical  net  works  can  be  transferred  to  the  i 


Mechanics  of  an  Elementary  Hydrostatic  Bearing. 

The  pump  between  terminals  0  (sump)  and  1  (pump  outlet)  of  the  hydraulic  system 
Fig  IVl.a,  raises  the  pressure  by  p^.  The  pressurized  oil  continues  through  the  mani¬ 
fold  from  1  to  2  with  a  negligible  pressure  drop  hence  p2/0  —  Pi/o-  Between  stations  2 
and  3  there  is  a  fluidic  resistor  causing  the  pressure  drop  p2/3.  The  fluid  continues 
into  the  pocket  identified  by  the  boundary  contour  4.  As  will  be  explained  later  the 
pocket  pressure  drop  p3/4  cs  0.  The  fluid  continues  from  boundary  4  through  the 
bearing  clearance  domain  to  the  exit  boundary  5.  From  5  the  fluid  continues  by  grav¬ 
ity  feed  back  to  the  sump  0.  The  flow  through  the  domain  bounded  by  contours  4  and 
5  is  through  a  passage  with  parallel  flat  surfaces  separated  by  a  film  thickness  h  of  the 
order  of  50  //m.  This  surface  domain  throughout  which  the  clearance  h  defines  the 
lubricating  film  is  alternatively  referred  to  as  the  clearance  zone  or  as  the  land  of  the 
bearing.  The  viscous  resistance  to  this  flow  and  the  flow  from  2  to  3  an*  the  control¬ 
ling  resistances  (R4_5  and  R2-3)  1°  the  total  flow. 

The  starting  point  to  the  quantitative  analysis  is  the  determination  of  the  lift  force 
W  generated  by  the  surface  with  boundaries  defined  by  the  contours  4  and  5.  The 
pocket  areas  yield  a  lift  Wp  --  p3/0  x  pocket  area.  The  pressure  at  terminal  0  is  the 
ambient  pressure.  Further  it  may  be  assumed  that  there  is  a  negligible  pressure  drop 
from  3  to  4  across  the  pocket.  The  land  domain  (4-5)  contributes  a  lift  W,  based  on 
the  integration  of  pressure  over  the  area  of  this  region.  Solution  of  E1V.3  will  provide 
p(x,y)  for  this  integration.  Solving  this  equation  for  a  domain  even  as  simple  as  that 
shown  in  Fig.  IV.  1  however  involves  complex  techniques.  For  exposition’s  sake  con¬ 


sider  initially  the  simpler  hydrostatic  bearing  of  Fig.  IV. 2  rather  than  the  model  of 


il  »l 


P  (r)  _  In  (b/r) 
p4/s  Xn(b/a) 


(c) 


IY.2  Circular  Hydrostatic  Bearing 
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The  Laplace  partial  differential  equation  1CIV.3  prevails  for  the  fields  of  Fig.  IV,  2 
but  polar  coordinates  are  more  convenient  for  the  case  of  a  circular  domain.  A 
transformation  such  as  described  in  the  complex  variable  reference^  yields  the 
equivalent  harmonic  equation  EIV.4a 


,»£lL  +  r4t  + 

ih2  ffr  (>02 


ErV.4a 

From  symmetry  of  the  pressure  map  about  axis  there  will  be  no  tangential  varia¬ 
tion  in  p  i.e. 

&SL  _  _  o 

no  ae2 

Hence  the  coordinate  0  disappears  as  a  variable.  The  partial  differential  equation 
above  now  becomes  the  ordinary  differential  equation  EIV.4b 

r2  +  r  =  0  ElV.4b 

dr2  dr 

From  Appendix  A.IV.2  the  solution  to  EIV.4b  for  the  specific  case  of  Fig.  FV.2  is 


p(r)  _  Ln(b/r) 
P4/o  Ln(b/a) 


ElV.Sa 


Lifting  Capacity  of  Hydrostatic  Bearings 


The  pressure  solution  integrated  over  the  land  area  4-5  yields  the  lift  .  The  lift 
contributed  by  the  pocket  area  is  Wp  =  7ra2  p3/0.  The  total  lift  W  —  +  Wp  is 

given  by  EIV5.b  (Appendix  AIV.2). 


W  =  p3/0  •  *b2  •  Kw  where  Kw  = 

L„(b/a)‘ 


EIV.5b 
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The  lift  force  W  for  other  bearing  plan  views  will  also  be  expressed  as  EIV.f>b  but 
of  course  with  Kw  taking  on  a  form  reflecting  the  actual  geometry  in  point. 

Most  practical  land  and  pocket  geometries  do  not  lend  themselves  to  closed  form 
solutions  as  was  done  for  the  case  of  Fig.  IV.2  in  the  illustrative  derivation.  Confor¬ 
mal  mapping*1)  could  be  used  to  analyze  a  few  more  plan  view  patterns.  Electrical 
conducting  sheet  analog  tests  is  an  experimental  alternate  to  conformal  mapping. 
Reference^  contains  a  directory  of  various  bearing  land  plan  views  for  which  p(x,y) 
were  determined  from  electrical  resistance  sheet  analogs.  was  subsequently  deter¬ 
mined  by  numerical  integration.  Specific  load  coefficients  Kw  could  be  abstracted 
from  the  appropriate  performance  figures  of  reference*2). 

The  analog  method  used  for  reference*2)  as  well  conformal  mapping  are  limited  to 
the  solution  of  the  Laplace  equation.  These  methods  have  been  largely  superceded  by 
numerical  approximation  methods.  The  Finite  Difference,  Finite  Element  and  the 
Boundary  Element  Methods  not  only  readily  solve  the  more  general  Reynolds  equation 
but  also  compute  various  related  subsidiary  data.  A  simplification  of  the  Finite  Ele¬ 
ment  Method  well  suited  to  the  lubrication  equation  is  described  in  Appendix  A1V.3. 
At  this  point  however  it  is  expedient  to  continue  with  the  special  case  of  the  circular 
hydrostatic  bearing  to  explore  some  of  the  general  behavior  characteristics  of  EX.P 
bearings. 
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Flow  Rate  Through  the  Bearing 


The  volumetric  flow  rate  Q4_5  through  the  clearance  domain  4-5  and  the 
corresponding  pressure  drop  due  to  laminar  flow  appears  from  Appendix  A IV. 2  as  a 
fluidic  equivalent  of  Ohm’s  law  namely 

•/*L.(b/a)  ^ 

PV5  “  ^3  * 

which  takes  on  the  more  familiar  Ohm’s  law  form 


P4/5  ~  R4-5Q4  5  101 V .5c 

where  R4_5  =  fluidic  resistance  to  flow  from  contour  4  to  contour  5.  The  resistance 
can  be  written  as 


R  —  •  s 

n4-5  —  3  °r 


EIV.5d 


Sr  is  a  function  of  the  shape  of  domain  4-5.  In  the  case  of  Fig.  IV. 2  Sr  is  easily  calcu- 
Ln(b/a) 

lated  to  be  - .  For  the  general  case  it  would  be  computed  as  a  part  of  the 

7T 

numerical  analysis.  Reference  (2)  plots  conductances  for  other  land  geometries  from 
which  specific  Sr  could  be  evaluated. 


The  important  observation  to  be  made  here  is  that  in  general  the  clearance  resis¬ 
tances  are  proportional  to  h  3  and  viscosity  n  and  that  Sr  is  function  of  geometric  pro¬ 
portions  and  not  the  actual  size  of  the  clearance  domain. 

The  cubic  form  of  the  clearance  is  a  significant  factor  in  the  attractive  behavior  of 
bearings.  The  h3  term  also  explains  the  earlier  assumptions  that  the  pressure  drop 
across  the  pocket  3-4  is  negligible.  Consider  a  normal  situation  where  h  in  region  4-5 


is  about  50  /tm  and  the  depth  of  the  pocket  is  10  x  larger  (  c  =  .5  mm)  then  the  .  esis- 
tance  R4-5  is  103  x  larger  than  113-4.  Hence  the  pressure  drop  within  the  pocket 
domain  3-4  is  insignificant  by  comparison  with  the  pressure  drop  in  domain  4-5. 


Static  Stiffness  k„  of  Linearly  Sliding  Hydrostatic  Bearings 

Let  ks  be  the  stiffness  resisting  changes  in  the  nominal  film  thickness  h  due  to 
changes  in  nominal  bearing  load  W.  Then 

k-  -  -it  or  k’  -  -  -w£^  E[Y-« 

where  p',/0  = 

3/0  dh 

The  nature  of  the  pressure  p3/0  is  defined  by  the  details  of  the  circulating  system. 
The  most  common  circulating  system  for  EX.P.  bearings  is  based  on  a  constant  pres¬ 
sure  source.  The  pressure  regulator  I\R.  maintains  P[/0  at  a  constant  value  while  the 
flow  rate  is  left  to  vary  with  variations  in  h.  A  less  common  system  is  based  on  a  con¬ 
stant  flow  source  without  restrictor  2-3.  In  this  case  p,/0  is  left  to  vary  with  variations 
in  h.  The  constant  flow  source  will  be  discussed  later.  With  a  constant  pressure 
source  an  intermediate  restrictor  2-3  is  necessary  if  the  stiffness  ks  is  not  to  vanish. 

(lauge  pressure  p3/0  is  related  to  the  pressures  at  the  other  terminals  by  the  loop 
equation  for  an  across  variable  defined  in  Section  I.  Applied  to  the  specifics  of  Fig. 
IV.  1,  the  loop  equation  is  p3/0  p3/2  +  p2/,  -4  p1/0.  Since  p2/1  ~  0  the  preceding 

equation  becomes  p3/0  =  p3/2  1  p,/0. 
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.  .  P2/3  P2/3  P2/3 

Let:  O  =  -  or  - or  - 

P3/0  P4/5  P4/0 


It  follows  that  p4/5  =  P|/0  - From  Appendix  A1V.2  lift  W  of  EIV.6  can  now  be 

'  '  1-K* 


written,  in  terms  of  the  fixed  supply  pressure  p1/0  as  W  =  A  p,/0(  — - )KW  where  </»  = 

'  </>+! 


The  flow  restrictor  2-3  acts  as  a  resistor  R2/3  to  create  the  pressure  drop  p2/3- 
There  are  three  types  of  passive  restrictors  in  common  use  for  hydrostatic  bearings  i.e. 

a.  Laminar  Flow  Restrictor  with  a  fixed  resistance  R2_3.  The  most 
usual  form  is  that  of  a  capillary  tube  of  length  Lc  and  an  ID  of  dc. 


The  constitutive  equation  in  that  case  is  Q2  3  =  - ——  .  p2^3 

128// Lc 


b.  Orifice  Restrictor  in  the  form  of  Fig.  IV. 3a  below.  The  orifice 


discharge  freely  into  the  pocket  domain  3-4.  The  flow  pressure  rela¬ 
tionship  is  the  basic  orifice  equation  of  incompressible  flow. 


Q23  =  cg- 


-(P2/3)  where 


Cp  =  coefficient  of  discharge 


tod 


dQ  =  orifice  diameter 
l>  =  mass  density  of  fluid 

Inherently  Compensated  Restrictor  or  the  name  preferred  hen  I’rox 


e 


imity  Restrictor.  This  is  a  variation  on  ihe  simple  Orifi> .  Rest  n.  tor 
as  in  the  Fig.  IV. 3l>  below.  The  orifice  is  plated  so  dose  to  1I1. 


ing  bearing  surface  that  the  restricting  area  to  How  e.  not  1 1.,  or  it, 


A0-A181  816  THE  SCIENCE  OE  AND  ADVANCED  TECHNOCOGV  EOS  2, 
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UNCLASSIFIED  DEC  8b  N00014-83-IL-0385  F/G  13/8  NL 
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-do  Trd* 

area  -  but  rather  the  area  7rd»h  as  shown  because  7rdah  <<  - 

4  1  4 


ie  —  >$>  d2h. 
4 


Fig.  IV.3a  Orifice  Restrictor 


Fig.  IV.3b  Inherently  Compensated  (Proximity)  Restrictor 
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The  constitutive  equation  now  is 

1 

The  various  restrictors  might  be  said  to  be  represented  by  the  generalized  con¬ 
stitutive  equation 

Q2.3  =  ^2-3  h  P2/3  EIV.7 

where  the  exponent  a  assumes  values  of  1,  — ,  —  respectively  for  the  restrictors 

2  2 

defined  above. 

From  the  series  arrangement  of  resistances  of  Fig.  IV.2, 

Q<  5  =  Q2.3  hence  K4^5  .  p4/5  =  C2_3  .  h'1  .  p2/3 
where  K4_5  =  (R4_s)_1  *  *  *  Conductance 


Differentiating  the  above  two  sides  w.  r.  to  h 


dh  ^4~5  ’  P4/s)  —  ^2-3  ^  ^  h  ^Pz/3^  +  ^p2/a)  *  •  (P2/3) ) 


or  Q4  5 

[»A  +  ^1 

1! 

<0 

ro 

U 

«  pEJIL 

p4/5 

1  1 

h  P2/3 

Since  p2/3  =  P2/1  +  P1/0  +  P0/5  +  P5/4  +  P4/3  or 

P2/3  =  0  +  P1/0  +  0  +  P5/4  +  0 

P2/3  =  0  +  0  +  0  —P'i/s 

which  when  substituted  into  the  r.h.s  of  the  above  equilibrium  equation  yields. 


hence 


r  ■ 
1 1 


(  ■  "i 

+A  i 

#'«  i 

•V ' 


7^(1  +  4 )  +  f  (3  -  a)  =  0 

P4/5  0  h 


Applying  EIV.6 


,  ,,,P  4/5  W  ,  3  —  Or  , 

k-  “  “  "T  'TTw1  °r 


EIV.8 


Optimizing  k„  w.r.  to  choice  of  <j>  — 

P4/5 


This  optimum  occurs  at  — —  =  0.  Hence 


dk.  —A 

d?  "  dj  (W ’  <*“"> 


— 2 — )  =  _(3_^  iW  .  ,_jL.  .  _W_ 

0  +  /?J  1  j  d0  (0+/?)+0+/? 


W</) 

(*«*)* 


where:  ^  K  •  p</‘)  -  ^  (Ak'p'/»  (TTW 


•  dw  . .  1 

••  d0  Ak.p1/0(  (]+^)J) 


- -  =  -W — - — ) 

-0)2)  (1+0J 


Substituting  dW/d0  into  the  previous  equation  for 


dks  _  —(3  —  o?)w  —0  0 

d0  0  +  fl  1+0  0  +  /? 

dk,  1  1 

for  70  =o  =  1“<A(TT0  +7T77)  =  0  Hence</,2  =  /y  0oPt  =  (jl/* 

The  final  tabulation  of  stiffnesses  associated  with  the  three  types  of  upstream  res- 


•Vv,  '■ 


>  -  :  ,-Mv 
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trictors  is  given  by  TABLE  II. 

TABLE  II 

Performance  Characteristics  by  Restrictor,  R2_3  Types 


Restrictor  2-3 

C2_a 

a 

P 

k« 

mm 

a.  Capillary 

7rd< 

0 

1 

-W  4> 

1 

1  5  w 

128/xLc 

h  °  <j>+l 

■L5  h 

b.  Orifice 

c^vT 

0 

1/2 

1 

P2 

w 

-1.757  — 

h 

c.  Proximity  * 

c„-VF 

1 

1/2 

1 

P2 

W 

-1.171  — 
h 

</,= 


P2/3 


W  =  A. 


Pi/o- 


1 


P4/5  1  +</> 

*  Inherently  compensated  restrictor 


_  _  w_ .  <!> 

h  ^  '^'{t+P) 


Q  —  C2_3h  p2/3 


Quantitative  Performance  of  a  Sample  EX.P.  Bearing. 


Given  the  nondimensional  load  coefficient  Kw  and  the  flow  resistance  shape 
factor  Sr,  the  quantitative  performance  figures  are  readily  generated  as  demon¬ 
strated  in  the  Appendix  ATV.2.  As  mentioned  earlier  Kw  and  Sr  however  must 
come  from  a  numerical  solution  of  EIV.3  for  each  EX.P.  bearing  plan  view.  To 
continue  a  study  of  the  magnitudes  of  performance  the  analysis  of  AIV.2  borrows 
the  Kw  and  Sr  for  a  circular  bearing.  Other  plan  views  will  have  their  own 
characteristic  values  but  most  likely  practical  plan  views  will  have  Kw  and  Sr 
coefficients  differing  by  less  than  60%  from  the  sample  circular  model.  This 
difference  can  not  be  ignored  for  practical  design  but  for  the  time  being  the  circu¬ 


lar  model  is  taken  in  AIV.2. 
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4  .1 

*\ 


The  bearing  chosen  for  the  demonstration  was  that  of  Fig.  IV. 2  with  the  fol¬ 
lowing  specification: 


Diameter  of  contour  5  —  100  mm 


Diameter  of  contour  4  =  50  mm 


Viscosity  /i  =  .07  Ns/m2  (light  oil) 
Clearance  h  =  40  nm 
Restrictor  2-3  a  capillary  with  (j>  =  1 
Supply  pressure,  Pi/0  —  3.45  x  105  pascals 
Sliding  speed  =  0.2  m/s 


From  these  given  data  the  following  performance  figures  were  determined: 


1)  Load  carrying  capacity  W  —  731N 

2)  Oil  flow  rate  Q  =  1.19  x  10~7  m 3 /h 

3)  Pump  power  required.  Power  =  0.164  watts 

4)  Restrictor  2-3  is  a  capillary  tube  0.75  mm  ID  and  length  of  164  mm 

5)  Static  stiffness  =  27.4  x  10®  N/m 

6)  Coefficient  of  friction  =  .0028 

7)  Reynolds  Number  at  the  inner  edge  of  the  4-5  field  was  3.8 

The  data  above  such  as  the  coefficient  of  friction,  low  power  requirements  etc. 
speak  adequately  for  themselves.  But  one  must  compare  the  bearing  stiffness  to  the 
flexural  rigidity  of  a  cast  iron  beam  to  adequately  appreciate  the  stiffness  of  a  hydros¬ 
tatic  bearing.  The  stiffness  ks  reported  (27.4  x  10®  N/m)  is  identical  to  that  of  a  cast 
iron  box  beam  100  mm  x  100  mm  x  1000  mm  with  s  wall  thickness  of  10  mm  and  cen¬ 
trally  loaded  in  bending.  The  major  dimension  of  the  bearing,  its  diameter,  is  no  more 


v,-  ; 

V  • 


* 


7  '  V**v\Y 

i, *)♦***»«*, ^  *'  '** 


r.2i 


than  the  width  of  the  cast  iron  beam  used  for  comparison. 


Lateral  Guidance  and  Bidirectional  Bearing  Performance. 


So  far  the  bearing  was  presented  as  a  device  for  supporting  undirectional  vertical 
loads.  Guideways  on  the  other  hand  must  restrain  in  transverse  as  well  as  vertical 
directions.  Hydrostatic  bearings  can  be  applied  to  all  four  faces  as  in  Fig.  IV. 4  to  do 
just  that.  Along  the  transverse  direction  (x  axis)  the  opposite  bearing  faces  balance 
each  other  hydrostatically.  However  external  disturbances  acting  to  cause  displace¬ 
ments  along  the  x  axis  will  be  resisted  by  a  hydraulic  stiffness.  This  arises  from  an 
extension  of  the  analysis  made  for  Fig.  IV.2  of  the  stiffriess  in  the  z  direction.  Earlier 

W 

the  stiffness  of  a  single  hydrostatic  bearing  film  was  described  as  ks  =  - — .3  (  '  ■  ) 

h  <p+ti 

or  ks  =  3  AKW  (-■■  -).  If  two  faces  oppose  each  other  they  combine  to  double 

the  stiffness  ie  k,  =  2  x  3  (  AKW  (~~~)) 

h0  <p+o 


Source 


Fig.  IV.4  Hydrostatic  Guideway  Two  Dimensional  Translatory  Restraint. 


For  the  case  of  the  undirectionai  bearing  of  Fig.  IV.2  the  supply  pressure  was 
related  to  the  applied  load  W.  For  bidirectional  bearing  arrangements  p,/0  may  be 
arbitrarily  set  to  any  higher  value.  The  higher  the  supply  pressure  the  greater  the 
static  stillness  attained.  For  instance  in  the  quantitative  example  given  earlier  a 
undirectionai  bearing  called  for  pj/0  =  3.45  x  105  pascals  to  support  a  load  of  731N 
and  yielded  a  stiffness  ks  =  27.4  x  106  N/m.  With  the  double  acting  system  of  Fig. 
IV. 4  not  only  is  the  stiffness  automatically  double  in  the  x  direction  but  also  if  the  sup¬ 
ply  pressure  p^o  is  increased  four  fold  the  overall  stiffness  in  the  x  direction  is  now  ks 
=  2  x  4  x  27.4  x  108  =  219.2  x  106  N/m.  If  this  eight  fold  increase  in  stiffness  were 
matched  again  by  a  cast  iron  beam  it  would  be  a  168  inm  x  168  mm  x  1000  mm  box 
beam,  while  the  bearing  geometry  retains  its  original  diameter. 

If  the  transverse  constraint  was  by  boundary  lubricated  surfaces,  sideways  play 
would  be  inevitable  even  if  mitigated  by  gibs  as  in  Fig.  II. 4a.  With  hydrostatic  bear¬ 
ings  on  these  surfaces  there  is  a  sideways  restraint  equivalent  to  an  enormously  stiff 
hydraulic  spring  with  negligible  imposed  drag. 

The  resistance  by  the  upper  and  lower  bearings  to  z  directed  deflections  is  slightly 
stronger  because  of  the  static  load  W  carried  by  this  pair.  As  the  static  load  W  is 
applied  the  upper  clearance  increases  from  ho,  its  zero  load  value,  while  the  lower 
bearing  clearance  decreases  by  the  same  amount.  The  two  stiffnesses  combine  to  pro- 
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k,  =  3  x  p1/0  AKW 


A . 


(-J !—  +  — L- ) 

'  h0-f  Ah  h0— Ah  ' 


or  k.  =  3  x  p.  /„  AK. 


‘S  Pl/o  "“VW  I 


4> 


MM^f) 


where: 


h0  =  Bearing  clearance  no  load  situation. 

Ah  =  Change  in  bearing  clearance  due  to  the  imposition  of  a  load  W. 


Pitching  and  Yawing  Restraint 

The  vertical  and  horizontal  guidance  as  described  above  is  very  effective.  But  a 
bearing  with  a  single  source  and  single  pocket  3-4  in  a  symmetrical  field  3-4-5  will 
develop  almost  no  resistance  to  pitch,  yaw  and  roll.  However  two  bearings  side  by 
side  as  in  Fig.  IV.5  will  use  their  moment  arms  and  kinematics  to  generate  moment 
stiffness.  If  the  bearing  faces  I  and  II  are  exactly  parallel  to  each  other  (fi  —  0  Fig. 
IV. 5a)  both  bearings  will  show  the  clearance  ho.  With  a  pitch  angle  (b)  (i  the  clear¬ 
ance  is  now  wedge  like  and  unequal  in  the  two  bearings.  For  the  left  hand  bearing  the 
clearance  is  squeezed  down  and  therefore  the  flow  resistance  R^_5'  increases  over  the 
neutral  case  of  i'i  =  0.  Because  of  this  increased  resistance  the  flow  rates  Q4»  5»  and 
Q2<3'  decrease.  Pressure  drop  p2y3i  correspondingly  reduce.  That  is  the  pressure 
level  p3y0  will  increase  beyond  the  neutral  level  p3/0  Fig.  lV.5c.  This  means  that  the 
lift  force  W;  of  the  left  hand  bearing  will  grow  as  the  clearance  h'  is  depressed.  On  the 
r.h.s.  the  action  on  surface  3"-4"-5"  is  the  converse  to  that  of  the  l.h.s.  The  two  bear¬ 
ings  together  then  create  a  correcting  moment  to  opposing  rotation  /f.  A  repeat  of  this 


argument  with  another  pair  of  bearings  in  the  transverse  direction  will  provide  yawing 
restraint.  These  same  'iffects  are  realized  by  a  bearing  with  multiple  sources  as  in  Fig. 

rv.i. 

To  compute  the  pitching  stiffness  k  y  one  must  solve  the  Laplace  equation  EIV.2b. 
A  simple  solution  is  no  longer  available  for  any  plan  form  including  the  circular  plan 
exploited  earlier  for  its  simplicity.  Numerical  analysis  is  practically  the  sole  tool  for 
quantitative  evaluation  of  angular  stiffness. 


-  pressure  /3  =  0  - pressure  £*0 

(c) 

Fig.  fV.5  Mechanics  of  Pitching  Restraint. 
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Constant  Flow  Source 

The  constant  pressure  source  described  above  is  the  dominant  choice  for  EX.P.  sys¬ 
tems.  A  constant  flow  source  is  an  alternative  occasionally  preferred  because  this 
yields  a  bearing  of  somewhat  superior  stiffness.  For  a  undirectional  EX.P.  bearing 

W  3  W 

with  a  fixed  flow  rate  the  stiffness  k,  =  3  —  rather  than  k,  =  —  — —  which  is  the 

3  h  3  2  h 

best  that  the  fixed  pressure  source  can  offer.  With  a  constant  flow  source  there  is  no 
control  restrictor  2-3  hence  some  simplification  is  gained  on  this  score.  On  the  other 
hand  a  regulator  capable  of  maintaining  a  fixed  flow  rate,  despite  pressure  variations, 
is  a  more  complex  mechanism  than  a  combination  of  pressure  regulator  and  any  of  the 
2-3  resistors  described  earlier.  Furthermore  with  a  fixed  pressure  source  several  bear¬ 
ings  can  be  fed  from  a  single  point.  On  the  other  hand  if  a  fixed  flow  source  system  is 
chosen  each  bearing  must  have  its  own  constant  flow  regulator.  This  adds  consider¬ 
ably  to  the  complexity  of  the  system. 

Adaptive  Control  of  Hydrostatic  Slideways 

Increasing  expectations  for  precision  guidance  puts  evermore  emphasis  on  support 
stiffness  for  slideways  and  more  error  free  guideways.  A  double  acting  guideway  as  in 
Fig.  IV.4  is  preferred  for  the  sake  of  stiffness  over  the  single  acting  hydrostatic  models 
of  Fig.  rV  1,2.  However  guideway  precision  is  still  limited  by  the  geometric  errors  of  a 
large  complex  structural  element  such  as  a  conventional  guideway.  An  adaptive  con¬ 
trol  whereby  the  slider  is  positioned  to  a  more  precise  reference  line  can  generate 
stiffness  and  precision  performance  beyond  the  inherent  limitations  of  a  structural 
guideway.  No  other  system  of  sliding  or  rolling  support  is  better  suited  to  adaptive 
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control  than  the  hydrostatic  floating  principle  because  the  support  pressure  is  so  easily 
manipulated. 

Reference  6  and  Fig.  IV. 6  from  that  reference  represent  a  basic  adaptive  control 
system  applied  to  the  hydrostatic  double  acting  guideway  principle.  The  passive  con¬ 
trol  restrictors  1-2  of  Figs.  IV  1,2,4  are  now  replaced  with  a  differential  control  valve 
(D.C.V.)  which  is  activated  by  reference  to  a  precision  straight  edge  (SE).  Tuning 
throttles  (TTU  and  TT\  )  replace  the  flow  conserving  and  pressure  generating  resis¬ 
tances  of  the  thin  films  of  Figs.  IV  1,2,4.  The  tuning  throttles  are  essentially  minia¬ 
ture  hydrostatic  pads,  surrogates  for  the  4-5  zone  of  Figs.  IV  1,2,4. 

The  operation  of  the  adaptive  control  system  is  as  follows.  The  tuning  throttles 
TT  have  their  clearances  hT  (upper)  and  hT  (lower)  initially  set  so  that,  at  the 
expected  nominal  load  W  carried  by  the  slider  II,  the  upper  and  lower  valve  opening 
gaps  (Gu,  )  of  the  DCV  Fig.  IV.6b  are  equal.  That  is  the  valve  is  in  neutral  posi¬ 
tion  at  nominal  load  W.  If  operating  conditions  add  a  dynamic  increment  in  load,  <Vw, 
slider  II  will  depress  <>z  relative  to  the  reference  straight  edge  SE.  The  upper  gap  Gu 
of  the  DCV  will  increase  by  <*>r  so  increasing  the  upper  valve  flow  rate.  Simultaneously 
the  lower  gap  Gu  will  decrease  reducing  the  lower  valve  flow  rate.  The  upper  DCV 
flow  feeds  into  the  lower  junction  and  onto  TT^.  The  increased  flow  rate  due  to  Az 
will  cause  an  increase  in  back  pressure  at  J (  hence  the  lifting  force  of  the  lower  pocket 
will  increase  to  counter  the  <W  imposed.  Similarly  the  lower  valve  gap  feeds  flow  into 
J„  and  thence  through  TTU.  The  reduced  flow  from  the  lower  valve  allows  the  back 
pressure  on  TTU  to  decrease.  Hence  bz  causes  a  decrease  in  upper  pocket  pressure 
adding  to  the  lift  force  countering  <Vw.  The  push-pull  pressure  variations  on  slider  II 
causes  it  to  stabilize  near  the  neutral  level.  With  the  pockets  sealed  around  their 
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peripheries  there  is  essentially  zero  flow  rate  from  junctions  J  to  the  pockets.  Hence 
the  pocket  depths  play  no  role  in  positioning  control.  It  then  follows  that  errors  in  the 
geometry  of  these  pockets  have  no  effect  on  positioning  precision.  Fidelity  is  deter¬ 
mined  by  the  precision  straight  edge  SE.  It  is  claimed^  that,  when  this  system  was 
set  for  a  nominal  load  of  1176  N,  a  variation  of  ±  980  N  altered  the  static  position  no 


a 


seat  h, 


a)  Control  System 


b)  D.C.V.  Enlarged 


Fig.  IV. 6  Adaptive  Control  of  Hydrostatic  Guideway.  Ref.  6. 
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APPENDIX:  AIV.I 

Reduction  of  Navier-Stokes  Equations 

The  analytics  of  specific  problems  of  fluid  flow  have  their  origins  in  the  general 
equations  of  Navier-Stokes3,1,5.  This  set  of  equations  is  based  on  the  equilibrium  of 
forces  acting  on  the  faces  of  an  infinitesimally  sma.l  fluidic  cube.  The  N-S  implicitly 
include  all  of  the  stress  terms  involved  in  the  equilibrium.  This  total  rigour  results  in 
equations  that  defy  closed  form  solutions  for  the  most  general  case.  Specific  problems 
however  allow  order-of-magnitude  comparison  of  terms  in  the  N-S  equations  and  a 
subsequent  reduction  of  these  equations  to  more  tractable  forms.  Such  is  the  case 
with  the  lubrication  problem  where  reductions  described  here  lead  to  the  Reynolds 
p.d.e.,  EIV.l. 

The  foundation  of  the  N-S  structure  is  the  equilibrium  of  intensive  forces  JF^N/in3) 
acting  on  cubic  element  dx,  dy,  dz.  Using  the  notation  of  Section  1  this  equation  is 

r,'  =  rEe  -h  :? 

where:  W » e  =  inertial  intensive  force  acting  on  the  cube  e. 


E.e  =-  Intensive  force  that  earth  exerts  on  cube  e.  That  is  the  gravita¬ 
tional  body  force  acting  on  e. 

s  e  -=  Intensive  force  that  medium  S  contacting  the  cube  e  exerts  on  it. 

Force  3*  Se  is  compounded  of  contributions  from  the  normal  stress  (pressure  p)  and 
parallel  stress  (shear  r).  The  latter  stress  arises  from  tin*  viscous  drag  on  the  faces  of 
e  hence 


»  SC  =  »  no  T  3  . 
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■  M 


I,1' 


1 


1 


The  complete  intensive  force  equation  becomes  E.AIV.la 


^  <•  '  =  ^  E.e  +  &  p.e  +  I'.e 


E.AfV.la 


The  necessary  Continuity  Equation  for  the  cube  e  is 


E.AIV.lb 


(>  =  mass  density  of  fluid 


u  =  velocity  vector  of  fluid  at  the  cube  e 


=  u  i  +  v  j  +  w  k 


<\a  =  mass  leakage  rate  from  e  to  a  sink  s  (per  m3  of  element  e). 


The  Through  terms  of  E.AIV.la  and  the  Across  terms  of  E.AIV.lb  are  related  by 


the  appropriate  Constitutive  Laws  E.AIV.lc,  c\  c". 


7  ■"■■■  Du  ...  r7\f  r 

^  +  (U-V)U 


E.AIV.lc 


where  U.V  =  u ■—  +  v~ — h  w— — 
ctx  cry  ch 


*  lie  =  V gfE/e 


E.AIV.lc’ 


&  p.e  =  VP  =  4 

F  dx  '7  <'Z 


E.ArV.lc" 


The  fourth  term  e  has  an  associated  constitutive  law  which  does  not  lend  itself 


to  abbreviated  expressions  as  those  above.  It  is  now  expedient  to  represent  W  e  as 


q  j  c 

the  matrix  E.AIV.ld  the  elements  of  which  are  taken  from  references  ’  ■  . 


OJ,  -~ 

^  lie 


(at  +  ii2  +  a3) 
(bj  +  b2  +  b3) 
(c,  +  c2  +  c3) 


E.AIV.ld 


_ _ _ -L  *  •  i'"  'l  'yA  *  %*  1*  •  »  'l  -  J#  '•  *1-  ' ■  v  ■'  *  *  '  *  A  /'  /*  "  «• 
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() 


&l  (hi 


K  _JL 

1  "  0y 


ci  — 


_d_ 

ch 


/i(2—  -  2/3  V.U) 
(ht 


n{2~—  -  2/3  V.U) 

()y 


n(2~~  -  2/3 V.U) 
dz 


E.AIV.ld’ 


v.u  =  +  4r  +  -v- 

i  /x  (>\  <>i 


i) 

*2  “  dy 


b2  —  ~~~ 

()z 


c) 


,  <‘>U  <7v 

"(7^  +  a 


.  (>v  ,  <iw  . 

,,(7^  +  7^T» 


.  (>w  ,  <7u  . 


a3  — 


dz 


)  l3 


b’  ~  fi 


,  flw  ,  <7u  . 

,,(7ir +  ar> 


,  (h  ,  <  >v  . 

,,(7v  r  TaT 


;  c3  = 


_c/_ 

dy 


,  civ  c/w  . 

+  V 


<7w 


E.AIV.ld" 


Simplification  by  Orders  of  Magnitude  Deletions 

Equation  E.AIV.ld  is  a  major  hurdle  to  meaningful  analysis.  The  operating  dimen¬ 
sions  and  materials  for  lubricating  systems  allow  the  following  simplifications. 


1.  The  film  thickness  h  (z  axis)  at  any  point  of  bearing  is  about  three  orders  of 
magnitude  smaller  than  the  transverse  flow  path  length  that  the  fluid  will 
traverse.  It  follows  then  that  the  velocity  gradients  --elate  as 


'Mi 

<h 


»> 


y-  »> 


dvr 

dz 


«<  — 


rMi 

(Mi 

/Me 

r7y 

dv 

ih 

(h 

dy 

<h\ 

(hr 

dz 

'  dz 

-  101  - 


AIV.l  4 


«»*  ' 

<S1 


•  I 


;v‘l 


Wl 


\  t'Ji 

\>v 


Hence:  V.u 


/<u  i/v 
t>z  '  i h 


a9  a-i 


b2  »>  b3 


Furthermore  any  mixed  second  derivative  is  orders  of  magnitude  smaller  than 
<Pu  <>2v 

ih2  °F  !h2 

Deleting  the  higher  order  of  magnitude  terms  reduces  E.AlV.ld  to  E.AIV.le 


<)  .  f^U  , 


W  — 

^  u.e 


0  ,  'Ov  . 

75-(,‘75T) 


E.AIV.le 


2.  ^~E.e  «<  ^  ,,.e  because  the  weight  of  a  thin  (li)  element  of  lubricant  is  negligi¬ 
bly  small  compared  to  the  viscous  terms. 

3.  Referring  back  to  the  flow  in  the  thin  film  is  essentially  laminar  so  that 

kinetic  effects  are  small.  Except  for  unusual  systems  the  inertial  stresses  are 
negligible  compared  to  shear  stresses  arising  from  viscous  drag  i.e.  &  e  *  «< 


4.  Finally  <<<  /Ip jih.  hence  also  - — ~{ll~ — )  terms. 

( 'z  <  h 


Hence  the  N-S  equations  for  purposes  of  more  lubrication  situations  is 


.  ..  -  •  •  •.%»  yV.'.fo v  v'.v  v  V  .h-y;  xytr/'/v-1.  ■, 
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(hp/ih. 

dpjdy 

0 


+ 


(f  .  <>u  V 

- (ll - ) 

ih  v  <h  ’ 

()  .  <hr  . 


E.AIV.lf 


Derivation  of  the  Reynolds  Equations  for  Thin  Film  Lubrication 

The  generalized  Reynolds  Equation  E.IV.l  is  based  on  continuity  of  flows  in  and 
out  of  the  sample  volume,  Fig.  AIV.l.  Surfaces  1  and  11  represent  the  two  bodies  of  a 
bearing  separated  at  points  X|/0,  y^0  by  the  distance  h  =  inn.  Let  M.  „  f  be  the  mass 
flow  rate  from  face  a  into  the  control  volume  f  using  the  subscript  notation  rules  of 
Section  I.  The  conservation  law  described  in  Section  1  becomes  in  this  specific  case 

w  f  +  Al  e.f  +  sf  +  *$  n  f  +  H.f  +  M  if  +  M  t  f  +  V||  f  +  hi  =  0  E.AIV.lg 


but  Jl  w  f  +  Jl  c  f  =  M  w  f  -  Jl  f  e  =  w-e  d* 

(be 


and  At  s{  +  A  nf  =  -~-M  s_ndy 

<>y 


A  nf  =  _/)WlI  dx  dy  +  />un  dy(-^)n  dx  +  /’v„  dx(-^-)„  dy 

(be  oy 


where:  un  =  xI1/0  ,  v„  =  yII/0  ,  w„  =  z,I/0 

pwjjdxdy  =  Mass  flow  out  of  f  due  to  the  velocity  of  element  dx  dy  at  II  moving 
away  from  f. 

pV||dy(— — )n  dx  =  Mass  flow  rate  into  volume  f  due  to  the  element  dx  dy  at  II  travel- 
ox 

e  ( )z 

ing  into  f  at  a  velocity  Uj|( - )«.  The  latter  being  a  component  of  Un 


in 


AJV.l  6 


along  the  outer  normal  to  the  surface  at  point  II. 


/»v„dx( 


ih 

(}y 


)n  dx 


Same  general  mechanics  as  above. 


A  If  =  — /'Wjdx  dy  —  /»u1(~)i  dx  dy  — />v,(~-)1  dx  dy 

<hi  ity 

ui  —  *1/0  vi  =  yi/o  wi  =  *1/0 


—  Mass  flux  rate  through  the  skin  (pores)  at  point  I  into  the  control  volume  f. 


,V  f  =  Mass  flux  rate  through  (he  skin  at  point  II  into  the  control  volume  f. 


A  .  f  —  -h  -Tp-  dx  dy 
(ft 


Entering  these  definitions  into  E.AIV.lg  produces  continuity  equation  E.AIV.lh 

Hit  ,  ,  i)A  . 


(hi 


<>y 


ih 

ih 

p 

"wII/(  +  uu 

ih. 

ll 

ih 

...  +  vn 

ih 

i)z 

i)y 

-  VI 

II 

;}y  i 

+ ... 


u 


Wx  dy  +  ... 


...  +  i^if  +  "i r — 

(  'l 


E.AIV.lh 


i  )A 
(hi 


dx  —  dx  dy 


'  )  r  t  . 
-  z'u.dz 

i>x  Jo 


E.AIV.li 


ih\  i)  i  hi 

From  K.AIV.lf  —  - —  /'  ~ — 
i  (x  <h  <h 


Integrating  twice  with  b.c  at  z  0,  u  u(l  at  z  -  h,  u  —  U(j  yields 


u  —  — 


2  ,2 


(h2  -  z2)  +  (z/h)  un/j  +  u, 
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Integrating  and  differentiating  as  per  E.AIV.i  yields 

dx  -  dx  dy  +  "") 


In  a  similar  fashion 


dy  =  dy  dx  +  vn) 

3  12  /i  dy  2<Jy 


Substituting  these  definitions  into  E.  AlV.lh  and  recognizing  that 

iXph) 

*\l/\/>  +  h<>p/«/t  =  ;)i 


it K  it  +  it K  %  “ "  "O  “ t ),! + ,A% )u I + 


...  +  -^-(/>l.U)  +  -£•  (/.hV)  -H 
/lx  4y 


^  -  Win  +  hu)  R.rv.i 

<  a 


where: 


K-^ 

12// 


u,  +  u„ 


v,  +  v„ 


; 

ft- 


'  •  ,  ■  .  ..  I  ■  m  •  -  -  ’  |*  E  ‘  j  -  ’  i  *  >  •*  ,5*  *  V  ■  J  *  i'  4«w’  i  *  <  *  •  .  *  3  v.  i  **  *  (  '  .  '_»•  *j.  *  ».*,1 

,  ■  ,*  .  ,  '(•*,  ,  ;<v  4“;  }«  i*  t*  i '  V  •'  ‘J.+  >  «  **<|,*“  ■  -  •  >  *  •  »  *  ?,*  *  *  .  -  *  ‘ 


V 

' 


i 


wn=  z*/0 


b*  n 


aK  V 

V' 

d  i 


w- — 1-  + 


Wl=Zl/0 


AIV.l  8 


-v 


Body  n 


un=xn^ 


,vi=y 


u«"x*b 


Body  I 


n,e,s,w  (Faces  of  Control  Volume) 


zvh 


Fig.  AIV.l  Element  of  Thick  Filin  Between  Moving  Surfaces  I  &  II. 


H 
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APPENDIX:  AIV.2 


Pressure  Distribution  in  Clearance  Zone  4-5  Fig.  IV. 5 


Differential  equation  for  zone  4-5  is 


i£. 


2  d2P  ,  dp 

F  d?  +  r  dr  =° 
dS 


EIV.4b 


Let  S  =  then  equation  EIV.4b  becomes  r2  +  rS  =  0.  Which  in  turn  yields 
dr  dr 

Ci  dp 

S  =  -  =  and  finally  p  =  -C.L.(r)  +  C».  Assuming  b.c.  at  r  =  a,  p  =  p4/0 

r  dr 

and  at  r  =  b,  p5/0  =  0  the  coefficients  C  become: 

Ci  =  T  (1/  \  >  C2  =  C1Ln(b) 

Ln(b/a; 

Applying  these  coefficients  the  solution  for  the  gradient  and  the  pressure  of  the 

dr 

circular  pad  bearing  are: 


dp  _  _  P4/0  1 

dr  r  Ln(b/a) 


Lifting  forces  Wf„W/-  and  W 


p  _  L„(b/r) 
P4/0  Mb/a) 


EIV.5a 


Wp:  lift  contributed  by  pocket  bordered  by  contours  3  &  4. 


Wp  =  p3/o  •  *ra  assuming  p3/4  =  0 
W*>:  lift  contributed  by  clearance  domain  4-5. 
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Let  (-J-)2  =  q  and  (-p)2  =  qa  then 
b  b 

W(  =  u  \  C  dq  hence 

-Ln(qa)  q* 

=  (1  _  L  (qJ  _  ) 

Ln(qa) 

Let  W  =  Wp  +  W,  (total  lift) 

W  =  p3/07ra2  +  -£^°  (1  -  qaLn(qa)  -  qa) 

Recall  that:  p3/0  =  p 4/0  and  a2  =  b2q  equation  above  becomes 


W 


qa 


=  P3/0  ^  qa  +  P3/on  (7-7—7  -  qa  ~  7 - )  or 

Ln(qa)  Lnqa 


W  —  P3/0  /T^2 


H^-La/b-fl 

Ln(b/a)2 


EIV.5b 


Therefore  W  =  p3/0  .  A  .  Kw  ;  Kw  is  the  lift  coefficient  for  the  wetted  field  3-4-5  of  Fig. 

fV.3 


Kw 


1  ~  (a/ b)2 
L«(-)! 

a 


Other  plan  views  with  there  own  shaped  contours  3-4-5  will  have  their  own  dimension¬ 
less  expressions  for  Kw.  These  coefficients  in  most  part  must  be  computed  by  numeri¬ 
cal  analysis. 


The  flow  entering  border  contour  4  and  exiting  the  border  contour  5  is  Q4  5  using 
the  subscript  notation  and  polarity  rules  for  through  variables  of  Section  I.  Flow  rate 
Q4  5  is  related  to  p4/s  with  a  constitutive  law  appropriate  for  the  domain  4-5.  The 
constitutive  law  can  be  borrowed  from  electrical  mechanics  by  way  of  the  analog  relat¬ 
ing  one  dimensional  laminar  flow  Q  to  electrical  current  flow  I.  The  current  flow  along 
a  conductor  t  units  thick,  B  units  wide  and  a  material  conductivity  k  is  described  by 

cIg 

the  constitute  equation  I  =  — ktB— — . 

dx 

Consider  the  sample  element  of  the  film  in  the  clearance  field  between  r  and  r  +  dr 
of  Fig.  IV.3.  The  fluidic  constitutive  equation  can  be  generated  by  the  following  ana¬ 
log  transformations: 

1  —*•  Q,  E  — *•  p,  dx  — ►  dr,  B  — ►  27rr,  k  — ►  cr 


The  analogous  constitutive  equation  becomes 


dp. 


H  =  viscosity 


Qr.(r+dr)  =  “Cf  h27rr  where  cr  =  — 

A  note  on  cj:  Traditional  fluid  mechanics  does  not  explicitly  designate  a  property 
cf  named  here  "fluidic  conductivity".  This  term  is  introduced  here,  however,  because  it 
provides  a  necessary  link  in  analysis  by  analog.  It  is  the  property  of  the  fluid  lamina 
bounded  by  the  parallel  walls  of  the  clearance  field. 


The  final  form  for  the  constitutive  equation  of  the  circular  field  4-5  becomes 
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dp  C]  _  P4/0 

Earlier  the  gradient  of  p  was  shown  to  be  — =  -  with  C,  =  - L— 

dr  r  M-) 

a 

ermore  for  this  illustration  p4/5  may  be  substituted  for  p4/0  since  p5/0  =  0. 


Furth- 


Hence: 


Qr.(r  +  dr)  =  Q4  s  =  or  in  the  format  of  a  fluidic  "Ohm’s  Law".  That 

) 

a 

is  Q4  s  =  7T-  —  p4/5  where  resistance  =  R4_5  =  Ln(b/a).  For  other  shaped  clear- 
R  -  '  ,r'' 


4-5 


7ThJ 


6ii 

ance  fields  with  uniform  clearance  h,  R4_5  =  Sr  where  Sr  is  a  resistance  shape  fac¬ 


tor  usually  determined  by  numerical  computation. 


Coefficient  of  Friction  f  — 


Let  D|  u  be  the  apparent  drag  force  vector  that  surface  I  exerts  on  surface  II.  This 
force  is  actually  that  exerted  by  the  fluid  in  the  clearance  on  surface  II.  As  such  it  is 
created  by  the  fluid  shear  and  pressure  stresses  acting  on  surface  II  that  is 


Dj  ||  =  f  rjid  A  +  J p(r,#)  .  .  dA  from  contour  4  to  contour  5. 


civ 

rjj  =  shear  stress  on  surface  II  =  ——(II)  x  /< 

dz 


The  fluid  velocity  pattern  within  the  clearance  h  is  a  composite  of  two  parts  (Fig. 


d  v 

AIV.2-1)  -7- 
dz 


*dv  wdv 
dz  dz 


Fig.  ATV.2-1  Thin  Bearing  Film  Kinematics 


Therefore:  J r(II)dA  =  /iUn/|  J  -j^-dA  +  ^  Jh  .  .  dA 

Hence:  D,  „  =  /xUu/I  /  -j^-dA  +  y  /h— ^-.dA  +  /  .  p  .  dA 

dA  —  r0dr 

The  sliding  velocity  is  taken  to  be  constant  in  magnitude  and  oriented  in  a  fixed 

direction.  The  pressure  gradient  on  the  other  hand  is  assumed  to  vary  with  r  and 

dr 

9.  If  h  and  the  plan  view  of  the  domains  3-4  and  4-5  are  symmetrical  about  the  z  axis, 
the  last  two  integrals  disappear.  Furthermore  if  the  clearance  h  is  constant  as  is  the 
usual  case  for  sliding  tables  the  frictional  drag  reduces  to 


Dui  =  /'U  ii/i 


A3~4  +  A4-5 
h3_4  h4_5 


E.AIV.2a 


Equation  E.AJV.2a  prevails  for  all  symmetrical  plan  views  not  just  the  circular  plan 


illustrated. 
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Quantitative  evaluation  of  hydrostatic  bearing  performance. 

Example: 

Given:  A  circular  pad  bearing  with  the  following  specifications. 
Light  oil  employed  with  a  viscosity:  //.  =  .07  Ns/m2 
Table  sliding  velocity:  Ujj/j  =  .2m/s 
a  —  50  mm,  b  —  100  mm,  (f>  —  1 
Pump  pressure:  pj/0  —  3.45  x  105  pascals  (N/m2) 
h  =  40  /xm  Pocket  depth  =  2.0  min 

Find:  Flow  rate:  Q 

Coefficient  of  Friction:  f 
Size  of  capillary  restrictor 
Pumping  power  required:  Powerp 
Static  stiffness:  ks 
Lifting  load  capability:  W 


Solution: 


Geometric  properties: 


A3_4  =  —  x  502  =  .1963  x  104min2  =  .196  x  10  2  m2 
3  4  4 

A3_5  =  —  x  1002  =  .785  x  104  mm2  —  .785  x  10~2  m2 
4 

A4_5  =  A3_5  -  A3_4  =  (.785  x  10~2  -  .196  x  10"2)  =  .589  x  10'2  m2 


^  1  ••  %  •  *■ . •  '-.I  V. 
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W  = 


‘  !+</» 


•  A3_5  .  kw 


W  =  3.45  x  105  .  -  .  0.54 

1+1 


1  +  (r)2 

Kw  - - j-r  =  0.54 

Ln(b/a)2 


W  =  731  N  (164  lbs) 


Static  Stiffness: 


|  ft 

i 


731  3 /2 

k,  =  - '-L-  k.  =  27.4  x  10®  N/m  (156,000  lbs/in) 

^  40  x  10~®  *  '  '  ’ 


This  is  identical  to  the  flexural  stiffness  of  a  1  m  long  cast 
iron  (E  =  117  x  109  N/m2)  box  beam  centrally  loaded.  The 
dimensions  of  the  beam  being  100  mm  x  100  mm  x  1000  m  with  a 
10  mm  wall  thickness. 


Viscous  Table  Drag: 


tt  A3-4  .  A4-5 

DI.ll  -  I1  UII/I  T  +  T 


r.  „„  „  „  .196  x  10“2  ,  .589  x  10"2 

D,  „  =  .07  x  0.2  - - - - — 

2  x  10-3  40  x  10'® 


E.AIV.2a 


D,  „  =  2.06  N 


Coefficient  of  Friction: 


(0.46  lbs) 


M 


I 


■ 


T..V 


f  =  .0028 


s? 

»*il 


KW1'*  W 
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Flow  Rate: 

P4/5 


Q4  5  = 


R4V 

^4.5  = 


Q4  5 


Pl/o  3.45  x  105 
PV5  ~  j  +  0  "  2 

6 11  Ln(b/a)  ^  6  x  .07 

h3  tt  (40  x  10~6)3 


=  1.73  x  10s  pascals 

7T 

1014  Ns/m5 


Q45  *  1.725  x  10s  x  68  x  10  14 

Q4  5  =  119  x  10-9m3/s 

(.436  in3/ min) 


Pump  Power: 


Powerp  Q,  0  x  p1/0  =  Q4  5  *  P1/0 

Powerp  =  119  x  10~9  x  3.45  x  105  =  .041  Watts 

(5.5  x  10-5  H.P.) 


Assume  that  the  combined  clRciency  of  motor  and  pump  is  about  0.5  x  0.5  =  .25  then: 
Net  Power  dissipation  --  .164  watts 

(.22  x  10  3  H.P.) 


Appropriate  Resistor  2-3. 


The  choice  of  resistor  ratio  4>  ==  1.0  anticipates  a  capillary  type  of  controlling  R2_3. 
This  might  be  a  capillary  tube  with  an  inside  diameter  dc.  Considering  the  probabili¬ 
ties  of  clogging  it  would  be  desirable  to  keep  d,.  to  more  than  0.5  mm  I.D.  From  Table 
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Mb/a) 


clearance  field  4-5  is  C4_5  = -  .  — .  For  the  specific  case  of  Fig.  IV.2  S.  = 

6  H  Sr  7 r 


From  continuity  Q23  =  Q45,  or  C2_3  =  C4_5  if  </>  =  1.0  as  given  hence 


^c4  h3  7T 


128K  c  6 n  Ln(b/a) 


therefore 


€  c  =  x  Ln(b/a)  where  dc  =  .5  x  10' 3  m 
h  l*o 

h  =  40  x  10-6  m 
b/a  =  2.0 

Yielding  €  c  =  .032  m  =  32  mm  (1.25  in) 

If  a  more  prudent  capillary  ID  was  chosen  such  that  dc  =  .75  mm  then  £  c  =  32 

(J-^-)4  —  162  mm  (6.3  in). 

.50 

While  the  capillary  tube  is  the  common  choice  for  laminar  restrictors  the  above 
illustration  does  suggest  some  limits  ou  performance.  It  would  not  be  unreasonable  to 
expect  precision  systems  to  run  with  clearances  h  =  20  // m  in  the  interest  of  increas¬ 
ing  stiffness  ks.  This  halving  of  h  will  call  for  an  eight  fold  increase  in  €  c.  Given  a  dc 
=  0.75  mm  the  corresponding  C  d  will  now  be  8  x  163  mm  or  1.296  m  (51  in)  a  length 
which  could  be  awkward  to  fit  into  assemblies  or  even  to  fabricate.  There  is  however 
an  alternate  not  often  recognized. 

Earlier  in  Section  IV  it  was  observed  that  the  flow  resistance  R4_5  was  independent 
of  the  actual  size  of  the  plan  view  so  that  another  bearing  geometrically  similar  but  of 
a  different  scale  will  have  the  same  resistance  providing  the  fluid  viscosity  and  clear¬ 
ance  h  were  similar.  This  recognition  suggests  that  a  scaled  down  facsimile  (except  for 
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Calculations  above  are  predicated  on  the  flow  in  the  clearance  h  being  laminar. 


That  is  Reynolds  Number  <  2000.  For  the  passage  formed  by  contour  4  of  the  pock¬ 


ets  the  Reynolds  Number  could  be  defined  as 


R*.S±s/h. 


Where:  L4  —  developed  length  of  contour  4 


/'  —  kinematic  viscosity  /i//> 


(oil  with  ii  =  .07,  i'  —  10  7) 


For  the  illustrated  bearing  L4  —  2  n  a  =  2  tv  50  mm,  Q4  5  =  119  x  10  9  and  v  =  10 


(m,s)  or 


#  _  H9  x  10" 


R*  =  ..  x  =  3.8 

4-5  27T .  50  x  10~3  10~7 


Using  a  —  size  facsimile  for  the  restrictor  2-3  increases  the  R/  four  fold  to  a  value  of 


- : — : - : — — - - — 1  - . .  '  '  ‘  •  ‘  ‘  ■  J  *  v  ‘ — — 
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APPENDIX:  AIV.3 


Linear  Triangular  Finite  Element  Method 


The  conductance  and  convection  of  currents  through  and  out  of  a  thin  electrical 
conducting  plate  is  described  by  the  p.d.e.  E.ATV.3a. 


(j  i'  f  J  t  Jp 

—  (kt  — )  +  —  (kt  — )  -  Ce  =  q(x.y) 
( be  i  be  ( fy  i  <y 


E.ArV.3a 


where:  e  =  potential 

k(x,y)  =  conductivity  of  the  plate  material 
t(x,y)  =  plate  thickness  at  point  x,y 

q(x,y)  =  current  flux  from  the  surface  to  an  external  sink 
C(x,y)  =  coefficient  of  convection  from  surface  to  external  sink. 

This  type  of  p.d.e.  prevades  much  of  physics  and  engineering.  The  Reynolds  equa- 

.3 

tion  is  obviously  from  this  family  where  pressure  p  replaced  voltage  e  and  - 

12// 

replaces  k  t.  By  this  substitution  where  h  is  equivalent  to  t  we  can  now  define  the 
fluidic  conductivity  of  a  thin  film  as  being  =  k. 


Reduction  of  E.AIV.3a  to  algebraic  form  by  the  finite  difference  method  is 

equivalent  to  setting  up  a  rectangular  network  of  conductors  as  in  Fig.  AIV.3a  The 

current  flow  from  w  to  c  actually  takes  place  through  the  shaded  band  of  Fig.  AIV.3a. 

In  the  F.D.  model  the  shaded  band  is  represented  by  the  discrete  conductor  Kw  _c  = 

kwtw  +  kct.  Hv  k.t,  +  krtr  dx 

( - )  — In  the  same  way  Ks  c  would  be  shown  to  be  ( - )  — — , 

2  '  dx  5  c  v  2  '  dy 


fii* 

tv 

K* 

>'4 


w 

V 

V 

'  V 
v 

•v  i 
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Fig.  AIV.3a  Finite  Difference  Model 

From  the  conservation  principle  Iw  c  +  Ie  c  +  I,  c  +  In  c  +  I*  c  =  q,.  •  dx  dy.  Apply¬ 
ing  constitutive  laws  (Ohm’s  law  in  this  illustration). 

Kw_c  -d  Ke_c  ee/c  +  K,_c  es/c  d  Kq_c  en^c  d*  Kc_*  g*/c  =  flc '  dx  dy  hence 

Kw-C  +  K,_cee  -  (Kw_c  +  Ke_e  +  K,.,.  -d  Kn_c  +  K._c]ec  +  ... 

-  +  Ks_c  es  +  Kn_c  en  =  qc .  d  x  dy 


In  the  rectangular  network  each  of  the  nodes  plays  the  role  of  the  center  node  of 
five  nodes.  Thus  a  field  with  N  interior  nodes  will  yield  N  simultaneous  equations 
represented  by  E.AJV.3b. 


dx  dy 


E.AIV.3b 


For  the  Reynolds  equation  pc  replaces  ec  and  the  r.h.s.  of  EIV.l  replaces  q^  Equa¬ 
tion  E.AIV.3b  is  solved  for  ec,  that  is  pc  which  integrated  over  the  field  yields  the  lift 


1 
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W  for  the  bearing  field. 


The  Linear  Triangular  Element 


A  network  of  uniform  rectangles,  on  which  the  F.D.M.  is  based,  is  not  the  best 


representation  if  the  rectangular  mesh  does  not  conform  well  to  the  borders  of  the  field 


and/or  if  the  potential  gradients  are  such  as  to  suggest  variable  mesh  spacings.  A 


mosaic  of  assorted  triangular  tiles  (finite  elements)  can  be  made  to  conform  more  suc¬ 


cessfully  to  the  actual  geometry.  As  in  the  net  for  the  rectangular  mesh  it  is  desirable 


to  establish  conductors  for  the  network  lines  between  each  node  of  triangular  elements 


as  indicated  by  Fig.  AIV.3b. 


Fig.  AJV.3b  Linear  Triangular  Finite  Element 


Derivation  of  the  values  of  K. 


The  equivalent  discrete  conductor  for  the  F.D.  model  was  easily  evaluated  from  the 


conduction  through  the  shaded  portion  of  the  rectangular  mesh.  The  triangular  mesh 


if 
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can  not  be  so  directly  analyzed.  An  energy  principle  of  mechanics  however  will  serve 
to  establish  equivalent  conductors.  This  is  the  principle  of  Minimum  Co— content7 
(watts)  which  is  an  extension  of  the  Minimum  Potential  Energy7  principle  (joules).  We 
shall  here  abstract  from  the  principle  of  minimum  Co-content  a  relationship  which  is 
sufficient  for  the  purposes  of  establishing  the  discrete  conductors  equivalent  to  the 
linear  continuum  of  the  triangular  element.  That  is,  for  any  arbitrary  combination  of 
the  potentials  at  nodes  i,  j,  k  the  Co-content,  ,  of  the  equivalent  system  of  discrete 
conductors  must  be  identical  to  the  Co-content  of  the  triangular  element  with  a  linear 
distribution  of  the  potentials  interior  to  it. 

The  determination  of  Co-content,  ^  ,  through  the  triangular  element  with  linear 
potential  distribution,  Fig.  AIV.3c,  can  be  expedited  if  that  figure  is  disassembled  into 
component  figures.  Fig.  AIV.3  ci,  cj,  ck.  The  Co-content  for  each  component  is  a 
simpler  analysis  for  which  we  shall  use  the  figure  ci  for  a  sample  derivation.  Fig. 


A IV. 3d  is  the  plan  view  of  the  ci  model.  The  strips  parallel  to  the  edge  i-j  are  lines  of 
constant  potential. 


’X'K 
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The  contribution  to  the  full  Co-content,  r4  is  d  ^  j  of  the  shaded  strip  of  Fig. 
AJV.3d  that  is 

d^  ,  =  l(s)  .  de  =  kt^s-^—  .  de  (watts) 

ds 


Lj  ;)e  e, 

where  (  ,  =  s  —  ,  -  =  —  — 

S  i) s  S 


sL:  — e:  „ 

Hence  d  €  ;  =  kt - ( - )2  ds 

s  s 


Then  ^  J  (kt)  zdz  ,  z  =  — 

b  O  S 


If  the  variable  kt  is  approximated  as  linearly  varying  from  the  corner  values  (kt);, 
(kt)j  and  (kt)k  the  kt  of  the  shaded  strip  of  Fig.  AIV.3d  is  (kt);  (1-z)  -(-  —  (kt)j  + 


—  (kt)k  hence 


/'(kt)zdz  =  l/2(kt)av  ,  (kt)av  = 


So  that  integration  above  yields: 


(kt),  +  (kt)j  +  (kt)k 


ktL;e 2 


'€  i  = - - 

2S 


SLi  .  1 

but  =  A  (area  of  the  triangle)  =  — 

2  2 


1  Xj  y; 

1  XJ  y; 

1  xk  Yk 


Finally:  %  ,  = 


kt 

4A 


L 2  e2  watts 


By  a  similar  process  the  power  expressions  for  the  cj  and  ck  figures  are 


mmmmi 
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kt  »  2  2 
TT  Lk  ek 


Invoking  the  strategy  mentioned  earlier  that  the  Co-content  of  the  discrete  conductors 
must  be  equal  to  the  Co-content  of  the  triangular  element 


Ki-k  e2  +  Kj_j  e2  =  ^  L?  e,2 


K,_k  ef  +  Kh  «,«  -  <€  ,  =  -£■  Lf  «/ 


Kk_j  ek  +  Kk_j  ek  =  ^  k  =  Lk  ek 


Solving  the  above  simultaneous  equations  for  the  three  K  values  yields. 


[Lj2  +  L l  -  Li2 


8 A  [L*+L*  L>. 


but  L2  =  (xk/i)2  +  (yk/j)2 

Lj2  —  (xk/i)2  +  (yk/L)2 

Lk  =  (xj/i)2  +  (yj/i)2 


L2  +  Lj2  -  Lk2  =  (xk/j)z  -I-  (xk/i)2  -  (Xj/j)2  +  ...  +  similar  y  terms 


but  (Xj/l)  =  xj/k  +  xk/i 


and  (Xj/l)2  =  x2/k  +  2xJ/k  xk/i  +  xkz/L 


Hence  L2  -I-  L2  -  Lk2  -  -2(xj/k  .  xk/i  +  yj/k  .  yk/.) 
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k  i 

Finally  (xj/k  .  xk/i  +  yj/k  .  yk/i)  or 

k  t 

Kj-j  -  +  4A  (xj/k  xi/k  +  yj/k  yi/k)  similarly 

E.AIV.3C’ 

kt 

Kj-k  -  +  4 A  (*j/i  •  xk/i  +  yj/i  •  yk/i) 

E.AIV.3c” 

k  t 

^k-i  —  (Xk/j  •  xi/j  +  y k/j  •  yi/j) 

E.AIV.3c”’ 

Equations  E.AIV.3c  then  are  the  equivalent  conductances  connecting  the  corners  of 
the  linear  triangular  element  of  Fig.  ATV.3b.  It  remains  to  establish  the  branches 
simulating  the  convection  and  external  source  flows. 

Equivalent  Lumped  Convective  Conductors  For  Element  i,  j,  k,  With  a  Fixed  Convec¬ 
tive  Coefficient  C  In  Its  Field 

Let  Q|  op,  Qj  Qkv  be  the  equivalent  lumped  flows  convected  out  of  nodes  i,  j  and 
k  towards  the  ambient  sink  x. 

Let  M  j_j  be  the  moment  of  Qj  v  around  the  edge  j-k  of  Fig.  AIV.3d.  M  ■l_J  is  also  the 
moment  of  convected  flux  integrated  about  the  edge  j-k. 

j-k  =  Qi  x  —  C/e^  ds(l— ~-).S  where  C  =  Li— 

S  S 

Let  s/S  =  u,  then  Qj  x  S  =  S2CLi  /  el/x.(l-u)2  u  +  ^  +  htA  u2(1_u)du 

2 
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Hence:  QI  V  =  (2el/v  +  ej/x.  +  ek/v) 


Qjv  =  —  (ej/v  +  2ej/x.  +  ek/ x) 


Qk  V  —  ^2  (ei/ v  ej/  ^  ^ek/?c) 


A IV. 3  9 


E.ArV.3d’ 


E.AJV.3d" 


E.AIY.3d’ 


Alternately  the  above  equations  can  be  stated  as 


x  '  3  12  12  ' 

q.  v  =  ca  +  M  +  fii) 

x  v  3  12  12  ' 

Qkx=cA(^  +  ^iA  +fiA) 

Xk"x  v  3  I2  12  > 


E.ArV.3e’ 


E.AIV.3e’ 


E.AIY.3e’ 


The  latter  equations  allow  a  pictorial  representation  of  -  Ce  of  E.AJV.3a  by  Fig. 
AIV.3e. 


- 1  CA 


Fig.  A.IV.3e  Network  Representation  of  Convective  Equations  E.AIV.3e. 


YtT 


‘rMsyi,** 


n 
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Equivalent  Lumped  Flow  Sources  for  Element  i,  j,  k  Given  a  Linear  Flux  Field  q. 

Let  Q,  *,  Qj  *,  Qk  •  be  the  equivalent  concentrated  source  flows  at  nodes  i,  j  and  k. 
Let  M  ,_j  be  the  moment  of  Q,  .  about  edge  j-k  as  well  as  the  integrated  moment  of 
flow  q  about  edge  j-k. 

M  k  =  Qj ,  S  =  j  (l-s/S)S  ds.  This  equation  is  the  same  form  as  M  j_k  for 
convected  flows  but  with  a  linear  q  distribution  rather  than  e  distributions  hence 

Qi  *  —  (2qi  +  qj  +  qk)  E.MV.3f 

1  £ 

Qj  ■  =  (Hi  +  2lj  +  Ik)  E.AIV.3r 

Qk..  -  -JJ  (q,  +  qj  +  2qk)  E.AJV.3f”' 

Fig.  AlV.3f  is  total  representation  of  E.AJV.3a  applied  to  the  finite  element  i-j-k. 
The  total  flow  towards  node  i  can  be  expressed  by  reference  to  the  subnetwork  of  Fig. 
ATV.3f. 


Fig.  A.lV.3f  Linear  Element  Network  Representing  E.AIV.3a  On  Element  i,j,k. 
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Conservation  of  Flows  Into  The  Central  Node  of  a  Cluster  of  Elements 


Figure  AIV.3g{a)  shows  a  field  of  linear  triangles  laid  out  to  solve  E.rV.a.  Figure 
AIV.3g(b)  is  a  subassembly  of  n  triangular  elements  such  as  discussed  above.  The  ele¬ 
ments  (in  this  case  n=6)  share  a  common  node  g.  Each  element  of  the  cluster  has  its 
own  properties  (A,  k,  t,  c,  q).  To  identify  a  property  with  a  particular  element,  it  will 
be  subscripted  by  the  addresses  of  the  nodes  making  up  the  corners  of  the  element.  As 
an  example  (CA)jmm+I  combines  C  and  A  of  the  Element  (g,  m,  m+1).  When  the  n 
elements  are  assembled  the  contiguous  edges  of  Fig.  AIV.3g(b)  and  nodes  of  neighbor¬ 
ing  elements  coalesce  their  conductances  and  sources  into  the  star  network  of  Fig. 
AIV.3.g(c)  where 


Q*  •  =  Qe  ‘  (gil»2)  +  Qg  .  (g.2,3)  ...  +  Qg  .  (g,m,m+l)  ...  +  Q, .  (g,n,l)  E.AIV.3g’ 

and  Qg  .  (g,m,m+l)  -  -Jy  (2q<  +  qm  +  qm+1) 

K,.v  =  (^-)(g,l,2)  +  (~~"){B(2,3)  ...(^•)(g,m,m+l)  ...  +  ^-(g.n.l)  E.AIV.3g” 


Elements  EL  (g,  m-1,  m)  and  EL  (g,  m,  m+l)  share  a  common  edge  g-m.  Hence  the 
equivalent  conductors  along  branch  g-m  combine  to  Kg  _m  of  E.AIV.3h. 


Kg_m  +  — 
K  m  l2 


g.m- l,m 


^  (^g-m  "l"  ^2  Vro.nifl 


E.AIV.3h 


From  E. AIV.3  K 


‘‘g-m 


g.m-l.m 


kt 

4A 


.  (xg/m- lxm/m- 1  ^  Ag/m-iym/m-l) 

'g.m-l.m 


A  similar  cycling  of  subscripts  will  generate  (Kg_m)g m m+).  Finally  applying  the 
conservation  principle  to  the  cluster  at  g,  using  E.AlV.3f,g,h,  yields  E.AlV.3i. 


-Q**  -K(g-x)eg/x  +  V)  (Kg_mem/g)  =  0  or 

m-1 

m*g 
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£  Kg_m  em 

m-1 

m*g 


Q 

V 

m«=l 

m*g 


(Kg-m)  +  *V 


K  — 


E.AIV.3i 


Where  em  =  emyx,,  eg  =  eg/>.'  an<^  m  are  the  global  addresses  of  the  modes  sur¬ 
rounding  the  central  node  g.  There  are  N  unknown  values  eg  in  the  entire  field  of  Fig. 
A.IV.3g(a)  hence  there  are  N  equations  of  the  form  E.AIV.3i  to  be  solved  simultane¬ 
ously  either  by  direct  methods  or  by  such  iterative  methods  as  successive  over  relaxa¬ 
tions  (SOR)  of  the  Gauss-Seidel  algorithm.  Direct  methods  such  as  Gauss  Elimination 
are  faster  than  iterative  methods  but  the  latter  are  easier  to  program,  adaptable  to 
non-linear  situations  and  require  less  storage. 

Interest  in  the  F.E.M.  has  generated  a  lot  of  references  for  handling  more  sophisti¬ 
cated  field  problems  than  posed  by  E.AIV.3a.  The  texts  of  Huebnerl8!,  Oden^,  Zienk- 
iewicgl10!  were  referred  to  frequently  to  develop  the  simplified  physical  treatment 
described  here. 


1 


I 

i 

> 

1 


i 
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V  Aerostatic  Bearing  Guidewaya 

An  alternate  to  the  liquid  EX.P  (hydrostatic)  guideway  is  the  gas  EX.P  (aerostatic) 
guideway.  An  obvious  attraction  of  an  air  bearing  is  the  ready  availability  of 
compressed  shop  air  for  the  pressure  source.  But  even  more  important  is  the 
elimination  of  an  oil  contaminated  environment  in  which  the  metrology  associated 
with  precision  systems  would  otherwise  have  to  operate.  The  low  viscosity  of  air 
(/i  air//i  oil  ~  1/3000)  speaks  for  friction  so  low  as  to  consider  such  bearings 
essentially  frictionless.  Not  only  is  any  suggestion  of  stick  slip  eliminated  but  also  the 
frictionless  operation  enhances  the  precision  of  work  piece  positioning  for  numerically 
controlled  machine  tools.  The  mechanics  of  air  bearings  is  closely  related  to  the 
mechanics  of  precision  air  gages.  The  two  principles  joined  provide  a  system  capable 
controlling  work  piece  positioning  repetitively  within  0.02  fim.  However  the 
compressibility  of  gases  in  the  supporting  domains  of  the  bearing  introduces  problems 
of  increas'd  response  time  to  feedback  controls  and  a  susceptibility  to  self  induced 
vibrations.  Considering  these  hazards,  the  design  of  aerostatic  bearings  must  be  tuned 
more  precisely  than  would  be  the  case  for  the  more  forgiving  hydrostatic  bearing. 
That  tigl  ter  tuning  begins  with  a  re-examination  of  EIV.l 

Liquh  bearings  imply  incompressibility  hence  the  density  term  being  constant  is 
cancelled  out  of  EIV.l.  For  gas  bearings  since  the  fluid  is  compressible  the  density 
variation  must  be  reckoned  with.  From  the  thermodynamics  of  perfect  gases  at  any 
instant  or  point  in  the  system  the  pressure  and  density  are  related  as  />(x,  y,  z,  t)  oc 
p1/'  where:  7  is  a  polytropic  constant.  It  is  infinite  for  an  incompressible  fluid, 


approximately  1.4  for  adiabatic  compression  of  air  and  equal  to  1.0  for  isothermal 
compression.  Because  the  thermal  mass  of  the  guideways  is  overwhelming  compared 
to  that  of  the  thermal  mass  of  the  gas  coursing  along  the  guiding  surfaces,  it  may  be 
usually  assumed  that  the  compression  -  expansion  experienced  by  the  film  is  isother- 

h3 

mal,  7  —  1.0.  Equation  EIV.l  therefore  transforms  to  V —  p  V  p  =  12  r.h.s.  (with 

absolute  pressure  p  replacing  p)  or  ’/2Vh3///Vp2  =  12r.h.s..  If  a  flat  EX.P  guideway  is 

considered  with  constant  h  and  r.h.s.  set  to  zero  the  gas  bearing  equation  reduces  to 
2  ^2  2 

—  P  +  — Y  =  0.  If  p2  is  replaced  by  then  the  flat  EX.P  gas  bearing  equation 

ihc  '}y2 

becomes 


cT2^  d2A> 

<h2  rTy2 


E.V.l 


This  forir  is  identical  to  the  liquid  EX.P  equation  E1V.3.  The  techniques  discussed 
earlier  for  evaluating  the  pressure  distribution  are  valid  for  evaluating  vl'.  The  work¬ 
ing  pressure  is  aAF  with  the  restriction  that  the  pressure  constituting  vl/  must  be 
defined  a«  absolute  whereas  the  pressure  previously  defined  as  the  lifting  pressure  was 
loosely  referred  to  as  gage  pressure.  For  given  boundary  pressures  p4/*  and  p5/»  the 
gas  bearing  will  give  slightly  higher  lift  load  W  than  the  liquid  bearing  because  of  the 
expansion  of  the  gas  as  it  approaches  the  border  contour  5.  There  will  however  be  a 
much  larger  mass  flow  with  a  gas  medium  primarily  due  to  the  much  lower  viscosity  of 
gases  when  compared  to  liquids. 


Perhaps  the  most  effective  way  to  compare  hydrostatic  and  aerostatic  performances 
is  by  an  analysis,  as  in  A.IV.2,  of  the  circular  bearing  Fig.  IV. 2  with  air  as  the  medium 


rather  than  oil. 
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The  circular  gas  bearing  equation  is  modified  from  EIV.4b  to 


d2^  .  d'l' 


dr2  +r  dr  ° 


The  solution  follows  exactly  that  given  in  AJV.2  so  that 


d!"  "  f'  Ln(b/a)  ’  =  V  “  V  “  E  V  2» 

Integrating  this  gradient  with  the  boundary  conditions  to  the  4-5  domain  yields 


Comparison  of  Mass  Flow,  A  Gas  &  Oil  Bearing 


E.V.2b 


Referr  ng  to  Fig.  IV. 2  the  mass  flow  from  the  contour  4  towards  contour  5  is 


M.  1  ZL  — 


4.5 - Pa 


h3 


Where: 


Let: 


12/, 

subscript  4  represent  the  state  at  contour  4, 
subscript  5  represent  the  state  at  contour  5. 

subscript  oo  represent  the  ambient  state, 


Isothermal  compressible  flow:  p4  = 


Poo/* 


P4f*  while  incompressible  flow:  p4  =  p^,. 


Let  subscripts  g  and  €  identify  gaseous  and  liquid  properties  respectively  hence 


4.5^  *  (p oo)e  7 —  nr  • 
o/V 


r  Ln(b/a) 


(•«  ,  b)<  = 


(>’«),  ■>,  i 

- s.  — 2—  — >—  .  - - — 

p^y.  12/ij  r  Ln(b/a) 


■M2 
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■>1 


Vl 

'*\;i 

v,  i 


*1**1 

v 


*>> 


■.•.•I 


•« 


$ 

*s 

:S 


Hence 


(•*  4.s)g  =  (Poo)g  _^g 

(•*4.s)*  (*>ocK  l*V 


4.5)1 


=  (^oo)g  ik 

'*  (PoeV  V 


U  [(P4/02  ~  (P5/02) 
.  *V  Pg  2Poo/*(P4/+  ~  P5/* 

3  P£_  (P4/*  +  Ps/*) 

.  I  Mg  J  2Poo/« 


E.V.3a 


E.V.3b 


[i(  —  .07  Ns/m2  as  in  original  oil  bearing  example 
/ig  =  //ajr  =  .02  x  10~3  Ns/m2 


—  =  3500  and 
Pg  (Poo). 


Let  p4/5  =  1.725  x  10“  N/m2  as  in  original  oil  bearing  example  AIV.2. 
let  p 5/»  —  p^/.  =  1.013  x  105  N/m2 


P4/*  +  Ps/*  =  1.725  +  1.013  4-  1.013 

2 p^/.  2  x  1.013 


=  1.85 


(■d  <  5h  1  h,  f  h„ 

=  if-  <*">«»■»-•••»  ^ 


Because  the  kinematic  viscosity  of  the  oil  is  about  5  x  greater  than  that  of  air  at 
ambient  conditions  the  overall  mass  flow  of  air  will  be  8.97  x  that  of  the  oil  example 
for  the  same  film  clearance.  It  is  therefore  the  practice  to  make  gas  bearing  clear¬ 
ances  considerably  smaller  to  conserve  on  the  mass  flow  amongst  other  things. 

In  the  earlier  example  of  the  oil  bearing  it  was  shown  that  the  volumetric  flow  was 
Q4  5  =  .436  in3/min.  It  then  follows  that  SCFM  flow  rate  for  the  air  bearing  would  be 


iSilk.  =  —  e  x  l>( 

{Qi.s)/  &  (  P%  x 


x  -  =  8.97  x  722  =  6476 


Since  (Qvs)f  =  -436  in3/m  from  AIV.2  then  (Q45),,  —  .436  x  6476  —  2823.5  in3/m  or 


r.5 


air  flow  --  1.63  SCFM 


If  the  original  Glia  thickness  of  40  /tin  is  halved  for  the  gas  bearing  the  flow  rate  would 
be  reduced  eight  fold  to  about  0.2  SCFM. 


Friction  Comparison 


Equation  E.AJV.2a  defined  the  viscous  drag  of  the  slide  bearing.  For  the  same 
geometry  and  velocities  the  drag  formulation  is  the  same  for  oil  or  gas  bearings  only 
the  difference  in  viscosity  plays  a  role.  The  coefficient  of  friction  for  the  oil  bearing 
was  .0028.  With  air  as  the  lubricant  the  coefficient  of  friction  would  then  be 
.0028/3500  —  8  x  10^7!  The  bearing  would  be  essentially  without  friction.  Actually  a 
pressurized  film  tapered  only  0.2  arc  seconds  would  create  more  of  a  disturbance. 


Lifting  Crpacity  of  the  Gas  Bearin 


Lift  contributed  by  the  h  region  (4-5)  is 


5 

W4— 5  J4  Pr/oo  •  271?  dr  ,  pr/oo  =  Pr/.  -  Poe/. 


Applying  E.V.2b  pr/.  -  |(pv.)2  -  (ps/.)2|  +  (pt/<)l 


W.-s  "  //  (!(P,/.)’  -  +  (Ps/O*  2"  dr  -//  inr  dr 


Let  r  )5/*  =  pxy* 


Vi  »>  «  *.  V.'-.Vi*  V, 


mm 
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from  the  fluid  state  at  contour  4.  The  actual  choice 
of  designations  is  a  matter  of  convenience  imposed 
by  a  particular  equation. 

The  dynamic  behavior  of  this  system  is  approximated  by  adding  dynamic  incre¬ 
ments  to  the  steady  state  conservation  equation.  For  this  purpose  let: 

m  =  tufl  ,  the  perturbation  from  steady  state  mass  flowrate. 

<;  =  tip,  the  perturbation  from  the  steady  state  pressure, 
w  =  8W,  the  perturbation  from  the  static  lift  of  the  bearing, 
z  =  £h,  the  increment  in  the  static  film  thickness. 

M||_*  =  mass  of  the  oscillating  body  =  Mjj. 

The  total  mass  flow  equation  focused  on  the  pocket  in  the  dynamic  state  becomes 

{M  2  3  +  rh2  3)  +  [M  5  3  +  m5  3)  4-  m»  3  +  rh,,  3  =  0 

or  m2  3  +  m5  3  +  m*  3  +  rii||  3  =  0  E.V.4 


From  Section  IV  the  flow  through  passive  types  of  restrictors  is  in  general  related 
to  pressure  drop  and  orifice  or  capillary  flow  areas.  Hence  perturbations  in  these  flows 
would  be  related  to  perturbations  in  these  parameters  so  that 


—  d.f?  9  9  — . 


()j& 


DM. 


2  3  .  ,  2.3  ,, 

~rt - op3/*  +  — .  oh  or 

(  P3/* 


(Ml 


rn2  3  — 


2  3 


^P3/‘ 


ft/*  + 


(Ml 


2  3 


Oh 


The  steady  state  flow  through  the  domain  5-4  and  into  the  pocket  envelop  is  Jl  5  34 
ex  'k5y3  and  h3.  It  follows  then  that 
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(\A%  c  o  c  o 

m53  =  -  :  .  Spy*  +  — <$h  or 

^3/‘  '  dh 


.  3  *  dJl  5  3 

ms  3  ”  0^7  ’ &  +  IT"  ’ z 


The  l.h.3.  of  EIV.l  contains  the  expression  that  is  h-—^-  +  p- The  first 

dt  d  t  <)t 


term  refers  to  the  compressibility  of  the  fluid  hence  the  capacitance  of  the  film.  Since 


the  film  h  in  domain  4-5  is  extremely  small  this  capacitance  will  be  neglected  initially. 


In  the  domain  of  the  pocket  its  depth  plays  the  role  of  a  very  large  h  hence  the  capaci¬ 


tance  here  must  now  be  considered.  From  Appendix  A.I  Table  A.I.c,  the  capacitance 


of  a  volume  is  C3_»  =  -  ie. 


vol  3p3 


,  '7  =  1  assuming  isothermal  conditions. 


f>  3  x  • 


Hence  m*  3  =  (vol3 - )  p./3  =  -C3_.  p3/.  =  -  C3_*£3 

Pa/* 


The  second  part  />— -  of  — f- —  above  describes  the  contribution  to  mass  flow  by 
at  at 


squeezing  the  surfaces  II  and  I  together.  As  far  as  contributing  to  the  mass  flow  into 


the  pocket  3  the  squeeze  effect  is  rhu  3  where  ih|j  3  — z  />3  A3  with  A3  =  area  of 


pocket  +  a  portion  of  the  area  4-5. 


The  perturbated  mass  flow  equation  E.V.4  becomes 


+  ..  <3  +  77  1  d — 7  Cj/*  H - 77  ^  +  •••• 

<^p3  »h  c?p3/»  rJh 

....— C3_«  <;3/*  —  /'3A3z  =  0 


E.V.Sa 


The  conservation  of  momentum  for  the  bearing  system  is  defined  by 


M|i  z  -  w4_5  +  w3_,,  =  $3/,vA3_4  +  <4/5  SwA4-5 


Is  t 


I 


& 


ft 


Pi 


Si¬ 

ft* 

I 


n 

•V 

*V. 

•V 

> 


& 
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a 

I 


where  A3  _4  is  the  area  of  the  pocket,  A4_s  is  the  area  of  the  domain  4-5  and  Sw  is  the 


J/pdA 


area  shape  lift  function,  Sw  =  Since  f4/5  =  f4/oo  and  C</3  =  0  let  c4/5  = 


6/j:  ‘hen  Mjj'z  =  <,3/oo  [A3_4  +  SWA4_5 


E.V.5b 


Substituting  E.V.5b  into  E.V.5a  yields 


*o  z  +  aj  z  +  a2  z  -(■  a3  i  =  0 


where:  a0  =  C3_«B 


|A3_4  +  SwA4_5 


(X4?  «  o  I  bit  c  i 

- — —  +  -r— B 

^3/*  f^P3/* 


a2  ~  /J3^3 


(bil  53  2  3 


From  Appendix  A.V 


f>M.  n  , 


•  ^  2-3  ^  2-3  = 


&h-ii  i/i+i) 


'^53  3.5 


3-5  ~  ^2* 

1_  P°o/* 


2.3  ^  2.3  , 

^  ^  j  (inherently  compensated)  otherwise  =  0 


Finally 


ao  —  C3_*  B 


4  =  +  ~ 


1  M  2; 

&  2-3  +  &  3-5 

J  Pi  /• 


[a3_4  +  swa4_5] 

1^23  „ 


a2  —  P3A3 


3  5  A  2  3  A  2  3 

a3  —  —  —3  +  - f-  H2_3 

h  n  n 

H2_3  =  2  Inherently  compensated  orifice 
3  otherwise 

From  the  Routh-Hurwitz  test  conditions,  stability  (no  self  induced  hammering)  is 
realized  if  all  aj  coefficients  are  positive  and  a^— a0a3  is  positive,  (a4a2  >  a0a3).  The 
term  ^35  —  ^23  °f  ai  could  conceivably  go  negative  but  the  tabulation  of  these 
terms  in  Appendix  A.V  indicates  that  this  expression  is  always  positive.  The  remain¬ 
ing  test  for  stability  is  that  aja2  >  a0a3  that  is 


(&  3-5  —  &  2-3) 


B/)3AJ>TC3.  .15- 


or  3  5  —  &  2-3)  A3  >  ?  C3  . 

P3/*  “ 


P  3  ^3 

where:  C3_«  =  vol3_4  -  ~  (A  +  h)  .  A3_4  - 

P3/*  P3 /* 


Hence  3  5  -  9  2_3)  A3  >  ? 


A3_4  H2_3 


E.V.6 
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A  =  pocket  depth 


&  3_s  —  &  j_3  =£  2  ,  A3  Total  wetted  lift  area 


A3  4  :  -  pocket  area  ,  Say  A3_4  =  fA3 


E.V.6  becomes  approximately 


2  A3  >  ?  —  +  1  fA3  H2_3 

n 


for  stability  |  —  +  1  j  fH2_3  must  be  less  than  2.  Choosing  an  inherently  compen¬ 
sated  orifice  sets  H2_3  =  2  otherwise  a  more  dangerous  H2_3  =  3  will  result.  The  sta¬ 
bility  condition  reduces  to  —  -|-  1  .  f  <  1.  If  =  10  f  must  be  less  than 

h  h  11 


Using  the  example  of  Fig.  IV. 2  it  would  require  —  < 

b 


VI 


and  if  h  =  20//m,  A 


.2  mm,  a  very  shallow  pocket.  A  more  practical  design  removes  the  pocket  completely 

from  Fig.  IV. 2  and  replaces  it  with  a  channel  along  the  contour  line  4  as  in  Fig.  V.4. 

The  channel  is  fed  by  N  equally  spaced  orifices  along  the  channel.  Since  there  is  no 

pressure  gradient  inside  contour  line  4  the  pressure  distribution  experienced  by  the 

pocketed  design  would  be  duplicated  and  so  would  its  lifting  capacity.  The  ratio  of 

the  volumes  of  the  two  recesses  is  a  measure  of  their  relative  influence  on  the  stability. 

....  ,.  vol.  channel  b2  x  2~a  b22  .  .  .,  ..  .  .  ,  .,  , 

Ratio  —  -  — -  r — -  =  — — .  Assuming  the  dimensions  of  example  cited 


v<  '.  pocket 


7ra  A 


above,  a  100  mm,  A  —  0.2  mm  and  b  assumed  at  2  mm  yields  a  Ratio  —  -  that 

2.5 

is  the  channeled  design  has  2.5  x  the  stability  margin  attained  by  the  .2  mm  deep  full 


pocket. 
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The  analytics  described  above  are  basically  similar  to  the  simplified  lumped  param¬ 
eter  analysis  of  Licht  and  Elrod^  of  circular  bearings  such  as  Fig.  IV.2.  Another 
analysis  by  Licht  and  Elrod^  was  refined  to  include  the  capacitance  of  the  h  film  in 
domain  4-5.  Tests  were  conducted  on  actual  bearings.  Confirmation  by  experiment 
was  somewhat  poor  but  it  did  provide  useful  design  guidelines.  That  is  that  the  simple 
lumped  parameter^1)  analysis  results  are  pessimistic  compared  to  the  refined  analysis^2) 
but  even  the  latter  results  are  pessimistic  compared  to  the  test  data^  when  supply 
pressure  is  less  than  85  psia.  In  the  test  domain  of  supply  pressures  between  85  psia 
and  95  psia  the  simplified  analysis^  gave  somewhat  optimistic  results  on  allowable 
pocket  depths  while  the  refined  analysis  remained  a  bit  pessimistic. 

A  circular  thrust  bearing  is  only  infrequently  used  for  linear  guideways.  More  com¬ 
monly  the  guideways  will  be  rectangular  as  suggested  by  Fig.  IV.  1.  However  for  the 
aerostatic  case  the  pocket  is  to  be  avoided.  Instead  a  series  of  orifices  arrayed  in 
parallel  lines  is  employed.  Stout  and  Sweeney^  have  prepared  design  charts  for  this 
type  of  field.  While  there  are  no  explicit  pockets  the  authors  did  reckon  with  the  una- 

! 

voidable  capacitance  associated  with  each  orifice  subassembly.  The  design  charts  offer 
dimensional  specifications  to  skirt  pneumatic  hammer  instability.  Other  charts  include 
load,  clearance,  supply  pressure  and  orifice  sizes  for  optimum  configurations. 

« 


Stabilizing  the  EX.P  Gas  Bearing  Slideway  With  Additional  R-C  Members 
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5,oo 


HI — 

C3-» 


HI — |i* 

c,-. 


(b) 


Fig.  V.l  Schematic  of  Aerostatic  Bearing  with  An  Added  R-C  Element  Ref.  4. 


Let:  Subscript  12  represent  the  envelop  of  the  chamber  containing  nodes  1 

and  2. 

Subscript  34  represent  the  envelop  of  the  chamber  containing  node  3  and 
the  contour  4. 

Subscript  45  represent  the  chamber  between  contours  4  and  5  occupied 
by  bearing  film  h. 


The  perturbated  mass  flow  equations  of  the  cascaded  system  of  Fig.  V.l  become 


m,  j2  "h  m*  12  ™  *3.12  ”  0 
0*2.34  4-  m*  34  +  m4  34  +  m]]  34 
*3  45  +  **.45  +  *oo45  +  *11.45 


0 
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Let:  q'  =  C3_«  =  VoI(34) 


s'  =  squeeze  flow  phA 


dM  3  5 

()M  4  5 

3 M 

c>h 

<9h 

h 

£  = 

dM  1 1 
dpi/* 

M  g.i 

P2/» 

&  8-1- 

dM  2  3 

M  2.3  . 

v  = 

^P2/‘ 

-  j 

P2/* 

2-3 

dM  2  3 

M  23 

c  = 

dpz/* 

P3/* 

^  2—3 

( 141  t 

* 

=  vol-12 

Pi 

d?i/> 

-  =  Cj_. 

Pi/* 

dM  3 

5  M  35 

2 

’  P3 /‘ 

1  - 

P5/* 

P3/* 

As  in  the  conventional  case  discussed  earlier  the  three  simultaneous  m  equations 
are  converted  with  the  above  influence  coefficients  into  simultaneous  equations  involv¬ 
ing  pertu  “bated  pressures  <5p  and  movements  z  at  the  various  points  of  the  system. 
The  equation  Mn  h  =  My  z  =  J  Ap  dA  is  used  to  express  perturbated  pressures  in 
terms  of  perturbated  motions.  A  fourth  order  differential  equation  is  generated  just  as 
a  third  o’  der  D.E.  was  generated  by  the  earlier  illustration  that  is: 

a0  z  4-  ai  z  +  a3z  +  a3z  +  a^  =  0 


£ 
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Where:  z  =  <$h  perturbated  clearance 

a0  =  qV  ,  Hi  =  (f  +  W  +  (f  +  rj) q' 

»s  -  C(r  +  ft  +  >iff  +  £-»v 

M„ 

a3  =  T7~  +  (£  + 

MII 

^  (C  +  ?/) 

Mu 

To  apply  the  Routh-Horowitz  stability  criterion  first  establish  the  following  deter¬ 
minates  21  . 


i  —  at 


®  2 


al  a0 
a3  a2 


2)  3 


aj  a0  0 

a3  a2  al 

0  a4  a3 


The  necessary  and  sufficient  set  of  conditions  for  stability  are  that  a„  >  0  and  2>  lf 
32i  are  each  positive.  Fig.  V.2  is  a  map  of  the  stable  and  unstable  zones  of 
operation  Mori,  Mori^  calculated  by  the  method  similar  to  that  described  here  and 
actual  test  values.  Their  results  are  for  a  circular  bearing  geometry  such  as  shown  in 
Fig.  IV.2  but  with  an  additional  R-C  component  placed  between  stations  2  and  1.  The 
good  agreement  between  calculated  and  test  values  establishes  the  creditability  of 
their  theory.  The  results  also  showed  that  capillary  restrictors  do  as  well  as  orifice 
restrictors  in  this  stabilizing  scheme  as  long  as  dimensions  are  chosen  from  the  res¬ 
tricted  stabilized  zone  shown. 
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z'  =  (C,_»),  IO"0cm*sec2 

Fig.  V.2  Zones  of  Stable  and  Unstable  Operation  From  Ref.  4. 

D.  Wilcock^  developed  the  dynamics  of  EX.P  thrust  bearings  somewhat  along  the 
lines  above^  using  servomechanism  analytics.  He  includes  more  parameters  such  as 
oscillatory  loading  and  inertance  (inertance  defined  in  Section  I,  Table  A.I.c  Element 
1)  into  the  1-2  restrictor  branch.  No  experimental  confirmation  is  given  however. 

Precision  Levels  of  Aerostatic  Linear  Guideways 

A  Kobayashi^  describes  ultra-precision  machine  tool  developments  in  Japan 
directed  towards  diamond  cutting  associated  with  the  fabrication  of  computer  memory 
disks,  laser  scanner  mirrors  etc..  Hydrostatic  and  aerostatic  guidance  play  significant 
roles  in  these  developments.  Performance  figures  for  linear  guideways  such  as  Fig.  V.3 


are  given  in  the  accompanying  table.  Exact  construction  details  are  not  given  but  a 
plausible  assumption  is  that  the  aerostatic  lands  contain  no  explicit  pockets  but  rather 


'  "l  -t*  ‘V  “ 
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each  land  is  fed  by  a  straight  line  array  of  orifices  essentially  as  described  in  reference 


(3).  The  table  I  indicates  that  the  guidance  fidelity  of  aerostatic  slide  ways  can  be 


made  good  to  an  order  of  20  nanometers!  This  ultra  precision  is  from  a  guideway 
which  can  nevertheless  support  loads  of  the  order  of  2000  N. 


Slid*- 


Hydrostatic  Linear  Slid* 

■TTSmiHi 

icEXJirn 

Naa*  of  Manufacturer 

T0Y0DA 

MACHINE 

NACHI- 

FUJXKOSHI 

HITACHI 

SEIKO 

TOSHIBA 

MACHINE 

Strok*  (mm) 

300 

350 

200 

200 

Stralghtri*** 

(ua/lOOaa) 

0.15 

0.029 

0.05(H) 

0.015(V) 

■HjBjQZTT 

Stiffnesa  (N/ua) 

— 

— 

49b 

294 

Allowable  Load 
Capacity  (N) 

2940 

1940 

H:Horlaontal  direction,  V: Vert l cal  direction 


Fig.  V.3  Aerostatic  Linear  Slide  (Hitachi  Seiko)  Performance.  Ref.  8 


A  High  Precision  Straight  Line  Positioning  System. 

In  section  IV,  Fig.  IV.8,  a  unique  hydrostatic  guideway  system  was  described  in 
which  the  conventional  passive  control  restrictor  1-2  was  replaced  by  a  valve  activated 
by  the  vertical  position  of  the  slide  table  relative  to  a  horizontal  straight  edge.  An 
important  aspect  of  this  concept  is  the  removal  of  the  bearings  film  thickness  h  as  the 
parameter  of  position  control.  This  therefore  removes  the  influence  of  manufactured 
errors  inherent  in  the  guideway.  Instead  the  more  precise  straightedge  becomes  the 
controlling  datum.  This  concept  is  expanded  by  A.  Shimokohbe  etal^  to  create  a  high 
precision  positioning  system  to  meet  the  needs  of  ultra  precision  tooling  such  as  might 
be  used  in  diamond  turning  machines.  Figures  V.4a,b  show  the  system  assembly.  The 
object  is  to  position  the  block  so  that  five  kinematic  errors  are  held  within  fractions  of 
microns  or  arc  seconds.  There  is  a  primary  table  (mechanical  slide  of  Fig.  V.4a)  run¬ 
ning  in  a  guide  with  conventional  thin  film  lubrication.  This  guideway  need  possess 
only  moderate  precision  since  it  provides  merely  nominal  positioning.  Fine  positioning 
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is  super  imposed  on  top  of  the  nominal  position  by  controlling  the  film  thickness  of 
nine  aerostatic  circular  thrust  bearings  Fig.  V.4b  which  support  the  block.  The  table 
below  demonstrates  the  position  errors  first  observed  without  the  feed  back  control 
and  then  with  the  feedback  controls  activated.  The  feedback  system  registers  approx¬ 
imately  a  30  fold  improvement  in  position  precision.  The  complexity  of  feedback  posi¬ 
tion  control  is  obviously  acceptable  when  high  precision  is  the  object. 


TABLE:  Measured  Position  Errors 


Error  Type 

Vertical 

Horizontal 

Roll 

Pitch 

Yaw 

Free  System 

8/im 

10/im 

10  sec 

10  sec 

10  sec 

Fe  .»d-Back  System 

0.1/xm 

0.4/im 

0.2  sec 

0.4  sec 

0.2  sec 

The  system  operates  as  follows.  A  pair  of  precision  straight  edges,  Fig.  V3.a,  is 
probed  by  a  pair  of  capacitance  micrometers  whose  readings  will  be  converted  to  sig¬ 
nals  registering  the  components  of  vertical  and  roll  errors  of  the  block  relative  to  the 
straight  edges.  An  autocollimator  registers  the  pitching  error  of  the  block  position. 
These  data  activate  three  servo  values.  Each  of  these  valves  controls  the  gas  flow  rate 
to  one  of  the  three  aerostatic  pads  bearing  on  the  underside  of  the  block.  The  princi¬ 
ple  involved  is  similar  to  that  described  for  the  adaptive  control  of  the  hydrostatic 
bearing  of  Fig.  IV.5.  In  that  case  the  probe  data  being  mechanical,  was  converted 
directly  t<>  a  valve  opening.  Thus  time  the  probe  signals  being  electrical  they  must  be 
converted  to  mechanical  valve  opening  by  solenoids.  While  the  electrical  feedback 
involves  a  more  complex  system  it  has  an  advantage  over  the  mechanical  probe-valve 
controllei  of  section  IV.  That  is,  computer  modification  of  system  characteristics  can 
be  readily  imposed  adding  to  versatility  of  performance. 


The  three  aerostatic  pads  acting  on  the  under  side  of  the  block  tilt  and  lift  the 
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block  to  negate  the  errors  reported  by  the  capacitance  probes  and  the  autocollimator. 
A  second  autocollimator  beams  across  the  top  of  the  block  to  monitor  the  yaw  errors. 
Its  signal  activates  the  servo  valves  associated  with  pads  ap  b^  a2,  b2.  The  laser  sys¬ 
tem  of  Fi"j.  V.4a  beams  onto  a  photodetector  to  register  horizontal  shift  errors.  This 
error  also  influences  the  servo  valves  associated  with  pads  a1(  b,,  a2,  b2.  The  combina¬ 
tion  of  yaw  and  horizontal  shift  signals  causes  pads  aj,  b1?  a2,  b2  to  expand  or  con¬ 
tract  to  negate  these  errors. 

The  two  pads  oriented  along  the  axis  of  table  travel  act  as  a  pneumatic  vise. 
Presumably  if  the  outward  travel  of  the  block  were  measured  by  an  interferometer 
these  pads  could  be  employed  to  add  the  sixth  dimension  of  control. 

The  details  of  the  air  pads  are  also  informative.  Fig.  V.5  shows  that  an  annular 
channel  :s  employed  obviously  to  reduce  the  pneumatic  capacitance  of  the  pad  in  the 
interests  of  removing  susceptibility  to  air  hammer.  The  channel  is  unusually  shallow 
suggesting  great  care  was  taken  to  develop  the  absolute  minimum  capacitance.  This 
observation  implies  also  that  feedback  control  was  not  relied  upon  to  suppress  air 
hammer.  The  plot  of  Fig.  V.5  also  offers  useful  data  regarding  practical  magnitudes  of 
film  thick-  ness  and  pressures.  The  real  character  of  the  static  load-clearance-  supply 
pressure  relationships  might  very  well  be  compared  with  analytical  results. 

The  example  of  the  adaptive  controlled  hydrostatic  slideway  Fig.  IV.5  shows  a  no- 
through-flow  bearing.  The  side  seals  reduced  the  oil  losses  which  with  the  high  pres¬ 
sures  involved  in  that  example  would  translate  both  to  extensive  contamination  of  the 
machine  area  reduced  effectiveness  of  the  adaptive  system  by  reintroducing  the 
influence  of  pocket  geometry  errors.  The  no-through-flow  concept  was  not  used  for  the 
system  of  Fig.  V.4a,b.  The  leakage  of  air  out  of  aerostatic  bearings  does  not  equate  to 
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contamination  of  the  metrology  environment  (except  for  causing  the  laser  beam  to 
wander  if  it  is  not  shielded  from  air  currents). 

Perimeter  seals  needed  for  no-through-flow  pads  would  add  complexity  without 
profit. 


Fig.  V.4b  Pads,  Servo  Valves  and  Block  Assembly.  Ref.  7. 
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ClMrincc,  jim 

Supply  pressure: 

(— C— )  0.098 MPa.  ( — S — )  0.196  MPa.  (—€>—) 

0.294  MPa.  (— ®-)  0.392  MPa.  (— •-) 

0490  MPa.  <-•-)  0.588  MPa.  (-3-) 

0.686  MPa 

Fig.  V.5  Load,  Clearance  and  Supply  Pressure  of  Pads.  Ref.  7. 
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APPENDIX:  AV 


Derivation  of  — - —  and  — - —  where  M  2  3  *s  mass  flow  through  a  thin  plate 

^P3/4  C*P2/* 

orifice.  Ipstream  pressure  is  p2/*,  downstream  pressure  is  p3/«.  The  flow  through  an 
orifice  is  customarily  described  as 


V 


A V  tk  (x=r)l(^)T  -  C* 1 : 

where:  9t  is  the  gas  constant  for  the  fluid  involved 


E.AV.l 


T2  is  the  temperature  at  state  2.  9ft  T2 


'1  is  the  polytropic  constant  —  1.4  (air) 


Let:  *  _  ft  Cd  =  coefficient  of  discharge  (assumed  invariant  in  this  analysis). 

P2/  * 

/Td„ 

Aq  =  orifice  area  =  - or  periphery  x  h  for  proximity  orifice. 

4 


a)  Evaluating  — 


otff  2  3 


Rewrite  E.AV.l  as 


(t  23  =  B.X  B  =  CdA 


dAoP2/‘ 


V 


—(-7) 

P2/*  l-1 


E.AV.2b 


X  =  ft'  -  ii 


Then 


<)M  2  3 _ B_  (IX 

< ;P3 /•  2X'"  tVp3/. 


-  154  - 


AV  3 


=  Ps/*  (-Pi/* )  =  —~ 
Op2/*  '  '  '  P2  /’ 


Therefore 


-/f /-  -(l+1)^1/'+1  — 
i  7  7  P2 /* 


t^p2/*  <^P2/* 


CdAoP2/ 


V 


Pi  7  _  B 

P2/*  7—1  2P2/’ 


Finally 


(lH  2  3  —M  2.3 


9  2-3  + 

P2/*  2p2  /• 


fill  2  3  *  1 1 


2.3  v"  2.3  *  r  _  I  <2f 

Of  “77 -  =  "2—3  >  "2—3  “  „  'r  2-3 

7p2/*  P2/‘  l  Z 


E.AV.4 


C<W?  53  t  5 

Derivation  of  — —  and  — — 
- <9p3/.  dh 


% 


a  *'<*  h  _L  n, 

a  5.3  =  ~  •  o  •  ^5/3 

Poo/*  14//  or 


DM  c-x  M  ,s 


w  3  5  *3/5 

*^3/5  r,p3/* 


()p3/.  7p3)/-* 


(P3V  "  P52/*)  =  2Pa/‘ 


( LiX  53 

Hence  — - —  =  —  H 

<>?3  r 


2p3/‘ 

'M  3.5 

2 

%/s 

p3/’ 

p3/’ 

■mwmwMmmmmmwmm. 
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